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Coulomb drag is a transport phenomenon whereby long-range Coulomb interaction be¬ 
tween charge carriers in two closely spaced but electrically isolated conductors induces 
a voltage (or, in a closed circuit, a current) in one of the conductors when an electrical 
current is passed through the other. The magnitude of the effect depends on the exact 
nature of the charge carriers and microscopic, many-body structure of the electronic 
systems in the two conductors. Drag measurements have become part of the standard 
toolbox in condensed matter physics that can be used to study fundamental properties of 
diverse physical systems including semiconductor heterostructures, graphene, quantum 
wires, quantum dots, and optical cavities. 
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I. FRICTIONAL DRAG 

Inner workings of solids are often studied with the 
help of transport measurements. Within linear response, 
the outcome of such measurements is determined by the 
properties of the unperturbed system, which are often 
the object of study. In a typical experiment a current is 
driven through a conductor and the voltage drop along 
the conductor is measured. In conventional conductors 
at low temperatures the resulting Ohmic resistance is 
mostly determined by disorder (which is always present in 
any sample) (Lifshitz and Pitaevskii, 1981; Ziman, 1965), 
while interactions between charge carriers lead to correc¬ 
tions that affect the temperature dependence of transport 
coefficients (Altshuler and Aronov, 1985). 

In his pioneering work, Pogrebinskii (1977) has sug¬ 
gested an alternative measurement that involves two 
closely spaced, but electrically isolated conductors (here¬ 
after referred to as “layers”). In such a system, an elec¬ 
tric current Ii flowing through one of the layers, known 
as the “active” layer, induces a current (or, in an open 
circuit, a voltage V 2 , see Fig. 1) in the other, “passive” 
layer by means of “mutual friction”. By this one typically 
understands scattering between charge carriers belonging 
to different layers due to long-range interactions. These 
scattering events are accompanied by energy and momen¬ 
tum transfer from the carriers in the active layer to the 
carriers in the passive layer, effectively “dragging” them 
along. At the simplest level, such friction effects can be 
described by introducing a phenomenological relaxation 
rate. In the case of frictional drag, the corresponding rate 
generally depends on the exact nature of the charge 
carriers, interlayer interaction, and microscopic structure 
of the electronic system. Thus, measurements of this re¬ 
laxation rate provide additional insight into microscopic 
properties of interacting many-body systems. 

A related phenomenon, where a quasiparticle flow in¬ 
stigates a partial transfer of energy and momentum be¬ 
tween separate, but interacting subsystems of quasiparti¬ 
cles, is known as “phonon drag” (Gurevich, 1946a,b; Her¬ 
ring, 1954) and manifests itself in a rising thermoelectric 
power in semiconductors at low temperatures (Fredrikse, 
1953a,b; Geballe, 1953; Geballe and Hull, 1954). In the 
presence of a temperature gradient, lattice vibrations be¬ 
come anisotropic since the phonons travel preferentially 
from hot to cold (providing a mechanism for thermal con¬ 
duction). Interacting with electrons, the phonons effec¬ 
tively drag them towards the cold end of the sample, 
creating an excess charge density (this process will con¬ 
tinue until the electrostatic field created by the accumu- 



FIG. 1 Schematic showing the drag signal V 2 induced by the 
current 7i. [From Price et al. (2007). Reprinted with permis¬ 
sion from AAAS.j 


lated charge will counterbalance the drag effect). In a 
nonequilibrium system of electrons and phonons, their 
mutual drag is intertwined with heating effects and af¬ 
fects charge transport (Gurevich and Mashkevich, 1989). 
The resulting correction to the standard transport theory 
is important in thermoelectric measurements. 

In contrast, frictional drag in double-layer systems is 
not a correction: in the absence of the interlayer inter¬ 
action charge carriers in two disjoined conductors are in¬ 
sensitive to each other (therefore, any drag effect should 
necessarily vanish in the limit of infinitely remote lay¬ 
ers). In other words, the drag phenomenon simply does 
not exist in noninteraeting systemsl Gonsequently, ini¬ 
tial experimental work on mutual drag was devoted to 
quantitative measurement of the strength of interactions 
between quasiparticle subsystems in various semiconduc¬ 
tor devices including p-modulation-doped GaAs quan¬ 
tum wells (Hopfel and Shah, 1988; Hdpfel et ai, 1986), 
capacitively coupled two- and three-dimensional (2D- 
3D) electron systems in AlGaAs/GaAs heterostructures 
(Solomon and Laikhtman, 1991; Solomon et al., 1989), 
2D electron systems in AlGaAs/GaAs double quantum 
wells (Eisenstein, 1992; Gramila et al., 1991, 1992, 1994; 
Solomon and Laikhtman, 1991), and electron-hole bilay¬ 
ers (Sivan et al., 1992). Drag between 3D systems was 
numerically simulated in Jacoboni and Price (1988). At 
low temperatures and for closely spaced layers, the inter¬ 
layer scattering rate appeared to be dominated by 
the Coulomb interaction (Price, 1983, 1988). 

Coulomb drag between spatially separated electron 
systems is ultimately caused by fluctuations (or inhomo¬ 
geneities) of the charge density in the two layers (Zheng 
and MacDonald, 1993). Indeed, an infinite layer with the 
uniformly distributed electric charge creates a uniform 
electric field in the normal direction that does not exert 
any lateral force upon the carriers in another layer. If 
both layers are in the Fermi-liquid state, then the usual 
phase-space argument (Gramila et al., 1991) yields the 
quadratic temperature dependence oc in qualita- 
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tive agreement with the observed behavior at low enough 
temperatures. More detailed analysis of the experimen¬ 
tal data revealed the presence of additional mechanisms 
leading to frictional drag, such as the indirect interlayer 
interaction mediated by phonons (Gramila et al., 1993; 
Noh et al., 1999; Rubel et al., 1995, 1996), plasmon effects 
(Hill et al., 1997; Noh et al., 1998), and thermoelectric 
phenomena (Solomon and Laikhtman, 1991). 

Theoretically, it was realized early on that mutual 
Coulomb scattering between electrons in the two layers 
results in the exchange of both energy and momentum 
(Boiko and Sirenko, 1988; Maslov, 1992; Price, 1983). 
Initial calculations aimed at energy and momentum re¬ 
laxation in a nondegenerate 2D electron gas (2DEG) 
due to proximity to a 3D conductor (Boiko and Sirenko, 
1988, 1990) were followed by investigation of transport 
properties in coupled 2D and 3D systems (Boiko et al., 
1992; Laikhtman and Solomon, 1990), ID systems cou¬ 
pled to conductors of arbitrary dimensionality (Sirenko 
and Vasilopoulos, 1992), coupled ID wires (Gurevich 
et al., 1998; Tanatar, 1996, 1998; Tso and Vasilopoulos, 
1992), and quantum Hall edge states (Orgad and Levit, 
1996). Following the groundbreaking experiments in 
AlGaAs/GaAs double quantum wells (Eisenstein, 1992; 
Gramila et al., 1991), a lot of work was devoted to 
drag between two degenerate 2DEGs. While the purely 
Coulomb mechanism (Flensberg et al., 1995; Jauho and 
Smith, 1993; Kamenev and Oreg, 1995; Zheng and Mac¬ 
Donald, 1993) does capture the most qualitative features 
of the effect, other mechanisms of momentum transfer 
may also contribute to the observed behavior. In samples 
with larger interlayer spacing {d ~ 50-500nm) as much 
as 30% of the measured signal was attributed to phonon- 
mediated interactions (Gramila et al., 1993). These mea¬ 
surements appeared to be consistent with the virtual- 
phonon exchange mechanism (Tso et al., 1992, 1994). 
Other suggested scattering mechanisms involved acoustic 
(Bpnsager et al., 1998a) and optical (Hu, 1998) phonons, 
plasmon effects (Flensberg and Hu, 1994) and coupled 
plasmon-phonon modes (Giiven and Tanatar, 1997a). Bi¬ 
layers subject to strong magnetic field were shown to 
form interlayer correlated states (Girvin and MacDonald, 
1997; Varma et al, 1994). For superconducting layers, in¬ 
terlayer magnetic interaction due to spontaneously cre¬ 
ated vortices has also been suggested (Shimshoni, 1995). 

Mutual Goulomb scattering has been studied also in 
a hybrid device (Huang et al., 1995) comprising normal 
(Au/Ti) and superconducting (AlOa,) 2D films separated 
by an insulating (AI 2 O 3 ) layer. In that case, as well as in 
“cross-talk” measurements in superconductor-insulator- 
normal-metal trilayers (Giordano and Monnier, 1994), 
the phenomenological Drude-like description of drag in 
terms of does not apply. The Drude description also 
fails when the system is subjected to a strong magnetic 
field: in contrast to the naive description, numerous ex¬ 
periments (Feng et al., 1998; Hill et al., 1996, 1998; Jorger 


et al., 2000c; Lok et al., 2001a,b; Patel et al., 1997; Rubel 
et al., 1997a,b) show significant dependence of the mea¬ 
sured drag resistivity pd on. the applied field, especially 
in the extreme quantum regime (Lilly et al., 1998). More 
sophisticated theoretical calculations on Coulomb drag 
in quantum Hall states (Shimshoni and Sondhi, 1994), 
snperflnid condensates in paired electron-hole layers (Vi- 
gnale and MacDonald, 1996), drag of composite fermions 
(Kim and Millis, 1999; Ussishkin and Stern, 1997, 1998; 
Zhou and Kim, 1999), vortex drag (Vitkalov, 1998), non- 
dissipative drag (Rojo and Mahan, 1992), supercurrent 
drag (Duan and Yip, 1993), as well as drag between 
charged Bose gases (Tanatar and Das, 1996) and meso¬ 
scopic rings (Baker et al., 1999; Shahbazyan and Ulloa, 
1997a,b) have confirmed the expectation that the drag re¬ 
sistivity reflects not only the exact character of interlayer 
interaction, but also the nature of elementary excitations 
in each layer and their fundamental properties. 

After the turn of the century, drag measurements be¬ 
came part of the standard toolbox in condensed matter 
physics. They have been used to investigate properties 
of electron-electron scattering in low-density 2D electron 
systems (An et al., 2006; Kellogg et al., 2002a); signatures 
of metal-insulator transition in dilute 2D hole systems 
(Jorger et al., 2000a,b; Pillarisetty et al., 2002, 2005a,b); 
quantum coherence of electrons (Kim et al, 2011; Price 
et al., 2008, 2007) and composite fermions (Price et al., 
2010); exciton effects in electron-hole bilayers (Groxall 
et al., 2008; Keogh et al., 2005; Morath et al., 2009; Sea- 
mons et al., 2009); exotic bilayer collective states (Eisen¬ 
stein, 2014), especially the quantum Hall effect (QHE) 
at the total filling factor = 1 (Finck et al., 2010; Kel¬ 
logg et al., 2003, 2002b; Schmnlt et al, 2010; Spielman 
et al., 2004; Tutuc et al., 2009); compressible quantum 
Hall (QH) states at half-integer filling factor (Muraki 
et al., 2004; Zelakiewicz et al., 2000); integer QH regime 
(Lok et al., 2002); Luttinger liquid effects (Debray et al., 
2001; Laroche et al., 2008, 2014); Wigner crystallization 
in quantum wires (Yamamoto et al., 2002, 2006, 2012); 
and one-dimensional (ID) sub-bands in quasi ID wires 
(Debray et al., 2000; Laroche et al., 2011). More gener¬ 
ally, interlayer interaction and corresponding transport 
properties have been studied in hybrid devices compris¬ 
ing a quantum wire and a quantum dot (Krishnaswamy 
et al., 1999); a SG film and a 2D electron gas (Fa¬ 
rina et al., 2004); Si metal-oxide-semiconductor systems 
(Laikhtman and Solomon, 2005); quantum point contacts 
(Khrapai et al., 2007); insulating a-SiNb films (Elsayad 
et al., 2008); ferromagnetic-antiferromagnetic-SC trilay¬ 
ers (Guoco et al., 2009); nanosize GdSe-CdS semiconduc¬ 
tor tetrapods (Mauser et al., 2010); electron-hole scatter¬ 
ing in quantum wells (Prunnila et al., 2008; Takashina 
et al., 2009; Yang et al., 2011); graphene monolayers 
(Gorbachev et al., 2012; Kim et al, 2011; Kim and Tu¬ 
tuc, 2012; Titov et al, 2013a); and hybrid graphene- 
semiconductor systems (Gamucci et al., 2014). 



4 


On the theory side, the variety of suggested extensions 
and generalizations of the original drag problem is even 
richer. The theory of Coulomb drag between two 2DEGs 
was extended to dilute 2D hole systems (Hwang et ai, 
2003) and to the cases where one allows for certain tun¬ 
neling processes between the layers (Oreg and Halperin, 
1999; Oreg and Kamenev, 1998), interlayer disorder cor¬ 
relations (Gornyi et al., 1999; Hu, 2000a), in-plane po¬ 
tential modulation (Alkauskas et al., 2002), and disor¬ 
der inhomogeneities (Apalkov and Raikh, 2005; Spivak 
and Kivelson, 2005; Zou et ai, 2009, 2010). Theory of 
Coulomb drag between composite fermions was gener¬ 
alized to include phonon-mediated coupling (Bpnsager 
et ai, 2000; Khveshchenko, 2000). Mutual friction was 
also suggested to occur between non-Fermi-Liquid phases 
including Luttinger liquids (Flensberg, 1998; Klesse and 
Stern, 2000; Nazarov and Averin, 1998), Wigner crys¬ 
tals (Baker and Rojo, 2001; Braude and Stern, 2001), 
and strongly localized electrons (Raikh and von Oppen, 
2002). Drag or similar measurements of interlayer in¬ 
teractions were also considered for composite (or hy¬ 
brid) systems comprising ballistic quantum wires (Gure¬ 
vich and Muradov, 2000, 2005; Raichev and Vasilopou- 
los, 2000; Wang et al, 2005), coupled 2D-1D systems 
(Lyo, 2003), nonequilibrium charged gases (Wang and 
da Cunha Lima, 2001), multi-wall nanotubes (Lunde 
et ai, 2005; Lunde and Jauho, 2004), quantum point 
contacts (Levchenko and Kamenev, 2008a), few level 
quantum dots (Moldoveanu and Tanatar, 2009), opti¬ 
cal cavities (Berman et al, 2010a, 2014), coupled meso¬ 
scopic rings (Yang and MacDonald, 2001), supercon¬ 
ductors (Levchenko and Norman, 2011), and normal- 
metal-ferromagnet-normal-metal structures (Zhang and 
Zhang, 2012). Other developments include mesoscopic 
fluctuations of Goulomb drag (Narozhny and Aleiner, 
2000; Narozhny et al, 2001), frictional drag mediated 
by virtual photons (Donarini et al, 2003) and plasmons 
(Badalyan et al, 2007), exciton effects in semiconductors 
(Laikhtman and Solomon, 2006) and topological insula¬ 
tors (Mink et al, 2012), interlayer Seebeck effect (Lung 
and Marinescu, 2011) and spin drag (Badalyan and Vi- 
gnale, 2009; D’Amico and Vignale, 2000; Duine et al, 
2011, 2010; Duine and Stoof, 2009; Glazov et al, 2011; 
Pustilnik et al, 2006; Tse and Das Sarma, 2007; Vig¬ 
nale, 2005). Recently, the focus of the theoretical work 
was shifted towards the drag effect in graphene-based 
devices (Narozhny, 2007; Narozhny et al, 2015; Song 
et al, 2013; Tse et al, 2007) and strongly interacting 
high-mobility double-layers with low-density carrier con¬ 
centration (Apostolov et al, 2014; Chen et al, 2015). 

Given the rather large amount of literature devoted to 
frictional drag, it seems unreasonable to cover all pos¬ 
sible angles in a single paper. Early work on frictional 
drag was reviewed by Rojo (1999). Various experimen¬ 
tal aspects were discussed in reviews on exciton conden¬ 
sates (Eisenstein, 2014; Snoke, 2002), electron-hole bilay¬ 


ers (Das Gupta et al, 2011), strongly-correlated 2D elec¬ 
tron fluids (Spivak et al, 2010), and ID ballistic electron 
systems (Debray et al, 2002). A discussion of drag in 
strong magnetic fields was included in a review of mag- 
netotransport in 2D electron systems (Dmitriev et al, 
2008). In the present review, we limit ourselves to the dis¬ 
cussion of standard (“electrical”) Coulomb drag. Spin- 
related phenomena and thermoelectric effects are beyond 
the scope of this review. 


II. COULOMB DRAG IN SEMICONDUCTOR 
HETEROSTRUCTURES 


In an idealized experiment, a constant (dc) current Ii is 
passed through the active layer, keeping the passive layer 
isolated at the same time (such that no current is allowed 
to flow in it), see Fig. 1. The voltage V 2 induced in the 
passive layer is proportional to Ii and the coefficient^ 

Rd = -V 2 /I 1 , (1) 


is a direct measure of interlayer interactions. 

In his original paper, Pogrebinskii (1977) derived the 
Drude-like description of transport in double-layer sys¬ 
tems comprising two coupled equations of motion 


dvi e e , 

— = — El H-Vi X B] 

dt mi mic 

Vl 

Vl - V2 

, (2a) 

n 

td 

dv 2 e ^ e , 

= — E 2 H- [V2 X B] 

dt m 2 m2C 

V2 

V2 - Vl 

, (2b) 

T2 

Td 

where e is the electric charge. 

Vi, mi, and Ei 

are the 


drift velocities, effective masses, and electric fields in the 
two layers, and the nonquantizing magnetic field B is 
assumed to be uniform. Intralayer impurity scattering 
processes yielding the Drude resistivity in the two layers 
are described by the mean free times Ti. The last term 
in each of Eqs. (2) describes the mutual friction between 
the charge carriers in the two layers that tends to equalize 
drift velocities. If treated phenomenologically, the model 
(2) describes two distinct types of carriers coupled by the 
friction term, but does not explicitly require them to be 
spatially separated (Gui et al, 1988; Hansch and Mahan, 
1983; Soderstrom et al, 1996). 

Solving the equations (2), one finds the resistivity ma¬ 
trix [hereafter the indices i,j = 1,2 denote the two 
layers and a, j5 = x,y - spatial coordinates orthogonal to 
B = BEz] the “layers” described by Eqs. (2) can repre¬ 
sent 2D or 3D conductors, see Sec. VI for the ID case]. 


^ The minus sign in Eq. (1) is motivated by Eq. (3a). An alterna- 
tive definition without the explicit minus sign is also widely used 
in literature. 
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The “drag resistivity” (also called the transresistivity or 
the drag coefficient) is given by the Drude-like formula 

PD = -P^xx'’ = (3a) 


The expression (3a) is independent of the magnetic field. 
This statement has the same status as the absence of 
the classical magnetoresistance^. Indeed, the single-layer 
longitudinal resistivity derived from Eqs. (2) is given by 


P 


( 11 ) = 


nil 

e^ni 




(3b) 


In most cases, drag is rather weak (td ^ r) and the usual 
Drude formula remains a good approximation for p^x'’ 
(Eisenstein, 1992; Rojo, 1999). The single-layer Hall co¬ 
efficient is unaffected by the presence of the second layer 
and is determined solely by the carrier density 

= Bfiniec). (3c) 

Within the applicability of the Drude model, frictional 
drag is purely longitudinal: “Hall drag” does not occur^ 

Pd = Pyx^ = 0- (3d) 

At the phenomenological level, the drag resistivity (3a) 
is independent of the disorder strength. Moreover, in the 
“clean” limit r —>■ oo the inter- and intralayer resistivities 
tend to the same value and the resistivity matrix becomes 
degenerate (the corresponding conductivities diverge): 

Pxx\t ^ = Pd{j ^ og). (4) 


Thus a system comprising two capacitively coupled, ideal 
conductors is characterized by non-zero resistivity and 
exhibits perfect drag! 


A. Interlayer Coulomb interaction 

The “Drude formula” (3a) for the drag resistivity be¬ 
comes falsifiable provided that something is known about 
the properties of the “drag rate” (i.e. its dependence 
on temperature, carrier density, interlayer separation, 
and other experimentally relevant parameters). To the 
leading order, the contribution of the interlayer Coulomb 
interaction to can be calculated within the Born ap¬ 
proximation (or equivalently, using the Fermi’s Golden 
Rule) (Jauho and Smith, 1993; Laikhtman and Solomon, 





FIG. 2 (Color online) Aslamazov-Larkin diagrams describing 
the lowest-order contribution to drag. The solid lines refer to 
quasiparticle Green’s functions and the wavy lines describe 
the interlayer interaction. The left and right triangles corre¬ 
spond to non-linear susceptibilities of the two layers. 


1990). In the language of Feynman diagrams, the corre¬ 
sponding process (Flensberg et al, 1995; Kamenev and 
Oreg, 1995; Zheng and MacDonald, 1993) is described by 
the Aslamazov-Larkin diagrams (Aslamazov and Larkin, 
1968) shown in Fig. 2. 

The effective interlayer interaction can be found as a 
solution to the Poisson equation for the potential of a 
point source belonging to one of the layers. In principle, 
this can be done for any system of coupled conductors. 
Coupling between a 2DEG and a 3DEG was considered in 
Laikhtman and Solomon (1990). A double quantum well 
system was discussed in Jauho and Smith (1993) where 
the finite width of the wells was taken into account by 
assuming a specific form of the electron wave function 
in the direction perpendicular to the layers. However, 
the obtained results are qualitatively the same as in the 
simplest case of purely two-dimensional layers. 

If electrons in each layer are confined to move in a 2D 
plane, the “bare” Goulomb potential has the form"* 

Vii = V 22 = ^-Ke^jq, Vi2{q) = j q)e~'^^. (5) 

Here e is the electron charge and d is the interlayer sep¬ 
aration, that determines the maximum value (or rather, 
the order of magnitude thereof) of the momentum q that 
can be transferred between the layers^: 

q < 1/d. (6) 


Taking into account dynamical screening within the usual 
Random Phase Approximation (RPA) modifies the inter¬ 
layer interaction propagator (Stern and Halperin, 1995), 
but does not change the exponential decay at large q: 


12 nfn|«^sinhgd-b[2:;^-bHf-bnf]e 9 rf' 

Here Hf is the single-layer retarded polarization opera- 
tor. It is quite common [see, e.g., Jauho and Smith (1993) 


^ In this Section we are discussing the simplest situation, where 
both Tjf) and r are unaffected by weak enough magnetic fields. 

^ Under the assumptions of the present Section, magnetic field 
has no effect on drag. Hence, up until Sec. II.G we focus on 
the zero-field, longitudinal transport. Drag in magnetic field is 
discussed in Secs. II.G, IV.E, and VILA. 


Although electrons are confined to move in two dimensions, they 
interact by means of “real, 3D” Coulomb interaction. 

^ While discussing the theory, we use the natural units with 
h = ks = !• We attempt to restore the Planck’s constant in final 
answers and while discussing experimental findings. 
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and Laikhtman and Solomon (1990)], to include the di¬ 
electric constant e of the insulating spacer into the “bare” 
potential. Since the same e should enter the expression 
for the inverse Thomas-Fermi screening length 

H = 2Tie^v = 27re^n^((7 < 2kF,oj = T = 0), (8) 

the dielectric environment can be taken into account by 
expressing the results in terms of k {v denotes the ther¬ 
modynamic density of states of the 2DEG). For high car¬ 
rier densities (Gramila et al., 1991), Eq. (7) can be sim¬ 
plified (Kamenev and Oreg, 1995) by assuming the small 
screening length, 1 [see Eqs. (20) and (33b) below]. 

The condition (6) allows one to distinguish the follow¬ 
ing two regimes (Kamenev and Oreg, 1995): 

(i) if the interlayer separation is large compared to the 
mean-free path £, then it follows from Eq. (6) 
that (7 <C l/£; in this case the motion of charge car¬ 
riers is diffusive] 

(ii) in the opposite case, d <^ £, transport is domi¬ 
nated by ballistic propagation of charge carriers 
with l/d^q^l/£, see Eq. (36) below. Most 
measurements (Gorbachev et at, 2012; Gramila 
et al., 1991) are performed on ballistic samples. 

The majority of analytic (Rojo, 1999) and numerical 
(Mosko et al., 1992) work on Goulomb drag in semicon¬ 
ductor heterostructures was performed treating the in¬ 
teraction (7) in the lowest order of perturbation theory. 
For generalizations see Secs. II.D and ILF. 


equilibrium Fermi-Dirac distribution function (as ap¬ 
propriate within linear response) are described by 

[i - /f]/i.=-r[a/fV5e]/i,. (lo) 

Here we only consider the steady state and uniform fields 

dfjdt = 0, V/, = 0. (11) 

The latter condition physically means that the sample 
size is large compared to the length scale of typical re¬ 
laxation processes in the system, see also Sec. IV.D. 

In the absence of interlayer interaction, the task of find¬ 
ing linear-response transport coefficients from Eq. (9) is a 
textbook problem (Seeger, 2002; Smith and Jensen, 1989; 
Ziman, 1965). Under the above assumptions, the theory 
is qualitatively equivalent to the Drude theory (2) yield¬ 
ing the standard results (3b) and (3c). Not surprisingly, 
taking into account the collision integral lij leads to the 
Drude-like description of the drag resistivity (3a) and 
(3d). The advantage of the present “microscopic” calcu¬ 
lation is that now we can determine the phenomenologi¬ 
cal relaxation time to in terms of the model parameters. 

The standard perturbative calculation (Boiko et al., 
1992; Jauho and Smith, 1993; Laikhtman and Solomon, 
1990; Lifshitz and Pitaevskii, 1981) amounts to finding 
the nonequilibrium distribution functions hi in the two 
layers to the leading order in the interlayer interaction 
and the electric field Ei applied to the active layer. Then 
one uses the definition of the electric current (here the 
sum runs over all of the single-particle states) 


B. Kinetic theory of ballistic drag 

Ballistic motion of charge carriers in semiconductors 
can be described by using the kinetic equation approach, 
where impurity scattering is taken into account within 
the simplest r-approximation (Jauho and Smith, 1993; 
Laikhtman and Solomon, 1990; Pogrebinskii, 1977). One 
starts with the generic Boltzmann equation 

^+v,Vf.+ (eE, + ^ K X B]) ^ (9) 

where fi is the distribution function (in layer i = 1,2), 
Tij is the collision integral due to interlayer Goulomb in¬ 
teraction, T is the transport impurity scattering time, 
and Sfi is the nonequilibrium correction to the distribu¬ 
tion function. Here we will only consider degenerate elec¬ 
tron systems [as realized in semiconductor heterostruc¬ 
tures (Gramila et al., 1991)]. Weak deviations from the 


3^ = e^v5fi, ( 12 ) 


and finds the current J2 in the passive layers. The coeffi¬ 
cient of proportionality between j 2 x and Ei^ defines the 
drag conductivity ao- The drag coefficient po can then 
be obtained by inverting the 2x2 conductivity matrix^ 


(To ero 

Pd ^ ~ ) 

(T1CT2 — cr cf\U2 


(13) 


where ai is the longitudinal conductivity in layer i; the 
latter relation follows from the smallness of the effect 


ctd < cn, (14) 

as observed in experiment (Eisenstein, 1992; Rojo, 1999). 
This way, one finds for the phenomenological “drag rate” 
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The nonlinear susceptibility (also known as the rectifica¬ 
tion function) Ti{u!,q) (in layer i) is a response function 
relating a voltage y(ri)e*“* to a dc current it induces by 
the quadratic response: 

J = J drij dr2T{uj;ri,r2)V{ri)V{r2), (16) 

with J being the induced dc current. From gauge invari¬ 
ance fdri r(w) = fdr2r( uj) = 0. 

The same result follows from the standard Kubo for¬ 
mula approach within the diagrammatic perturbation 
theory (Flensberg et al., 1995; Kamenev and Oreg, 1995), 
memory function formalism (Zheng and MacDonald, 
1993), and more general Boltzmann-Langevin theory of 
stochastic kinetic equation (Chen et al., 2015). 


1 . Electron-hole asymmetry and rectification 


The rectification function r(a;, q) is the central object 
in the perturbative theory of Coulomb drag. The ex¬ 
pression (15) of the interlayer relaxation rate in terms of 
r(a;,q) explicitly demonstrates the key role of electron- 
hole asymmetry in the leading-order drag effect®. 

Indeed, in order to induce a voltage (or generate a cur¬ 
rent) in the passive layer, one needs to somehow move the 
charge carriers. This is achieved by transferring momen¬ 
tum from the active layer. The macroscopic state of the 
electronic system in the active layer is characterized by 
the finite electric current driven by an external source. 
In a typical electron gas, there are two kinds of excita¬ 
tions - “electron-like”, with energies e > Ep above the 
Fermi energy (i.e. the occupied states outside the Fermi 
surface), and “hole-like”, with e < Ep. These quasipar¬ 
ticles are oppositely charged. As the current is driven 
through the active layer they move in opposite directions, 
see Fig. 3. Then the active layer can be characterized by a 
nonzero total momentum only if there is some asymmetry 
between electron-like and hole-like quasiparticles. Like¬ 
wise, in the passive layer the momentum is transferred 
equally to electrons and holes, such that the resulting 
state can carry current only in the case of electron-hole 
asymmetry. In conventional semiconductors (Kamenev 
and Oreg, 1995), the electron-hole asymmetry appears 
due to curvature of the conduction band spectrum [lead¬ 
ing to the energy dependence of the density of states 
(DoS) and/or diffusion coefficient]. Consequently, in the 
Fermi-liquid theory the electron-hole asymmetry can be 
expressed (Narozhny et al., 2001) as a derivative of the 


® Another known effect of the electron-hole asymmetry in elec¬ 
tronic systems is the thermopower described by the Mott formula 
(Lunde et al., 2006, 2007; Mott and Jones, 1936). 



FIG. 3 Schematic illustration of the momentum transfer due 
to interlayer interaction. As the current Ji is driven through 
the active layer, electrons and holes are moving in the oppo¬ 
site directions since they carry the opposite charge. Such a 
state has non-zero total momentum only due to electron-hole 
asymmetry. Once the momentum is transferred to the passive 
layer, the electrons and holes there are pushed in the same di¬ 
rection. This process can induce a voltage again only due to 
electron-hole asymmetry. 


single-layer conductivity cri(2) with respect to the chem¬ 
ical potential (assuming either a constant impurity scat¬ 
tering time or diffusive transport). The simple estimate 
9 tTi(2)/9/J. o'i(2)//r then explains the typical smallness 

of the effect (Gramila et al., 1991; Sivan et al., 1992; 
Solomon et al., 1989), see Eq. (14). 

Same arguments can be applied to any system contain¬ 
ing carriers with opposite signs of the electric charge. For 
instance, one can consider semimetals (or even band in¬ 
sulators at high enough temperature) where the electric 
current can be carried by electrons from the conduction 
band and holes from the valence band. A particularly 
interesting example is graphene (see Sec. IV), which ex¬ 
hibits exact particle-hole symmetry at the charge neu¬ 
trality point (Katsnelson, 2012). At that point, the non¬ 
linear susceptibility of graphene (73) vanishes (Narozhny, 
2007; Tse et al., 2007) implying the absence of the drag 
effect. In contrast, experiment (Gorbachev et al., 2012) 
shows nonzero drag resistivity at charge neutrality, which 
in addition is greatly enhanced by the external magnetic 
field (Titov et al., 2013a). 

Indeed, the outlined physical picture is not universal. 
In fact, it only describes a particular (although often 
dominant) scattering process, where momentum is trans¬ 
ferred from an electron-hole pair in the active layer to an¬ 
other electron-hole pair in the passive layer. Technically, 
this process is described by the leading-order perturba¬ 
tion theory, see Fig. 2, yielding Eq. (15). Higher-order 
processes may result in additional contributions which 
are less sensitive to electron-hole symmetry. These in¬ 
clude the so-called “third-order” drag (Levchenko and 
Kamenev, 2008b; Schiitt et al., 2013) and the effect of 
the correlated disorder (Gornyi et al., 2000, 1999; Schiitt 
et al, 2013; Song et al., 2013), see Secs. 11.D and IV.F. 















In conventional heterostructures, higher-order pro¬ 
cesses remain subleading at least within the temperature 
range where most of the experiments are performed, see 
Sec. II.D. Specifically in the ballistic regime, the dom¬ 
inant contribution to drag is indeed given by Eq. (15) 
[with the corresponding drag resistivity (3a)] and is de¬ 
termined by the nonlinear susceptibility, which in the 
simplest case of energy-independent impurity-scattering 
time r is given by (Kamenev and Oreg, 1995) 

F{q,uj) = — erq^^ 9{vFq — oj). (17) 

TT vpq 

As shown in Kamenev and Oreg (1995), the resulting 
expression (17) for the nonlinear susceptibility is propor¬ 
tional to the imaginary part of the single-layer polariza¬ 
tion operator 

r(q,w) = ^Imn«(9,a;), (18) 

m 

where (for two-dimensional, noninteracting electron gas 
in the ballistic regime) 

lran^{q,U!) = 6{vpq — (jj). (19) 

vpq 

Within the kinetic theory, one can observe Eq. (18) al¬ 
ready at the level of the collision integral (Giuliani and 
Quinn, 1982); hence many authors [see, e.g. Jauho and 
Smith (1993); Shimshoni and Sondhi (1994); Ussishkin 
and Stern (1997); and Zheng and MacDonald (1993)] 
proceed to express Eq. (15) in terms of Imn^(q,w) in¬ 
stead of the nonlinear susceptibility. Under the assump¬ 
tion of energy-independent impurity-scattering time r 
and neglecting intralayer correlations (Flensberg et al., 
1995; Kamenev and Oreg, 1995), such calculations lead 
to the correct result [see Eq. (21) below]. At the same 
time, within such an approach the physics of electron- 
hole asymmetry remains hidden. Generalization to more 
general settings is also nontrivial: the relation (18) is 
by no means a general theorem (Flensberg et al., 1995; 
Kamenev and Oreg, 1995; Narozhny and Aleiner, 2000; 
Narozhny et al, 2012); for explicit examples of the two 
quantities being inequivalent see Secs. Ill and IV. 


2 . Drag resistivity in ballistic samples. 


In the limit of strong screening, >cd ^ 1, one can ap¬ 
proximate (Kamenev and Oreg, 1995) the interlayer in¬ 
teraction propagator (7) by the expression 


V 


R 

12 


Tte^ q 
>ci Ki sinh qd 


( 20 ) 


Gombining the nonlinear susceptibility (18) with the in¬ 
terlayer relaxation rate (15) and Eq. (3a), one finds the 
following expression for the drag resistivity (Flensberg 


and Hu, 1994; Jauho and Smith, 1993; Kamenev and 
Oreg, 1995; Zheng and MacDonald, 1993) 


h 7r2C(3) T2 1 

16 EpiEp2 xiX2kpikp2d'^ 


(21a) 


The same result can be also expressed^ in terms of the 
interlayer relaxation rate (15) [e.g., using Eq. (3a)] 


_-i 71^(3) m _ 1 /21b) 

^ 16 771.2 EpiEp2 X\X2kp\kp2d^ 

Physically, these expressions'^ can be understood based 
on the Fermi Golden Rule [which was explicitly used in 
the solution of the kinetic equation (Jauho and Smith, 
1993; Laikhtman and Solomon, 1990)]. Indeed, there are 
three basic elements that combine into the result (21): 
(i) the phase space available for electron-hole pairs in 
the two layers, which is limited by temperature, hence 
oc r^; (ii) the electron-hole asymmetry, which results 
in the overall smallness of the effect, oc {EpiEp 2 )~^; 
and (iii) the matrix element of the interlayer interaction, 
determining the dependence on the interlayer separation; 
in the ballistic case the matrix element is dominated by 
small-angle scattering (Gramila et al., 1991). 

The drag resistivity (21) - and especially the quadratic 
temperature dependence - is often quoted as the “Fermi- 
liquid” result. However, Eq. (21) was obtained under 
a number of assumptions: (i) xd ^ 1, (ii) d^ £, and 
(iii) T <^Td ^ vp/d ^ Ep/{kpd). The latter assumption 
appears only implicitly and is often overlooked. 

Indeed, substituting the interaction propagator (20) 
and the nonlinear susceptibility (17) into Eq. (15), one 
finds that except for the ^-function in Eq. (17) the fre¬ 
quency and momentum integrals factorize. The exponen¬ 
tial decay of the corresponding integrands allows one to 
estimate the typical values of transferred energy w ~ T 
and momentum g ~ 1/d. Assuming T T^, this yields 
uj < Vpq, which automatically satisfies the 0-function. 
Based on this observation, one may omit the 0-function 
and subsequently extend the integration limits in both 
integrals in Eq. (15) to infinity. The remaining integra¬ 
tion is straightforward and yields Eq. (21). 

At higher temperatures, T Td, the 0-function in 
Eq. (17) is not satisfied automatically. Physically, it rep¬ 
resents kinematic restrictions on the phase space avail¬ 
able to electron-hole pairs associated with predominantly 
small-angle scattering (Gramila et al., 1991). The fre¬ 
quency integration is now cut off at vpq (or Td), rather 


^ Most expressions for pp, (Flensberg and Hu, 1994; Flensberg 
et al., 1995; Jauho and Smith, 1993; Rojo, 1999; Zheng and Mac¬ 
Donald, 1993) can be reduced to Eqs. (21) using the following 
simple relations, valid under the assumptions of this Section: 
Ep = TTn/m, n = Epu, uD = Eprjp, where D = u|,t/2 is the 
diffusion constant and v = m/vr is the density of states. 
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than T, which leads to the linear temperature depen¬ 
dence [first reported in Gramila et at (1991) and Solomon 
and Laikhtman (1991), see also Jauho and Smith (1993), 
and recently rediscovered in Chen et al. (2015)], 


Pd{T » Td) 


h TT^ 1 T 
e^ 360 {kFd)^{>t:dY Ep 


( 22 ) 


This behavior may be observable in samples with either 
nondegenerate 2DEGs or large interlayer separation. In 
the latter case, Td <^T Ep, both layers are perfectly 
described by the Fermi liquid theory which is not synony¬ 
mous with quadratic temperature dependence of trans¬ 
port coefficients. 


3. Plasmon contribution 


Pd 

Pmax plasmons 



FIG. 4 (color online) Schematic illustration for the drag re¬ 
sistivity at high temperatures showing the plasmon peak at 
T ~ Th. The asymptotic dependences are exaggerated for 
clarity. Definitions of the three crossover scales are given in 
the main text. [Reproduced from Chen et al. (2015).] 


The approximate form of the interlayer Coulomb in¬ 
teraction (20) appears justified in the “ballistic” regime 
where the dominant interlayer relaxation processes are 
characterized by relatively large momentum transfers 
oj < vpq. The imaginary part of the single-particle po¬ 
larization operator (19) vanishes at smaller momenta (or 
larger frequencies) making the above calculations con¬ 
sistent. At the same time, approximating the interlayer 
interaction propagator (7) by Eq. (20) one completely 
neglects a possible contribution of plasmon modes, that 
(within the simplest RPA approach) can be found by 
setting the denominator of Eq. (7) to zero. At zero tem¬ 
perature and for w vpq (where Imll^ = 0), the polar¬ 
ization operator is known to be given by (Stern, 1967) 

n(q,tLi) ~ —ne^jinUiP'). 

Using this expression and expanding the bare Coulomb 
potential in small momenta, yields the acoustic (“—”) 
and optical (“-I-”) plasmon modes with dispersions 

LO- = eq\j2'imd/m, oj+ = esj TKnqjm. 

Both of these modes lie outside of the particle-hole con¬ 
tinuum and in the parameter region, where the nonlinear 
susceptibility (17) vanishes. Hence, one may conclude 
that the plasmons cannot contribute to frictional drag. 

However (Flensberg and Hu, 1994), at finite tempera¬ 
tures thermally excited quasiparticles and plasmons may 
coexist in the same parameter region, which may result 
in an additional contribution to drag. In order to ac¬ 
curately describe the plasmon contribution to pp, one 
has to consider intralayer equilibration due to electron- 
electron collisions (Chen et al., 2015) which gives rise to 
two important features: (i) the polarization operator ac¬ 
quires nonvanishing spectral weight within the high fre¬ 
quency part of the spectrum at a; > vpq (Flensberg and 
Hu, 1994); and (ii) the plasmons acquire finite life-time 


(Hruska and Spivak, 2002; Mishchenko et al., 2004) that 
regularizes the pole in the interaction propagator. 

The theory discussed in the preceding Sections in based 
on the implicit assumption that the intralayer equilibra¬ 
tion is the fastest process in the system. Characteriz¬ 
ing inelastic electron-electron scattering by the quasi¬ 
particle lifetime Tee, one finds that the standard theory 
- and hence Eq. (22) - is valid as long as the length 
Tee is much smaller than the interlayer scattering time. 
Tee ^ Tp and at temperatures below the corresponding 
threshold T ^ Tc ^ Epyjkp/{>c‘^d). 

At higher temperatures, T > Te, the system enters the 
collision-dominated regime, where Coulomb drag is dom¬ 
inated by plasmons. In this regime, Chen et al. (2015) 
find a stronger temperature dependence 


The rise of the plasmon contribution to the drag resis¬ 
tivity persists so long as the quasiparticle decay rate re¬ 
mains small compared to the plasma frequency (at the 
wave vector I/d), i.e. up to the third crossover temper¬ 
ature Th Epxjkp/x^\/xd. At temperatures above 
the crossover, T > Th, the electronic system enters the 
hydrodynamic regime of density fluctuations, where 


n 1 Ep 

~ ^ {kpdPixdP^' 


(24) 


This limit can be understood on the basis of the classical 
Navier-Stokes hydrodynamics (Apostolov et al., 2014). 

The resulting temperature dependence of the drag co¬ 
efficient is summarized in Fig. 4. The nonmonotonicity 
oi pp originates from the delicate interplay of various 
scattering channels in the electronic system. Perhaps the 
most striking feature of the theory of Chen et al. (2015) is 
that intralayer collisions promote stronger drag. Indeed, 
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should one naively continue Eq. (22) up to the tempera¬ 
tures of the order Th one would underestimate the actual 
maximum value of pc by ^Jkpd ^ 1 . 


C. Effects of potential disorder 



In ballistic samples, potential disorder played a very 
limited role. In fact, the resulting drag resistivity (21) 
is independent of the impurity scattering time r. In dif¬ 
fusive samples with i only small momenta q^Xjl 
can be transferred between the layers. Typically this re¬ 
sults in a small contribution to the drag resistivity, which 
in ballistic samples can be neglected. This is not always 
the case - at low enough temperatures drag is dominated 
by mesoscopic fluctuations which are mostly due to pro¬ 
cesses with small momentum transfers, see Sec. III. 

Coulomb drag in diffusive systems was considered in 
Zheng and MacDonald (1993) using the memory func¬ 
tion formalism and in Kamenev and Oreg (1995) using 
the diagrammatic technique. To the lowest order in inter¬ 
layer interaction, one can use the Kubo formula analysis 
(Kamenev and Oreg, 1995; Narozhny and Aleiner, 2000) 
to derive the expression for the drag conductivity (S is 
the area of the sample) 

where numerical subscripts indicate spatial coordinates, 
and are implied to be integrated over. Averaging over 
disorder restores translational invariance. In the absence 
of interlayer disorder correlations [this special case was 
considered in Gornyi et al. (1999)], the nonlinear suscep¬ 
tibilities in each layer have to be averaged independently 
of each other. Then one recovers the drag relaxation 
rate (15), where each quantity should be understood as 
disorder-averaged, i.e. (/Ifg) —> T^{q). 


1. Drag resistivity in diffusive regime 


In the diffusive regime, the nonlinear susceptibility F 
can be found from the Ohm’s law (Landau et al, 1984), 

j = dE — eD'Vn, (26) 


where d is the conductivity matrix and D is the diffusion 
coefficient (in two dimensions D = v^pT 12). Combining 
Eq. (26) with the continuity equation, one finds the linear 
response of the carrier density n to the electric field E 


{n{q,oj)) 


1 iq‘^CF°‘^E^{q, uj) 
e —ioj + Dq^ 


(27) 


where (...) indicates averaging over disorder. Nonlinear 
response follows from the density dependence of the con- 


FIG. 5 (Color online) Disorder averaging of the nonlinear 
susceptibility (Kamenev and Oreg, 1995). The dotted lines 
represent the diffuson ladder (Altshuler and Aronov, 1985). 


ductivity jdc = Re {da/dn) n{q,uj)E{—q, —oj), and yields 


in 


2v d{a^^) s 
e dn ^ Lu^ + D'^q'^ 


(28) 


In the absence of a magnetic field (cr“^) = crS°^^, and 
the nonlinear susceptibility (28) is parallel to q. The 
disorder-averaged conductivity is linear in the carrier 
density, da°‘^fdn « jn. As result. 


(F) = 2q 


evD 

Ep 


UjDq^ 
uP' -I- D'^q^ 


(29) 


This expression can be recast into two equivalent forms. 
Noting the similarity between Eq. (29) and the standard 
diffusive form of the polarization operator (Altshuler and 
Aronov, 1985; Smith and Jensen, 1989) 

one finds (Kamenev and Oreg, 1995) 

eD 

(E) = 2q—Imn^(q,w). (31) 


Eurthermore, one can emphasize the fact, that the den¬ 
sity dependence of the conductivity cr is a manifestation 
of electron-hole asymmetry by rewriting the fraction in 
Eq. (31) as (Narozhny and Aleiner, 2000) 


(r) = 2eqD ^^^j'’^^\ mIl^{q,uj). (32) 

a/r 

In the simplest case^, this expression can be obtained 
directly from Eq. (28) by noticing that da°‘^/dn = 
{da°‘^/d^){dfj,/dn) = (9cr“^/i9^)(l/i/) and using the 
Einstein relation. The same result can be found eval¬ 
uating diagrams shown in Fig. 5. 

The diffusive approximation for the interlayer interac¬ 
tion follows from Eqs. (7) and (30). Focusing on small 
momenta 9 ^ 1 /d, one can obtain alternative expressions 
for the interaction propagator by either expanding the 
bare matrix element (5) in small qd (and subsequently 
limiting the momentum integration from above) or keep¬ 
ing the exponential in Eq. (5) intact, leaving the mo¬ 
mentum integral converging in the ultraviolet. The for¬ 
mer approach was taken in Narozhny and Aleiner (2000). 
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Generalizing to inequivalent layers one finds 
1 {—iui + Diq^){—iuj-\-D2q^) 


_ _ 
■^12 — 


where 


{viDi + V 2 D 2 )[-iu: + (1 + >c*d)D*q^] 


(33a) 


= 47re^ 


ViV2 


D* = 


{Vi + V2 )DiD2 


Vl + V 2 ViDi + V2D2 

The latter alternative was taken in Kamenev and Oreg 
(1995), where in addition (just as in the ballistic case) 
the limit xd ^ 1 was used. As a result, the interaction 
propagator takes the form 

—iuj + Diq^ —iu) + D2q^ 


_ __ 
■^12 — 


Tre^q 


>£■1 >£’2 sinh qd Diq^ T* 2 <Z^ ^ ^ 

With logarithmic accuracy, the resulting drag coeffi¬ 
cient is independent of the distinction between the two 
and can be written as (Kamenev and Oreg, 1995) 

h 


(Pd) = 




1 


■In 


2T' 


(34a) 


12EpiEp2 xiX2kpikp2ti(-2d?‘ 

The only difference between using the two expressions for 
the interaction propagator in Eq. (33) is the exact value 
of Tq. Using Eq. (33a) in the limit x*d^ 1, one finds 

^ 1^101^202 
° {viDi + V2D2)d’ 

while Eq. (33b) leads to (Kamenev and Oreg, 1995) 

To = min{>iri£)i, >ir2Tl2}/'^- 

Both expressions are of the same order of magnitude and 
coincide for the case of identical layers. 

The result can be expressed also in terms of the inter¬ 
layer relaxation rate (Zheng and MacDonald, 1993) 


1 

td 


TT^ ni 


T2 


Tn 


1 


12 777.2 EpiEp2 2T H\H2kp\kp2^i^2d‘^ 
Equivalently, one can use Eq. (32) and express the drag 
conductivity (Narozhny and Aleiner, 2000) as (here the 
layers are assumed to be identical for simplicity) 


an = 


[hTf 




To 


(35) 


h 3 g'^{xdY 

where the derivative highlights the crucial role of the 
electron-hole asymmetry in the leading-order drag effect. 

The diffusive result for the drag resistivity (34) appears 
to be rather similar to its ballistic counterpart Eq. (21). 
Indeed, disregarding the numerical prefactors and the 
logarithm in Eq. (34), one finds 


diff / I 

Pd IP 


bal 

D 


dV(^i^ 2 ) 


(36) 


This relation may serve as an a posteriori justification 
for the statement that the drag effect in samples with 
d ^ £ is dominated by ballistic propagation of carriers 
with momenta q d~^. Carriers with small mo¬ 

menta q £~^ also participate in drag, but their contri¬ 
bution is small [according to Eq. (36)] and is typically 
neglected. 



FIG. 6 (Color online) Leading weak localization corrections 
to the nonlinear susceptibility (Kamenev and Oreg, 1995). 
The black, parallel dotted lines represent the diffuson ladder 
(Aleiner et al, 1999; Altshuler and Aronov, 1985). The blue, 
crossing lines represent the Cooperon (Gorkov et al, 1979). 


2. Weak localization corrections 


The nonlinear susceptibility (29) and drag coefficient 
(34) were obtained as the leading approximation in the 
standard perturbation theory of disordered metals (Alt¬ 
shuler and Aronov, 1985), controlled by the large param¬ 
eter g = 25.8kf2/i?n representing the dimensionless con¬ 
ductance of the layers (with £?□ being the layer (sheet) 
resistance). Within the assumptions adopted in this Sec¬ 
tion"^ g ^ vD ^ kpi ^ EpT ^ 1. 

The next-order terms in the perturbation theory are 
known as quantum corrections to transport coefficients 
(Aleiner et al, 1999; Altshuler and Aronov, 1985). Phys¬ 
ically, they describe leading interference processes that 
arise in the course of subsequent scattering events. Al¬ 
though the resulting contribution to transport is propor¬ 
tional to a small factor 1/g, quantum corrections domi¬ 
nate the temperature and magnetic field dependence of 
transport coefficients at low temperatures. 

To the leading order in 1/g, one may distinguish 
three types of corrections: (i) interference between self- 
intersecting, time-reversed scattering paths leads to a 
positive correction to resistivity, known as the weak lo¬ 
calization correction (Abrahams et al., 1979; Altshuler 
et al, 1980; Gorkov et al., 1979); (ii) coherent scattering 
off Friedel oscillations yields the Altshuler-Aronov cor¬ 
rection (Altshuler and Aronov, 1979; Einkelstein, 1983, 
1984; Zala et al., 2001); and (iii) in small, mesoscopic 
samples interference between scattering paths gives rise 
to universal conductance fluctuations (Altshuler, 1985; 
Lee and Stone, 1985). The latter effect has a direct 
counterpart in double-layer systems, namely mesoscopic 
fluctuations of Coulomb drag discussed in Sec. III. At 
the time of writing, no qualitative interference effect due 
to electron-electron interaction has been identified for 
drag measurements. At the technical level, the third- 
order drag effect (see the following Section) bears certain 
resemblance to the Altshuler-Aronov diagrams (Gornyi 
and Narozhny, 2014). Here we discuss the weak localiza¬ 
tion correction to Coulomb drag (Flensberg et al, 1995; 
Kamenev and Oreg, 1995). 
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In the absence of interlayer disorder correlations [such 
effects were discussed in Gornyi et al. (1999)], impurity 
scattering is confined to each individual layer. It should 
come as no surprise, that the same mechanism behind the 
weak localization correction to single-layer conductivity 
(i.e., interference between time-reversed, self-intersecting 
paths) yields a correction to the nonlinear susceptibility. 
Technically, this interference mechanism is described by 
a “maximally crossed” element of the diagram technique 
known as the Cooperon (Gorkov et al, 1979). Diagrams 
for the corresponding corrections to the nonlinear suscep¬ 
tibility are shown in Fig. 6 [further corrections, e.g. two- 
Gooperon diagrams, considered in Flensberg et al. (1995) 
and Kamenev and Oreg (1995) were found to be sublead¬ 
ing]. The resulting nonlinear susceptibility is given by 


(r) = 2q 


Ep 


ujD{oj, q)q^ 
up- + q)q'^ ’ 


(37) 


where the renormalized diffusion coefficient in two dimen¬ 
sions in (Gorkov et al., 1979) 


D{uj,q) = D 


1 

nkpi 



(38) 


and is the dephasing time (Altshuler et al., 1980). The 
result (37) is valid in the first order in 5D = D{uj, q) — D. 

The resulting leading-order weak localization correc¬ 
tion to Goulomb drag is (Kamenev and Oreg, 1995) 


SpD ^ 1 ^ _ 1 1 

Pd irkpiii 2 Tti 'KkF2f-2 ‘2 Tt2' 


(39) 


where po is given by Eq. (34). The result (39) is simi¬ 
lar to the weak localization corrections in 2D (Altshuler 
et al, 1980; Gorkov et al, 1979), except that in Eq. (39) 
the logarithmic singularity is cut by temperature rather 
than by the dephasing time. 

In conventional 2DEG, weak localization effects re¬ 
sult in a dependence on a weak magnetic field (Altshuler 
et al, 1980). Here, the characteristic scale of the mag¬ 
netic field would be He j (eD). Similar scale describes 
intralayer interaction corrections to magnetoresistance 
(Altshuler and Aronov, 1985), making the weak local¬ 
ization corrections to the drag coefficient hard to observe 
experimentally (Kamenev and Oreg, 1995). 


D. Third-order drag effect 

The leading contribution to Coulomb drag, Eqs. (15), 
(25), describes the effect to the lowest order in the inter¬ 
layer Coulomb interaction, see Fig. 2. Since the particles 
belonging to different layers interact through a layer of 
an insulating material, certain weakness of the effective 
interaction is intuitively expected. In many-body elec¬ 
tron systems the Coulomb interaction is usually screened 
and the perturbative analysis gives reasonable account of 




FIG. 7 (color online) Typical diagrams describing higher- 
order drag effects. Left: third-order drag (Levchenko and 
Kamenev, 2008b). Right: the effect of interlayer disorder 
correlations (Gornyi et al., 1999; Hu, 2000a). 


most basic observable quantities (Altshuler and Aronov, 
1985; Ziman, 1965). Consequently, the vast majority of 
theoretical studies of Coulomb drag are devoted to the 
investigation of the lowest-order effect. Notable excep¬ 
tions are given by the studies of the interlayer corre¬ 
lated states, either in the context of the quantum Hall 
devices (Girvin and MacDonald, 1997; Kim et al, 2001; 
Stern and Halperin, 2002; Stern et al, 2000; Yang, 1998; 
Yang and MacDonald, 2001) or quantum wires (Klesse 
and Stern, 2000; Nazarov and Averin, 1998), as well as 
strongly correlated intralayer states, such as Wigner crys¬ 
tals (Baker and Rojo, 2001; Braude and Stern, 2001) or 
Anderson insulators (Raikh and von Oppen, 2002). 

The “single-particle” drag resistivity, Eqs. (21), (34), 
is determined (besides the interlayer interaction) by the 
quasiparticle phase space, electron-hole asymmetry (see 
Sec. H.B), and disorder effects (see Secs. H.C and HI). At 
T = 0 or at a point of exact electron-hole symmetry (e.g., 
in neutral graphene, see Sec. IV), these factors may con¬ 
spire to nullify the effect. Then pD may be determined by 
higher orders of the perturbation theory, implying that 
saturation of drag resistivity at low temperatures should 
not necessarily point towards a strongly correlated state. 

To the third order in interlayer interaction (see Fig. 7 
for the “skeleton” diagram). Coulomb drag was first dis¬ 
cussed in Levchenko and Kamenev (2008b) in the diffu¬ 
sive regime, Tt 1. It was shown that the third-order 
drag contribution remains finite at zero temperature®: 

pg^(T < h/r) = {).27{h/e^)g-^{>cd)-^. (40) 

This surprising result was attributed (Levchenko and 
Kamenev, 2008b) to the singular behavior of matrix el¬ 
ements in the diffusive regime. In single-layer systems, 
similar enhancement of the matrix elements leads to sin¬ 
gular interaction effects (Altshuler and Aronov, 1985). 
Here, the divergence of the matrix elements is compen¬ 
sated by the smallness of the phase space yielding the 
T-independent contribution to the drag resistivity. 


More precisely, the result is valid down to lowest temperatures 
T ~ . Below this scale the diffusive approximation 

breaks down. 
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The third-order effect (40) does not rely on electron- 
hole asymmetry (technically, the third-order diagram in 
Fig. 7 contains four-point vertices instead of the trian¬ 
gular vertices in Fig. 2). Hence, is independent of 
Ep. This provides an additional explanation of the T- 
independent result (40): in the diffusive regime there is 
no other scale for a temperature dependence. 

Another contribution to drag that is insensitive to 
electron-hole symmetry is due to interlayer disorder cor¬ 
relations (Gornyi et al, 1999; Hu, 2000a). For tempera¬ 
tures higher than the inverse interlayer coherence time, 
but still in the diffusive regime, <^T <^t, one hnds 

Pd{t~^C T < Tt”^) ~ {hle^)(kpd^xtj~'^ In(TTg), (41) 

which might dominate over Eq. (34). 

While the above higher-order effect have not been ob¬ 
served in semiconductor samples, they may provide an 
explanation of the observed nonzero drag resistivity in 
neutral graphene (Gorbachev et al., 2012), see Sec. IV. 



FIG. 8 (a) Two diagrams contributing to the transconduc¬ 
tance that are second order in tunneling matrix element de¬ 
noted by a cross. Full lines with arrows are electron Green 
functions, dashed lines represent diffusons, and wavy lines 
screened interactions. Two additional diagrams with arrows 
in the opposite direction should be included. The numbers 
indicate the layer index, (b) Examples of diagrams contribut¬ 
ing to the transconductance that are fourth order in tunnel¬ 
ing. Diagrams with interaction lines connecting “upper” and 
“lower” Green functions, as well as diagrams with an opposite 
direction of electron lines are also implicit. [Reproduced from 
Oreg and Kamenev (1998).] 


E. Transconductance due to tunneling bridges 

A qualitatively different mechanism of transconduc¬ 
tance takes place in the double layer systems with point¬ 
like shortages (or bridges) - places where electrons may 
tunnel between the two layers (Oreg and Halperin, 1999; 
Oreg and Kamenev, 1998). Such bridges can be present 
in metallic double-layer systems due to device fabrica¬ 
tion imperfections, or they can be introduced on purpose 
(Giordano and Monnier, 1994). In sharp contrast to pure 
Goulomb drag, the interplay between interlayer tunnel¬ 
ing and Goulomb interaction leads to a transconductance 
that is nonvanishing at zero temperature and negative for 
carriers of the same charge: the induced current flows in 
the direction opposite to the driving one. The origin of 
the effect is in the interlayer exchange interaction, which 
is possible due to wave functions overlap at the bridges. 

To illustrate the physics of the effect consider a bunch 
of electrons injected into the active layer. Some of them 
tunnel through a bridge to the passive layer, while the 
others continue to propagate in the initial direction. 
Without interlayer electron-electron interaction, the tun¬ 
nelled electrons “forget” about the direction of their ini¬ 
tial momentum: the wave front propagating in the pas¬ 
sive layer is spherically symmetric; thus no net current 
is induced. The repulsive interaction between the tun¬ 
neled electrons and the electrons remaining in the active 
layer makes the former move in a direction opposite to 
the injected current. This leads to a negative transcon¬ 
ductance. This lowest order in tunneling mechanism, 
which is illustrated diagrammatically in Fig. (8a) is dom¬ 
inant only at temperatures that are not too small. At 
lower temperatures the leading mechanism involves co¬ 
herent tunneling to the second layer and back to the first 


one accompanied by intralayer Goulomb interactions, see 
Fig. (8b). Since the exchange contribution is not pro¬ 
portional to a small electron-hole asymmetry factor, its 
absolute value may well overcome the standard Goulomb 
drag even for a small tunneling rate. 

The interlayer tunneling rate may be described by a 
mean intralayer lifetime ri 2 , which is related to the inter¬ 
layer tunneling conductance per unit area crj_ = It\ 2 . 

Assuming weak tunneling, ri2 3> r (t is the mean elas¬ 
tic scattering time within each layer), one finds the fol¬ 
lowing regimes determined by the three energy scales 
Tf2^ < Kdr^.^ < T-i. 

At high temperatures, T > , the motion of the wave 

packets is ballistic and the exchange contribution to the 
transconductance is negative and independent of T 


TT 1 VpT^ 

h 32 xd dTi 2 


(42) 


At temperatures T < the diffusive character of the 
electron motion should be taken into account. To the 
leading order in the tunneling amplitude, one finds the 
following temperature dependent contribution: 


1 In(xd) 1 
h 247r xd Tti 2 


(43) 


Singular temperature dependences are not uncommon in 
diffusive systems. What is less common is the fact that 
the divergence is so pronounced, ~ 1/T, instead of being 
logarithmic (cf. the Altshuler-Aronov corrections to the 
conductivity of 2D systems). 

The temperature dependence becomes even more sin- 
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gular at smaller temperatures, T < >€dTi.^, where 


e" 3C(3) ln(rri2) 

h StT^ (Tti2)2 


(44) 


The divergences in Eqs. (43) and (44) are cut off by the 
Thouless energy Et = DjLp' {L is the system size). For 
large systems, L ^ \/Dti 2 , there is an additional tem¬ 
perature range Et <T < ■ Here the exchange con¬ 

tribution to the cr_D is dominated by multiple tunneling 
processes, leading to the form 


In — 
h StT^ Tti2 


(45) 


This expression differs from the Altshuler-Aronov cor¬ 
rection to the 2D conductivity only by a factor of 1/4, 
which reflects the peculiarity of the drag measurement 
setup: the current flows only in one half of the system 
and the potential is measured in the other half. 

To reveal above quantum corrections ajj in an actual 
experimental situation one should take into account a 
classical contribution to the transresistance. It originates 
from currents which may flow in the second layer in the 
presence of interlayer tunneling. These currents lead to 
an additional voltage drop in the second layer, which 
modifies the measured transresistance. Such currents ex¬ 
ist either in the presence of many point-like bridges, or 
in the case of a uniform tunneling through a sufficiently 
thin barrier. In the latter case the tunneling rate is sen¬ 
sitive to a Fermi-surface mismatch between the layers, 
and thus may be affected by a gate voltage or an in¬ 
plane magnetic field (Berk et al., 1995; Boebinger et al., 
1991). The case of many point-like bridges may be vi¬ 
sualized by a network of classical resistors. Calculations 
done in the continuous limit (Raichev, 1997) show that 
for relatively long samples, L > y/D ti 2 , approximately 
half of the current applied to the first layer eventually 
leaks out to the second one. In this case the transresis¬ 
tance practically reduces to the resistance of a single layer 
of a doubled width. Such classical effect however leads to 
a weakly temperature dependent transresistance. Thus 
the exchange contribution, although possibly smaller in 
its absolute value, may be distinguished due to its singu¬ 
lar temperature dependence. 


F. Comparison to experiment 

The theory outlined in the preceding sections describes 
an idealized phenomenon of mutual friction between two 
two-dimensional electron systems. The electrons were 
assumed to belong to a parabolic band, with energy- 
independent impurity-scattering time and negligible in¬ 
tralayer correlations. Clearly, such assumptions can be 
realized in any experimental sample only approximately. 

Coulomb drag between two two-dimensional electron 
gases was first observed by the group of J. Fisenstein 
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FIG. 9 Left panel: temperature dependence of the drag resis¬ 
tivity GaAs double-quantum-wells. The additional scale on 
the right provides the corresponding values of the momentum- 
transfer rate (see main text for more details). The inset shows 
an idealized energy diagram for a double-quantum-well struc¬ 
ture indicating the ground subband energy Eq and the Fermi 
energy Ep. [Reproduced from Gramila et al. (1991).] Right 
panel: Temperature dependence of the interlayer momentum 
transfer rate divided by . The three sets of data were mea¬ 
sured in samples with interwell barrier widths of 175A, 225A, 
and 500A. [Reproduced from Gramila et al. (1992).] 


(Fisenstein, 1992; Gramila et al., 1991) in GaAs double¬ 
quantum-wells, see Fig. 9. Detailed comparison of the 
experimental data to the quantitative predictions of the 
Goulomb drag theory showed that the latter accounts for 
about 50% of the measured values®. This was judged 
as sufficient evidence of the relevance of the Goulomb 
mechanism of frictional drag. Also, the overall reduction 
of the drag resistance with the increase of the interwell 
barrier width (see Fig. 9) was in rough agreement with 
Fq. (21b). At the same time, the data (see the right 
panel in Fig. 9) show noticeable deviations from the T® 
behavior predicted by Fqs. (21) and (34), indicating that 
other scattering mechanisms might also be important. 

One additional mechanism (Gramila et al., 1991) is 
due to electron-phonon interaction. This suggestion was 
developed theoretically in Badalyan and Rossler (1999); 
Bpnsager et al. (1998a,b); and Tso et al. (1992) and ex¬ 
perimentally in Jorger et al. (2000a); Noh et al. (1999); 
and Rubel et al. (1995), see Fig. 10. 

Further corrections to the single-particle Goulomb 
mechanism are associated with the plasmon contribu¬ 
tion. As shown in Flensberg and Hu (1994), plasmons 
are expected to be most important at intermediate tem¬ 
peratures, T ~ O.STp. This prediction was tested experi¬ 
mentally in Hill et al. (1997), see the left panel in Fig. 11, 


® The momentum relaxation time reported in Gramila et al. (1991) 
is twice smaller than Eq. (21b). In addition, the paper cited 
unpublished calculations of MacDonald, Gramila, and Eisenstein 
involving a more realistic modeling of finite-width quantum wells. 
In particular, these calculations were reported to include vertex 
corrections to the RPA interaction propagator (7). Hence, it is 
difficult to judge whether that factor of 2 has played any role in 
the actual analysis of Gramila et al. (1991). 
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FIG. 10 Left panel: calculated jT'^ (solid curves) com¬ 
pared to the data of Gramila et al. (1991). A less opti¬ 
mal choice of a fitting parameter yields resnlts shown by 
the dashed cnrves. Inset: calculated contribution of virtual- 
phonon exchange processes to r^j^/T^). [Reproduced from 
Tso et al. (1992).] Right panel: calculated po/T^ for vari¬ 
ous values of the phonon mean free path and d = 500A (solid 
curves). The dots show the data of Gramila et al. (1992). The 
dotted line represents the contribution of the modified plas- 
mon pole. Inset: the crossover of the T® to T temperature 
dependence. [Reproduced from Bpnsager et al. (1998a).] 


and in Noh et al. (1998). While the theoretical results 
show qualitative agreement with the data, discrepancies 
persist. Taking into account many-body correlations [see, 
e.g., Swierkowski et al. (1995)] improves the agreement, 
but further advances in many-body theory are necessary 
before a more precise quantitative description of the cor¬ 
relation effects in double-layer structures is achieved. 

The discrepancies between the simple single-particle 
description and laboratory experiments are by no means 
universal, especially since many measurements were per¬ 
formed in very different systems. One of the first drag 
experiments (Solomon et at, 1989) has been performed 
on a hybrid 2D-3D system. This device showed consider¬ 
able thermoelectric effects masking the purely Coulomb 
contribution to drag. Experiments on electron-hole sys¬ 
tems (Sivan et ai, 1992) showed behavior that could not 
be accounted by neither the phonon, nor plasmon cor¬ 
rections. Instead, generalized RPA (taking into account 
exchange processes to all orders) (Tso et al., 1993) ap¬ 
pears to yield satisfactory agreement with observations of 
Sivan et al. (1992) at low temperatures. Apparently, the 
traditional RPA overestimates screening which results in 
the underestimated drag resistivity. 

Experiments on dilute 2D hole systems (Pillarisetty 
et ai, 2002, 2004) show marked enhancement of the 
drag resistivity, along with the stronger temperature de¬ 
pendence (empirically, pu oc at low temperatures, 
followed by a crossover towards a sublinear tempera¬ 
ture dependence at T ~ Ep). These systems are char¬ 
acterized by rather high values of the interaction pa¬ 
rameter Ts — 20-40 and also exhibit signs of a metal- 
insulator transition in single-layer measurements (Pil¬ 
larisetty et al., 2005b). The data obtained in Pillarisetty 
et al. (2002) are not explained by taking into account 
corrections due to phonons (Bpnsager et al., 1998a), plas- 
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FIG. 11 Left panel: measured po/T^ for various values of the 
carrier density, ni = n 2 . Dashed lines represent the results of 
Flensberg and Hu (1994) adjusted for the sample parameters 
of the experiment. Solid lines show the results of additional 
calculations taking into account intralayer many-body corre¬ 
lations within the Hubbard approximation. (Reproduced from 
Hill et al. (1997).] Right panel: measured po/T^ for differ¬ 
ent carrier densities. Inset: (a) peak position temperature vs 
matched layer density; (b) po vs density ratio for T — 860, 
730, and 600mK. [Reproduced from Pillarisetty et al. (2002).] 


mons (Flensberg and Hu, 1994), or many-body effects 
(Swierkowski et ai, 1995), as follows from the density 
dependence of the measured drag, illustrated in the right 
panel of Fig. 11. The lack of adequate theoretical descrip¬ 
tion of these experiments is not surprising, given that the 
regime of relatively high rg remains an unsolved problem 
in single-layer (bulk) systems as well. 

Croxall et al. (2008) and Das Gupta et al. (2008) re¬ 
port anomalous drag in electron-hole bilayers. Below 
T = IK, the measured drag resistivity exhibits an up¬ 
turn that may be followed by a downturn, although pp 
does not seem to vanish for T —)■ 0. The observed upturn 
may indicate exciton formation (Hu, 2000b; Vignale and 
MacDonald, 1996), however neither the observed viola¬ 
tion of Onsager reciprocity, nor the apparent downturn 
at lower temperatures are anticipated by the theory. The 
effect of density imbalance on the drag upturn was stud¬ 
ied in Morath et al. (2009). The data were interpreted in 
terms of a pairing-fluctuations mechanism based on the 
theory of Hwang and Das Sarma (2008b). The theory 
accounts for most qualitative features of the effect, how¬ 
ever the predicted peak in at equal layer densities was 
not observed in experiment (Morath et al, 2009). 

Further experiments demonstrate interesting correla¬ 
tion effects such as Wigner crystallization in quantum 
wires (Yamamoto et ai, 2006, 2012), exciton forma¬ 
tion in electron-hole bilayers (Seamons et al., 2009), or 
quantum Hall effect (Girvin and MacDonald, 1997; Lilly 
et al., 1998), see Sec. VILA. Glearly these phenom¬ 
ena cannot be described by the simple theory presented 
in this Section. At the same time, single-particle ef¬ 
fects are still important at relatively low temperatures 
{T < 0.2Tf) in traditional semiconductor heterostruc- 














16 


tures hosting two-dimensional electron systems and even 
more so in graphene (see Sec. IV), where interlayer sep¬ 
aration can be as small as several interatomic distances 
(Gorbachev et al, 2012). 


1. Phonon effects 


Electrical resistivity due to electron-phonon scattering 
is a standard topic in condensed matter physics (Ziman, 
1965). At temperatures higher than the Debye frequency 
T ^ ludj it exhibits linear behavior p (x T, that is ob¬ 
served in a wide class of materials including high-mobility 
2DEG (Stormer et al., 1990) and graphene (Efetov and 
Kim, 2010). At low temperatures T ^ uo [in low density 
electron systems the crossover occurs at a lower scale, the 
so-called Bloch-Griineisen temperature T Tbg < 
(Stormer et al., 1990)] the phonon contribution is rapidly 
decreasing as p oc T® in metals (Bloch, 1930; Griineisen, 
1933) and heterostructures (Price, 1984; Stormer et al., 
1990) and as p (xT* in graphene (Efetov and Kim, 2010; 
Hwang and Das Sarma, 2008a). 

Qualitative physics of electron-phonon interaction in 
semiconductor double-quantum-well heterostructures is 
captured by the following interaction Hamiltonian 

= ^ E E^AV(Q)^’A,v(gQ 

^ X,X']k Q;7] 


bU-Q) + KiQ)\ci{k)cy{k + q). (46) 


Here Q = {q, q^) is the 3D wave vector of a phonon with 
polarization q, k is the 2D electron wave vector, M’^y 
is the bulk electron-phonon matrix element corrected by 
the subband form-factor 


Fx.,y{q.) = I dzUz)CyF'^^\ ( 47 ) 

— OO 

where ^x{z) is the bound state wave function associated 
with the quantized motion in the subband A. This Hamil¬ 
tonian was used to study effects of interaction between 
electrons and longitudinal optical phonons in Das Sarma 
and Mason (1985) and to calculate quasiparticle proper¬ 
ties in weakly polar 2DEG in Jalabert and Das Sarma 
(1989). In double-layer systems, the Hamiltonian (46) 
was used to describe interlayer interaction mediated by 
acoustic phonons in Bpnsager et al. (1998a) and Zhang 
and Takahashi (1993) and by optical phonons in Hu 
(1998). 

Electrons experience the phonon-mediated interaction 
(46) alongside the Coulomb interaction. The propaga¬ 
tor of the effective interlayer interaction can be obtained 
within the RPA (Bpnsager et al., 1998a; Jalabert and 
Das Sarma, 1989; Zhang and Takahashi, 1993) similarly 


to Eq. (7). The result can be represented in the form 
Di 2 = (Vi 2 + Di 2 )/e{q,oj), where D 12 is the propagator 
of the phonon-mediated interaction and e{q, oj) is the ef¬ 
fective dielectric function for interlayer interactions that 
is also determined by the sum of the Coulomb interac¬ 
tion (5) and the phonon propagator. Thus the phonon 
and Coulomb mechanisms are generally not independent 
of each other. However, the Coulomb interaction con¬ 
tributes only to small momentum transfers (6), while the 
phonon contribution peaks at q ^ 2kp (Bpnsager et al., 
1998a). Neglecting interference between the two, one 
can estimate the effect of phonon-mediated interaction 
by considering only the phonon part D 12 —> D 12 /e{q,ui). 

A simple analytical estimate for the strength of the 
phonon-mediated interaction in GaAs/AlGaAs systems 
was suggested in Bpnsager et al. (1998a). In this ma¬ 
terial, electron-phonon interaction is due to the defor¬ 
mation potential and piezoelectric effect. It turns out, 
that the deformation mechanism dominates (except for 
very low electron densities). Assuming infinite phonon 
mean free path, the corresponding (unscreened) effective 
interaction has the form 


vkFCix/c^q^ - 


(48) 


where c; is the velocity of longitudinal acoustic phonons 
and Cdp ~ 2.7 x 10“^A:F/(I0®cm“^). The smallness of 
electron-phonon coupling constants implies weakness of 
the phonon-mediated interlayer interaction as compared 
to the Coulomb interaction. However, the effective inter¬ 
action (48) diverges near u) « ciq leading to a logarithmic 
divergence in the drag resistivity. Although this diver¬ 
gence is removed by either dynamic screening or phonon 
relaxation, the above argument illustrates the reason be¬ 
hind the relative strength of the phonon-mediated inter¬ 
layer interaction. 

Detailed calculations of the phonon-mediated drag re¬ 
sistivity have been performed numerically by several au¬ 
thors. Tso et al. (1992) showed that combining the 
phonon and Coulomb mechanisms of mutual friction ac¬ 
counts for the nonparabolic temperature dependence ob¬ 
served in GaAs/AlGaAs devices (Gramila et al., 1991), 
see the left panel of Eig. 10. A refined discussion of the 
phonon mechanism was given in Bpnsager et al. (1998a), 
see the right panel of Eig. 10. It was shown, that the tem¬ 
perature dependence of the phonon contribution to drag 
exhibits a crossover from linear to T® behavior around the 
Bloch-Griineisen temperature (see the inset in Fig. 10), 
explaining the peak in the drag resistivity, Fig. 9. In ad¬ 
dition, it was shown that there exists a collective mode 
that can be found setting e{q, uj) = 0. This mode is simi¬ 
lar to the usual plasmon and results from coupling of the 
electrons from both layers to the phonons with oj ~ ciq. 
A similar mode resulting from interaction between elec¬ 
trons and optical phonons was discussed in Giiven and 
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Tanatar (1997a,b). A detailed analysis of the mutual 
friction due to optical phonons is given in Hu (1998). 

2. Interlayer interaction beyond RPA 

The expression (7) for the dynamically screened in¬ 
terlayer Coulomb interaction has been obtained within 
the RPA. While capturing the qualitative physics of the 
effect, this representation is by no means exact. In par¬ 
ticular, RPA-based calculations seem to underestimate 
the value of pn as compared to experimental data (Sivan 
et al, 1992). A pedagogical discussion of the RPA and 
possible approaches to interacting many-body systems 
that go “beyond” the RPA can be found in Giuliani and 
Vignale (2005). Most of these approaches are not para¬ 
metrically justihed. The results of the calculations are 
typically compared to either experimental data or com¬ 
puter simulations. 

Coulomb drag between electron and hole layers within 
the generalized RPA approach was considered in Tso 
et al. (1993). The resulting po is about twice larger 
than that calculated within RPA, but still about twice 
smaller than the experimental data. Furthermore, it 
was understood in Swierkowski et al. (1995) that the 
true temperature dependence oi pn should exhibit a 
crossover from the dependence at low temperatures 
to a power-law at higher temperatures. However, the 
local field approach [or the Singwi-Tosi-Land-Sjdlander 
method (Singwi et al., 1968)] used in this work still fails 
to reproduce Pd{T) measured in Sivan et al. (1992), al¬ 
though yields roughly the same magnitude of the effect 
(in contrast the above RPA and generalized RPA cal¬ 
culations). This approach was further extended to drag 
between two 2DEG in Swierkowski et al. (1996, 1997). 
The results of that work suggest that many-body corre¬ 
lations enhance interlayer interaction and improve agree¬ 
ment with experiments. Nevertheless, experiments [see, 
e.g.. Hill et al. (1997) and Fig. 11] show, that existing 
theoretical methods are still incapable of providing pre¬ 
cise quantitative description of real systems. 

A detailed consideration of Coulomb drag resistivity 
based on an extrapolation of Fermi-liquid-based formu¬ 
las to the region where intralayer correlations are strong 
has been carried out by Hwang et al. (2003) in an at¬ 
tempt to address the striking data of Pillarisetty et al. 
(2002) in low density and high mobility hole bilayers. 
The observed drag was two to three orders of magni¬ 
tude larger than previously reported values. The calcu¬ 
lations of Hwang et al. (2003) were different from that 
leading to Eq. (21a) in several points, all of them lead¬ 
ing to an increase of the drag resistivity: (i) Hubbard 
approximation was employed to obtain polarization op¬ 
erator, which accounts for the exchange-driven local field 
corrections; (ii) experimentally measured dependence of 


conductivity on density was used to extract electron-hole 
asymmetry factor; (hi) large-momentum transfer compo¬ 
nent was included to calculate drag; (iv) Hnite thickness 
of quantum wells was included to calculate form-factors 
of Coulomb matrix elements; (v) lastly, phonon contri¬ 
bution was added. Combining all these factors, Hwang 
et al. (2003) were able to account for most of the results 
of the measurements within a Eermi liquid approach. 


G. Single-particle drag in magnetic field 

The semiclassical Drude model described by Eqs.(2) 
predicts that the drag resistivity is independent of the 
magnetic field. Moreover, there is no “Hall drag”: the 
direction of the induced motion of charge carriers in the 
passive layer is expected to coincide with that of the driv¬ 
ing current. These predictions contradict numerous ex¬ 
periments [see, e.g., Eisenstein and MacDonald (2004); 
Einck et al. (2010); Lilly et al. (1998); Muraki et al. 
(2004); Nandi et al. (2012); and Rubel et al. (1997b)] 
showing that Coulomb drag is not only sensitive to mag¬ 
netic field, but in fact the drag resistivity can be greatly 
enhanced, once the held is applied. 

In single-layer measurements, magnetoresistance is 
usually associated with either (i) multi-band systems, or 
(ii) quantum effects. A close analog of the former can be 
found in graphene-based systems, see Gorbachev et al. 
(2012), Titov et al. (2013a) and Sec. IV. The latter ef¬ 
fects are manifest in strong, quantizing magnetic Helds 
leading to emergence of a qualitatively different behavior 
(Eisenstein and MacDonald, 2004; Girvin and MacDon¬ 
ald, 1997) discussed in Sec. VH. 

The situation somewhat simplihes if the Held is tuned 
close enough to the point where the Landau level Hll- 
ing factor in the two layers is about = 1/2. In this 
case, the many-body state in each layer can be viewed 
as a Eermi liquid of composite fermions (Halperin et al., 
1993). In some sense composite fermions are similar to 
usual electronic quasiparticle in metals. Long-range, in¬ 
terlayer interaction between such excitations could lead 
either to a “single-particle” drag effect (Kim and Millis, 
1999; Sakhi, 1997; Ussishkin and Stern, 1997, 1998), or 
to novel correlated states, see Sec. VH. Alternative ap¬ 
proaches include magnetodrag due to electron-phonon in¬ 
teraction (Badalyan and Kim, 2003), semiclassical theory 
(Brener and Metzner, 2005), diagrammatic theory in high 
Landau levels (Bpnsager et al., 1996, 1997; Gornyi et al., 
2004; von Oppen et al., 2001), self-consistent Hartree ap¬ 
proximation (Tso et al., 1998), and the effect of magneto- 
plasmons (Khaetskii and Nazarov, 1999; Manolescu and 
Tanatar, 2002). 
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1. Hall drag in weak (classical) magnetic field 


Recall that the standard single-band Drude theory (2) 
does not allow for any dependence of the drag resistiv¬ 
ity on the magnetic held and in particular predicts zero 
Hall drag, see Eq. (3d). Same conclusion can be reached 
using the diagrammatic perturbation theory (Kamenev 
and Oreg, 1995). This result is justihed by the assump¬ 
tion of energy-independent impurity scattering time r. 
Lifting this assumption (Hu, 1997), one can show that a 
weak Hall drag signal may appear 


pg oc sT^, 


dT{e) Ep 

de t{Ep) ' 


(49) 


As argued in Hu (1997), this effect is hard to observe 
in conventional semiconductor heterostructures where in¬ 
tralayer relaxation processes are dominated by electron- 
electron interaction: in this case the nonequilibrium dis¬ 
tribution function quickly relaxes to a drifted Fermi- 
Dirac distribution and hence the impurity scattering time 
is effectively almost independent of energy, i.e. s <C 1. 

Hall drag in weak magnetic helds was studied in Pa¬ 
tel et al. (1997). The experimental device comprised two 
180A-wide quantum wells separated by 100A and exhib¬ 
ited measurable tunneling between the layers, contrary 
to the assumptions of Hu (1997). Hall drag in graphene 
(Titov et al., 2013a) was attributed to a different mech¬ 
anism, see Sec. IV. Other observations of Hall drag were 
performed in the quantum Hall regime (see Sec. VH), 
where the effect is much stronger (von Oppen et al., 2001) 
than Eq. (49). 


2. Coulomb drag of composite fermions 


All of the previous discussion was based on the un¬ 
derlying physical picture of weakly interacting fermions. 
Typically, this picture becomes invalid in a strong, quan¬ 
tizing magnetic field. The only exception to this state¬ 
ment is the peculiar state at the half-filled Landau level. 
This state can be described as a Fermi liquid of composite 
fermions (Halperin et al., 1993). Each composite fermion 
is an electron with two attached flux quanta (Jain, 1989), 
that interacts with the others both electrostatically and 
by means of a Chern-Simons interaction. 

Composite fermions can be characterized by linear re¬ 
sponse functions similar to those of electrons. In partic¬ 
ular, their respective single-layer resistivities related to 
each other by (Halperin et al., 1993) 

2h f 0 

Pel = Pcf +■ (50) 

If one is interested in the relation between conductivities 
of the electrons and composite fermions, then one has 


to invert the resistivity matrices in Eq. (50). Clearly, 
the electronic conductivity is not identical to that of the 
composite fermions. 

Extending Eq. (50) to the case of a double-layer sys¬ 
tem, one obtains a similar relation for the 4x4 resistivity 
matrices (Ussishkin and Stern, 1997). If interlayer inter¬ 
action is weak enough, so that composite fermions in a 
given layer are not sensitive to the Chern-Simons field of 
the other layer, then similarly to Eq. (50), longitudinal 
resistivities of the electrons and composite fermions are 
the same and hence 


Pd = Pd- ( 51 ) 

Again, conductivities (in particular, drag conductivities) 
of electrons and composite fermions are not equivalent. 

The quantity measured in drag experiments is the elec¬ 
tronic drag resistivity pf), Eq. (13). Given the equality 
(51), one can calculate either or p|(. The former ap¬ 
proach was developed in Kim and Millis (1999), while 
the latter was considered in Ussishkin and Stern (1997). 
Both calculations are based on the standard lowest-order 
perturbation theory and yield similar results (albeit with 
a rather different interpretation^^). The calculation of 
Kim and Millis (1999) consists evaluating Eq. (15) for 
composite fermions and using the correspondence (51). 
Alternatively (Ussishkin and Stern, 1997), one can treat 
the problem in purely electronic terms assuming that in¬ 
terlayer interaction is dominated by the direct Coulomb 
coupling [the assumption which justifies Eq. (51)]. At 
the same time, single-layer electronic response functions 
(such as ImH^) can be calculated within the composite- 
fermion approach of Halperin et al. (1993). 

Within RPA (including the response of composite 
fermions to the external. Coulomb, and Chern-Simons 
potentials) and in the limit q kp, uj ^ vpq, the elec¬ 
tronic density-density response function (the polarization 
operator) is given by (Halperin et al., 1993) 

^ dp q^ — 8TTiu}kp{dn/diE)’ ^ ^ 

where dn/dp is the thermodynamic compressibility of the 
V = 1/2 state. At large momenta, ImH”^ ex q~^; conse¬ 
quently, the momentum integration in Eq. (15) is domi¬ 
nated by the region q « kpiT/ToY^^ [i.e., determined by 
poles of the interlayer interaction, rather than Eq. (6)]. 
As a result, the temperature dependence of the drag re¬ 
sistivity is weaker than in the absence of magnetic field 


The subquadratic temperature dependence (53) of the drag re¬ 
sistivity at = 1/2 was interpreted in Kim and Millis (1999) as a 
signature of the non-Fermi-Liquid nature of composite fermions. 
In particular it was related to the similar power law in the self¬ 
energy of the composite fermions leading to the ui ^ scaling 
of the typical frequencies. 
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FIG. 12 (Color online) Coulomb drag measurements at 
V = 1/2. Left panel: the top plot shows the experimental 
Pd{T) (solid line) compared to the theory of Stern and Us- 
sishkin (1997); the bottom plot shows the field dependence. 
The inset shows pn at B = 11.45T for two values of the driv¬ 
ing current. [Reproduced from Lilly et al. (1998).] Right 
panel: (A) circuit schematic; (B) Pd{T) (dots) vs Eq. (53); 
(C) Pd{B) for different temperatures, T = 0.05 — 5.6K. The 
vertical line corresponds to the B field at which the points 
plotted in panel (B) were measured. [Reproduced from Price 
et al. (2010).] 


(Stern and Ussishkin, 1997; Ussishkin and Stern, 1997) 

PD = 0.825(/i/e2)(r/To)4/3^ (53) 

where the characteristic temperature depends on the car¬ 
rier density n, interlayer spacing d, dielectric constant e, 
and thermodynamic compressibility 

e /fin\ 

2TTe^d \dpj 

The same temperature dependence was reported in Kim 
and Millis (1999) and Sakhi (1997). 

For realistic parameter values similar to those of the 
experiment of Lilly et al. (1998), the drag resistivity (53) 
is much larger than the zero-field result (21). This fact 
is associated with the smallness of the typical momenta 
involved in the interlayer scattering processes and slow 
relaxation of density fluctuations in the u = 1/2 state. 

The role of disorder in the drag effect in the z/ = 1/2 
state was considered in Stern and Ussishkin (1997). In 
the diffusive regime, the polarization operator is given by 
the standard form (30) and hence the drag resistivity is 
given by Eq. (34), albeit with a different diffusion coeffi¬ 
cient than the same system would have in the absence of 
magnetic field. The result is much larger than at B = 0. 
In the diffusive regime, this follows from the observation 
that the longitudinal conductivity (or the diffusive con¬ 
stant, which encodes all microscopic details) at v = 1/2 
is much smaller than at B = 0. 

Although the above theory is qualitatively similar to 
the experimental observations (Lilly et ai, 1998) (e.g. 


Tn = 


Tie^nd 


drag at = 1/2 is much larger that at i? = 0; the tem¬ 
perature dependence in clean samples is subquadratic), 
theoretical calculations significantly underestimate the 
overall value of po as compared to the experiment of 
Lilly et al. (1998). Yang (1998) suggested, that the rea¬ 
son for the discrepancy is that the interlayer separation in 
the samples of Lilly et al. (1998) was close to the critical 
value, where the system forms an incompressible inter¬ 
layer state (for a detailed discussion of correlated states 
see Sec. VII). An alternative suggestion by Ussishkin 
and Stern (1998) attributes the unexplained features of 
the experiment (including the extrapolated nonvanish¬ 
ing drag at r = 0) to pairing fluctuations of composite 
fermions. The two scenarios can be distinguished by mea¬ 
suring Hall drag, which vanishes in the latter theory. A 
later experiment (Price et al., 2010) reported both the 
magnitude and temperature dependence of Eq. (53) to 
be in good agreement with the measured data. 


III. MESOSCOPIC FLUCTUATIONS OF COULOMB 
DRAG 

Universal conductance fluctuations (Altshuler et al., 
1991) is a quantum interference effect which is a manifes¬ 
tation of the wave nature of electrons. As the same elec¬ 
trons are responsible for Coulomb drag, it is natural to 
expect that the drag resistivity also exhibits mesoscopic 
fluctuations. The drag fluctuations were first studied the¬ 
oretically in Narozhny and Aleiner (2000) and Mortensen 
et al. (2001, 2002a) and then observed experimentally 
(Kim et al., 2011; Price et al., 2008, 2007, 2010). 

In a disordered system, it is impossible to track each 
individual impurity and one uses a statistical approach. 
Impurities are described by a distribution function and 
each physical quantity is treated as being random. Ob¬ 
servables correspond to average values of the random 
physical quantities with respect to the distribution of im¬ 
purities. If a system is large enough, it can be viewed 
as a combination of smaller parts, which become sta¬ 
tistically independent if they are separated by distances 
larger than any relaxation length. Then instead of av¬ 
eraging over a statistical ensemble, one can average over 
the volume of one large system. 

In the problem of electronic transport, averaging over 
the system volume can be understood as averaging over 
all possible paths that an electron can take moving be¬ 
tween points A and B (Aleiner et al., 1999; Altshuler 
and Aronov, 1985; Altshuler et ai, 1991). Such paths 
can cover all of the system volume and thus experience 
all possible local impurity configurations, making average 
over the system equivalent to ensemble averaging. 

Consider two paths between the points A and B. The 
total transmission probability is determined by the abso¬ 
lute value of the sum of the corresponding quantum am- 
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plitudes (Altshuler and Lee, 1988; Altshuler et al, 1991): 

W = |Ai + A 2 P = |Aip + IA 2 P + 2|Ai|IA 2 I cos((pi — (P 2 ), 

where (pi( 2 ) are the quantum-mechanical phases that an 
electron accumulates along the paths. Typically, the 
phases v?i( 2 ) are random (or incoherent). As a result, the 
interference term vanishes upon averaging over all pos¬ 
sible paths (or impurity configurations), leading to the 
semiclassical sum of transition probabilities 

(cOs((/?i — (P 2 )} = 0 ^ W = Wi + W 2 , Wi( 2 ) = 1 ^ 1 ( 2 ) P- 

Random quantities can be characterized not only by 
their average value, but also by higher moments of their 
statistical distribution, which may be sensitive to the in¬ 
terference term even if the phases ^pi( 2 ) are still random. 
Indeed, fluctuations of the transition probability 

{[W - {W)f) = 4ITiIT2(cos2((^i -ip2)) = 2 W 1 W 2 . 

are completely determined by the interference term. 

Fluctuations of the transmission probability result in 
fluctuations of transport coefficients. The effect of such 
fluctuations can be observed only in small enough sam¬ 
ples (Altshuler et al., 1991). Indeed, in order justify the 
concept of the phase associated with a given electronic 
path, the length of the path should be less than a typical 
inelastic relaxation length otherwise coherence would 
be lost before the electron reaches point B (Anderson 
et al., 1979). At the same time, the path length should 
be larger than the mean free path in the system (other¬ 
wise electron motion along the path would be determin¬ 
istic). Therefore, typical paths (and hence, the sample 
sizes) should be characterized by intermediate lengths L 

£ < A < 

Fluctuations observed at such length scales are known as 
“mesoscopic fluctuations” (Altshuler et al, 1991). 


A. Drag fluctuations in conventional diffusive samples 

Mesoscopic fluctuations of the usual conductance (Alt¬ 
shuler et al., 1991) are known as the “universal conduc¬ 
tance fluctuations” (UCF). The universality is manifest 
when T Et, where Et is the Thouless energy of the 
sample (i.e. in small samples or at low temperatures; in 
the diffusive regime, Et = D j = g/ {2 'kvLF'), with g 


In special cases of coherent paths (for instance, time-reversed 
paths) the phase difference is exactly zero. Then the interfer¬ 
ence term does not vanish and leads to quantum corrections to 
semiclassical transport properties, such as the weak localization 
correction (Altshuler and Aronov, 1985; Altshuler et al., 1991). 


being the dimensionless conductance). In this case, the 
fluctuations are characterized by the universal value 




1 


(54) 


where G = ge^/h is the conductance of the system. The 
latter equality emphasizes the fact that the dimensionless 
conductance is a function of the system size. 

In larger samples, (JG^) is a function of temperature 
and the sample size. Arguments leading to Eq. (54) are 
valid only for coherent samples (Altshuler et al., 1991). 
At larger length scales, L ^ the coherence is lost, and 
the disorder averaging should be performed by dividing 
the sample into patches of the size L,p. Individual self- 
coherent patches (54) can be combined as a network of 
random conductors. This yields (in dimension d) 


{6G^{L))c^{SG\L^))iLjLf. (55) 


The patches of the size L,p remain self-coherent as long 
as T ^ Et{L,p). At higher temperatures, thermal aver¬ 
aging should be performed up to energies of order T, 
suppressing the conductance fluctuations 

((5G2(A^)) cs {e^/hfET{L^)/T. 


The conductance fluctuations of the sample become 

{5G^[L-,T>ET{L^)])^{e^lhf{L,,/LYh/{TT^), (56) 

where = hEt^"^{L,p) is the dephasing time (Altshuler 
et al, 1980). 

The fluctuations (56) are only observable in mesoscopic 
samples. Assuming the samples to be “metallic”, 5 1, 

the UCF (54) yield only a small correction to the average 
value of conductance. For example, in the experiment of 
Price et al. (2007) the single-layer resistance fluctuates 
by about 200mn around the average of about 50017. 

Now, we apply the above arguments to Coulomb drag 
(Narozhny and Aleiner, 2000). The drag conductivity 
depends on (i) the phase space available to electron- 
hole excitations; (ii) matrix elements of the interlayer 
interaction; and (iii) electron-hole asymmetry, expressed 
through the energy dependence of the density of states 
(or, the density dependence of the single-layer Drude con¬ 
ductivity). This can be schematically summarized by 




/ phase 
I volume 


X 


/ matrix 
I element 


(57) 


The average drag conductivity [cf. Eqs. (34) and (35)] 
can then be understood (up to the logarithmic factor) by 
estimating the phase volume by T^, the matrix element 
by (coming from static screening), and the factor 

of the electron-hole asymmetry by Ep"^. 
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Fluctuations of the drag conductivity can also be esti¬ 
mated with the help of Eq. (57). Consider first the low¬ 
est temperatures T <C Et, where the sample is effectively 
zero-dimensional (OD). The phase space is then only lim¬ 
ited by temperature, yielding the usual factor of T^. The 
factor of the electron-hole asymmetry in Eq. (57) is a 
random quantity with the typical value ~ E^^, since the 
Thouless energy is the typical scale of mesoscopic effects. 
Interaction matrix elements in OD are independent of en¬ 
ergy (Aleiner et ai, 1999); fluctuations are determined 
by off-diagonal elements that contain a small factor of 
g~^. As a result, one finds the variance of the drag con¬ 
ductivity that strongly exceeds the average 


6a D 


e2 7^2 ^ (Sal) 

h g'^E^ ’ (ao) 




(58) 


The quadratic temperature dependence of the variance 
of the drag conductivity (Narozhny and Aleiner, 2000) 
for mesoscopic samples {L L^, T Et) was also ob¬ 
tained in the context of quantum circuits [see Sec. V and 
Levchenko and Kamenev (2008a)] and within the random 
matrix theory (Mortensen et ai, 2001, 2002b). 

In order to extend the OD argument to larger sam¬ 
ples, L 3> Lip, we again divide the system into patches of 
the size Lp. Since the patches are largely uncorrelated 
(due to the loss of phase coherence), they can be com¬ 
bined as a network of random conductors, see Eq. (55). 
Each patch can be analyzed similarly to the OD case. 
However, now the interaction matrix elements become 
energy-dependent on the scales larger than Et, decreas¬ 
ing with the transmitted energy w as |Mp ^ a;“^. Thus 
the energy transfer is limited by the Thouless energy of 
the patch uj ~ Et{Lp) = rather than temperature. 
As a result, the phase space is limited by Tt~^, rather 
than the usual . The fluctuations of the density of 
states (which determine the factor of electron-hole asym¬ 
metry) should now be calculated on the scale of tempera¬ 
ture rather than the Thouless energy. This suppresses the 
fluctuations in each layer by the factor of y/EriEpY/T. 
Combining the above estimates, we find 


e2 Tr-i Et{Lp) 
h g^E^Lp) 



(59) 


which is T-independent, in contrast to the OD result (58). 

Final the Coulomb drag fluctuations in 2D samples can 
be estimated by combining Eqs. (55) and (59): 


{5al{L)) 





Et{L)tp 


oc 


1 

T' 


(60) 


The temperature dependence of the fluctuations (60) is 
contained in the dephasing time Tp ^ g/T (Altshuler and 
Aronov, 1985). At high enough temperatures, T T*, 
the fluctuations are small [the average value oi a d (35) 
is representative], but for T <^T* fluctuations dominate. 




FIG. 13 (Color online) Left: Qualitative picture of the typical 
measured drag signal. At high enough temperatures T > T* 
the average drag conductivity (35) is representative, an oc T^, 
with positive coefficient. Below T*, fluctuations dominate and 
the sign of the measured signal becomes random, i.e. depen¬ 
dent on a particular configuration of disorder. For T < T* the 
temperature dependence weakens to ctd oc Ify/T. At very 
low (most likely, experimentally inaccessible) temperatures 
T < Et T*, the quadratic temperature dependence is re¬ 
stored, but with a random coefficient, as fluctuations in the 
effectively OD system are much stronger than the average, see 
Eq. (58). Right: Measured drag resistance as a function of 
carrier density in the passive layer for T — 1,0.4,0.24K. In¬ 
set: the temperature dependence of the same data for the 
two values of n 2 denoted by vertical dotted lines in the main 
plot. The line indicates the dependence. [From Price et al. 
(2007). Reprinted with permission from AAAS.] 


see Fig. 13. The crossover temperature T* can be found 
by setting the relative fluctuation to unity 

T* ~ Epig^nL^)-^/^. 

The fluctuation-dominated regime is characterized by 
typical values of a t> determined by Eq. (60) rather than 
the average. In particular, the temperature dependence 
of the measured drag conductivity in this regime appears 
almost saturating as ao oc l/VT. The value of the pref¬ 
actor in this expression is sample-dependent and has a 
random sign. If temperature is decreased further, then 
eventually (although probably only in theory) one may 
reach the regime where T < Et- Then the sample will 
become effectively zero-dimensional and the quadratic 
temperature dependence an oc will be restored. In 
this regime of lowest temperatures, fluctuations greatly 
exceed the average [see Eq. (58)] and therefore the coef¬ 
ficient in the quadratic temperature dependence will be 
random (with random sign). The temperature depen¬ 
dence of a typical drag signal is sketched in the left panel 
of Fig. 13 (cf. the inset in the right panel of Fig. 13; see 
also the inset in the lower left panel of Fig. 12). 

The above qualitative picture is in full agreement with 
microscopic calculations (Narozhny and Aleiner, 2000). 
The average square of the drag conductivity has the form 
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e 7 
IStt^ 



In xd 
g^{xd)^ 


(61b) 


where 7 = 1.0086. Comparing Eq. (61) with the average 
drag conductivity in the diffusive regime (35), one finds 
the cross-over temperature T* = i?i?( 167 rg^uL^)“^/^. 

For heterostructures used in Gramila et al. (1991) and 
Lilly et al. (1998), the value of T* can be estimated as 
T* « 0.2K, which is below the temperature range of these 
experiments. Hence, the average drag coefficients (34) 
and (35) were sufficient to account for the observed effect 
with no trace of the random sign predicted by Eqs. (61). 

The result (61) is valid for homogeneous 2D diffusive 
samples in the absence of magnetic field. The random¬ 
ness (i.e. the sample to sample variation) of the sign 
of the effect should be contrasted with the determinis¬ 
tic sign change of the drag resistivity suggested for bi¬ 
layer systems with in-plane periodic potential modula¬ 
tion (Alkauskas et al., 2002). Drag signals of both signs 
have been observed in vertically integrated ID quantum 
wires (Laroche et al., 2011). While the observed effect 
has been argued (Biittiker and Sanchez, 2011) to have 
a mesoscopic origin (Mortensen et al., 2001) dominated 
by charge fluctuations (Levchenko and Kamenev, 2008a; 
Sanchez et al., 2010), the data appear to be not ran¬ 
dom, but reproducible. Very similar data were obtained 
in the subsequent experiment (Laroche et al., 2014) and 
interpreted with the help of the Luttinger Liquid theory 
(Pustilnik et al., 2003) (see Sec. VI). 

More recently, drag fluctuations were observed in dif¬ 
fusive graphene-based double-layer samples (Kim et al., 
2011; Kim and Tutuc, 2012), see Fig. 14. The temper¬ 
ature dependence Sajj cx [following from Eq. (61) 

and the assumption that the main phase-breaking mech¬ 
anism in the device is electron-electron scattering (Alt¬ 
shuler and Aronov, 1985)] appears to be in agreement 
with the experimental data. Other aspects of this exper¬ 
iments are specific to graphene. The fluctuations appear 
to be more pronounced in the vicinity of the charge neu¬ 
trality point. However, at the time of writing, a theory 
of drag fluctuations in graphene has not been developed. 
There is also no explanation for the most puzzling feature 
of the data reported in Kim and Tutuc (2012) showing 
an apparent violation of Onsager reciprocity as the drag 
fluctuations depend only on the charge density in the 
passive layer and not in the active layer. 

Closing this Section, we note that strong fluctuations 
of Coulomb drag ultimately follow from strong fluctua¬ 
tions of the nonlinear susceptibility. The fact that the 
fluctuations of the drag resistivity exceed the average is 
related to the overall smallness of the drag effect due 
to electron-hole symmetry. A related phenomenon is the 



FIG. 14 (Color online) Mesoscopic fluctuations of Coulomb 
drag in graphene. At low temperatures the fluctuations fully 
obscure the average drag. The curves are shifted for clarity; 
the horizontal dashed lines indicate OH for each curve. [Re¬ 
produced from Kim et al. (2011).] 


fluctuations of the electro-acoustic current®, which is also 
determined by the same nonlinear susceptibility. 


B. Giant fluctuations of Coulomb drag 

The predictions of the fluctuation theory (Narozhny 
and Aleiner, 2000) were put to the test in the dedicated 
experiment (Price et al., 2007). Both the UCF and drag 
fluctuations were measured in the same structure. The 
observed UCF have shown the usual behavior (Altshuler 
et al., 1991). A direct comparison of the correlation fields 
for the UCF and drag fluctuations confirmed that both 
effects depend on the same coherence length (Aleiner 
et al., 1999) and have the same quantum origin. Surpris¬ 
ingly, the observed giant drag fluctuations (Price et al., 
2007) greatly exceeded the original prediction (Narozhny 
and Aleiner, 2000). This discrepancy was attributed to 
the fact that the experiment was performed in the bal¬ 
listic regime (Narozhny et al., 2002; Zala et al., 2001). 

Let us remind the reader (see Sec. H) that the in drag 
measurements difference between “diffusive” and “ballis¬ 
tic” samples is in the relation between the mean-free path 
I and interlayer separation d. The latter sets the upper 
limit for the interlayer momentum transfer due to the ex¬ 
ponential decay of the Coulomb interaction (5). Thus, if 
the mean free path is small i d, then q d~^ 
and the effect is dominated by the diffusive motion of 
charge carriers. In “cleaner” samples with £ ^ d, 
both small q £~^ and large £~^ q d~^ momen¬ 
tum transfers are possible. The conventional statement, 
that in such samples Coulomb drag is dominated by bal- 
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listically moving carriers (Kamenev and Oreg, 1995), fol¬ 
lows from observing that processes with large momentum 
transfers yield a much larger drag resistivity Eq. (21) 
compared with the diffusive result (34), see Eq. (36). 

Coherence properties of electrons are also sensitive to 
the nature of their motion. The dephasing time is 
a manifestation of inelastic electron-electron scattering 
(Aleiner et al, 1999; Altshuler and Aronov, 1985). The 
conventional theory of interaction effects in electronic 
systems (Altshuler and Aronov, 1985) yields the follow¬ 
ing estimate for the dephasing time in diffusive systems 


r^^(rr < 1 ) - (rin 5 )/g. (62a) 

At higher temperatures, transport is dominated by pro¬ 
cesses with one or few successive impurity scatterings. 
In this “ballistic” regime (Zala et ai, 2001), the dephas¬ 
ing time exhibits somewhat stronger temperature depen¬ 
dence (Narozhny et ai, 2002) 

t-\Tt » 1) ~ (T^/Ef) Iii{2Ef/T). (62b) 

In Eqs. (62a) and (62b) the parameter distinguishing the 
diffusive and ballistic regimes is Tr which is independent 
of the inter-layer separation. This is to be expected since 
the theory leading to Eqs. (62a) and (62b) was devoted 
to two-dimensional systems and not bilayers. 

The effect of the external magnetic field on the single¬ 
layer conductance fluctuations analyzed in Price et al. 
(2007) demonstrates the expected crossover between the 
ballistic and diffusive results : 


-1 


oc 



Tt < 1 , 
Tt > 1 . 


(62c) 


The same sample where Coulomb drag is dominated by 
the ballistic motion of electrons with large interlayer mo¬ 
mentum transfers, <C <7 ^ d~^, may exhibit both the 
diffusive and ballistic behavior of single-layer transport 
properties, e.g. of the dephasing time (62). 

Similar crossover was observed also in the drag fluc¬ 
tuations that exhibited strikingly different temperature 
dependence at large and small Tr (Price et al., 2007): 


((5cr|,) OC 



Tt< 1, 
Tt> 1 . 


(63) 


The crossover temperature in Eq. (63) was found to be 
about the same as in Eq. (62). This coincidence raised 
the question of whether the large magnitude of the ob¬ 
served drag fluctuations and their unexpected temper¬ 
ature dependence (63) had the same origin that would 
involve large momentum transfers <C 9 «C [given 
that the small momentum transfers lead to Eq. (61)]. 

Interaction processes characterized by large momen¬ 
tum transfers q ^ £~^ involve two electrons separated 



FIG. 15 Drag fluctuations in ballistic samples. The lines rep¬ 
resent the asymptotic power laws, see Eq. (63). The inset 
shows the measured dephasing time. The lines in the inset 
represent the power laws from Eq. (62). [Adapted from Price 
et al. (2007). Reprinted with permission from AAAS.] 


by a distance that is smaller than the average impurity 
separation. Therefore the effect should be determined by 
local electron properties. Local properties, such as the 
local density of states, are known to exhibit mesoscopic 
fluctuations stronger than those of the global properties 
(responsible for drag fluctuations in the diffusive regime). 
In particular, fluctuations of the local density of states 
are given by (Lerner, 1988) 

( 64 ) 

9 ^ 

where Lt = \JD/T is the thermal length. 

Contribution of processes with large momentum trans¬ 
fers to drag fluctuations can be estimated using Eqs. (57) 
and (55). Electron-electron scattering can be described 
with the help of “ballistic” expressions discussed in 
Sec. II.B. As small angle scattering plays the dominant 
role, the matrix element of the interlayer interaction is 
proportional to the ratio of the mean-free path to the 
interlayer separation |M|^ ^ g~'^P' j(>c^d^'). As this in¬ 
teraction is static, the phase space is only limited by T. 
Assuming that the fluctuations of the nonlinear suscepti¬ 
bility are dominated by fluctuations of the local density 
of states (64), one finds (up to a logarithmic factor) 


where the thermal smearing was taken into account sim¬ 
ilarly to Eq. (59). This leads to the estimate for the drag 
fluctuations in the whole sample (Price et al., 2007) 




e4 £4 rj.2 7^2 
^2 52^4^8 EI { L ^) T2 


(X TV3. 


( 66 ) 


The result ( 66 ) contains two falsifiable predictions: (i) 
the magnitude and (ii) the temperature dependence of 
the drag fluctuations. In comparison with Eq. (61), the 
prefactor in Eq. (65) contains the large factor £‘^/d‘^ and 
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moreover, g'^ instead of in the denominator. Con¬ 
sequently, the drag fluctuations (66) are much stronger 
than the diffusive prediction. At the same time, using the 
temperature dependence of the dephasing time (62), one 
immediately recovers the measured temperature depen¬ 
dence of the drag coefficient (63). The crossover between 
the two temperature regimes in Eqs. (63) and (66) is il¬ 
lustrated in Fig. 15. 


C. Drag fluctuations at the half-filled Landau level 


constant can be taken at their average values (their fluc¬ 
tuations are much smaller than the averages). To es¬ 
timate fluctuations of da/dn, one can express the con¬ 
ductivity in terms of the response functions of composite 
fermions using Eq. (50). On average, the conductivity 
matrix of composite fermions is diagonal. Assuming the 
large dimensionless conductance of composite fermions, 
gcf 1, the electronic longitudinal conductivity is in¬ 
versely proportional to gcf, meaning smallness of the 
electronic dimensionless conductance 

g « l/{4gcf) « 1. (67) 


Mesoscopic fluctuations of Coulomb drag of composite 
fermions were studied theoretically in (Narozhny et al., 
2001) and experimentally in (Price et al., 2010). De¬ 
spite the significant increase in the magnitude of drag of 
composite fermions relative to that of normal electrons 
(Jdrger et al, 2000c; Lilly et al., 1998; Muraki et al, 2004; 
Zelakiewicz et al, 2000) the fluctuations of the drag re¬ 
sistivity can still exceed the average, resulting in an al¬ 
ternating sign of the measured drag resistivity. 

Qualitatively, one can estimate the fluctuation effects 
using Eq. (25). Similarly to the B = 0 case, drag fluctua¬ 
tions stem from the fluctuations of the nonlinear suscep¬ 
tibility. In the diffusive regime, (F) is given by Eq. (28) 
with the polarization operator having the standard form 
(30), although with a different diffusion constant (Stern 
and Ussishkin, 1997). In contrast to the B — 0 case, the 

= 1/2 state is characterized by a large Hall conduc¬ 
tivity. This leads to the nonlinear susceptibility being 
approximately orthogonal to the transferred momentum 
q [unlike Eq. (31)]. 

Fluctuations of F (and thus of the drag resistivity) 
result from mesoscopic fluctuations of da jdn. Other pa¬ 
rameters, such as the compressibility and the diffusion 


This is the reason one needs to perform calculations in the 
composite-fermion basis: the B = 0 theory of Sec. III.A 
is justified by the small parameter 1/g. 

Adapting the B = 0 theory to the case of composite 
fermions, Narozhny et al. (2001) found the fluctuations 
of the nonlinear susceptibility (28) of a coherent sample 
of size L in the = 1/2 state to be large 
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similarly to Eq. (58). This is already an observable con¬ 
clusion: in a fully coherent sample in the diffusive regime, 
fluctuations of the acoustoelectric current (determined by 
the same nonlinear susceptibility) are much larger that 
its average. The result (68) is justified as long as the 
thermal /T and phase breaking length 

scales of composite fermions are much larger than L. 

For larger samples, the global phase coherence is lost 
and one has to employ the averaging procedure described 
in Sec. III.A. The system can be divided into L/^ 
self-coherent patches of the size of the phase-breaking 
length of composite fermions . Summing up contri¬ 
butions of all patches according to Eq. (55), one finds 
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(69) 


where as « 0.2(32/97r) = 0.23 and the coefficients ai ,2 
are of order unity (Narozhny et at, 2001). 

The magnitude of the mesoscopic fluctuations depends 
on the precise source of phase breaking, but their temper¬ 
ature dependence is robust: all generic models of phase 
breaking in two dimensions (Altshuler and Aronov, 1985) 
lead to \|T^p oc T in the diffusive regime. In the v = 1/2 
state phase breaking comes from the quasi-elastic scat¬ 
tering of composite fermions off the thermal quasi-static 
fluctuations of the Chern-Simons magnetic field. This 
mechanism can be illustrated using a cartoon shown in 
Fig. 16. Consider a density fluctuation where the excess 


charges in the two layers have opposite signs. Such a 
fluctuation is accompanied by a random flux that inter¬ 
acts with the composite fermions leading to the loss of 
coherence. The energy of this fluctuation is of order T. It 
can also be estimated as the energy of a simple capacitor, 
27re^d/[£(L^-^)^] ~ T, where £ is the bnlk dielectric con¬ 
stant and is the typical size of the density fluctuation 
with the electron number of the order of unity creating 
the random flux of approximately $o. As a result, 

1/T^f c^gcfT/i^d). (70) 

Substituting the above estimate into Eq. (69), one finds 
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FIG. 16 (Color online) A cartoon illustration of the phase¬ 
breaking mechanism a.t 1 / — 1/2. A random flux in the system 
can be generated by charge density fluctuations with the op¬ 
posite signs of excess local charges in the two layers. 


(assuming xrd) 

^ 2 \ 2'Ke^d 

{Pd) - ^ TeL-^g^^j:' 

Using realistic parameters (Lilly et ai, 1998) (i.e. L ~ 
100/rm; d = SOOA; T = 0.6iL; R = SfcU/D leading to 
gcf ~ 8; and {pjj) = ISn/D), the magnitude of the drag 
fluctuations can be estimated as dpo ~ 0.312. For lower 
temperatures and smaller samples, the theory predicts 
stronger fluctuations (i.a. exceeding the average). Such 
strong fluctuations were observed in Price et al. (2010), 
albeit with a substantially larger magnitude that follows 
from the above estimate. 

The dephasing time due to above mechanism of quasi¬ 
elastic scattering of composite fermions on thermal fluc¬ 
tuations of the Chern-Simons field appears to be shorter 
than the temperature scale <C 1. This does not cre¬ 
ate any additional complication since most of the phase 
breaking results from scattering off the Chern-Simons 
field fluctuations whose dynamics (with characteristic 
frequency T/gcf) is very slow compared to but fast 
compared to the time of the experiment. Field fluctua¬ 
tions which are static on the scale of the experiment time 
affect the mesoscopic fluctuations only by affecting gcf. 
Field fluctuations that are faster than that scale make the 
potential landscape seen by the composite fermions time 
dependent, and lead to a suppression of the mesoscopic 
fluctuations by partial ensemble averaging. 

Consider the correlation function 

{pD{B)pD{B+5B))-{pDm{pD{B+5B)) = F, , 

with the field B near the ;/ = 1/2 value. An experimental 
study of the decay of this correlation function is a way 
to measure : the characteristic magnetic field of the 

decay is B* ~ ^ $ 0 - The decay of this correlator 

as a function of density 

{pD{n)pD{n + 5n)) - {pD{n)){pD{n + 5n)) = F 2 

also yields L^J: the characteristic density change n* at 
which it decays is expected to correspond to half of an 
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FIG. 17 (Color online) Mesoscopic fluctuations of Coulomb 
drag a.t v — 1/2, T = 50mK. Upper panel: comparison of Spn 
as a fnnction of the filling fraction v obtained by varying ei¬ 
ther the carrier density or magnetic field. The red curve shows 
a different measurement run demonstrating the reproducibil¬ 
ity of the fluctuations. Similarity of the periods of pn(n) 
and pd{B) is the proof of composite-fermion drag. Lower 
panel: autocorrelation function of the fluctuations shown in 
the upper panel. Squares represent Pd{B) and circles - poin). 
[Reproduced from Price et al. (2010).] 


electron in a phase coherent region, i.e. n* = ^ / 2 . 

This statement holds as long as the composite fermion 
cyclotron radius is much larger than its mean free path, 
i.e., for \v - 1 / 2 | < {2gcf)~^. 

The ratio of the above characteristic field B* to the 
characteristic density n* yields two flux quanta 


B* 

— - 2$o. 
n* 


(71) 


This should be contrasted to the zero field case, where 


FI* ^ 

^ ^ TT’ 
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B* ^ $0 
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At B = 0 the electrons do not carry any attached flux. 
Therefore the characteristic density n* corresponds to 
a change in the chemical potential of order t~^. Con¬ 
sequently, observation of the ratio (71) in a laboratory 
experiment serves as a verification of the concept of the 
flux attachment and the fact that charge carriers in the 
system are indeed composite fermions. 

In single-layer measurements of mesoscopic fluctua¬ 
tions, the ratio (71) has been reported in Kvon et al. 
(1997). In double-layer systems, mesoscopic fluctuations 
of Coulomb drag were investigated in Price et al. (2010), 
where it was shown, that the fluctuations of drag resis¬ 
tivity obtained either by varying of the magnetic field 
(with n = const) or by varying the carrier density (hold¬ 
ing B constant) exhibit the same characteristic scale (or 
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a “period”), if plotted as a function of the filling fac¬ 
tor V = nh/eB, see Fig. 17. The similarity of the two 
“periods” is equivalent to the ratio (71). 


IV. DRAG IN GRAPHENE-BASED DOUBLE-LAYER 
DEVICES 

The physical picture of frictional drag outlined in the 
preceding Sections is based on the following assumptions: 
(i) each of the layers is in a Fermi-liquid state, which at 
the very least means /ii( 2 ) 3> T; (ii) electron-electron in¬ 
teraction does not contribute to the intralayer transport 
scattering time; (iii) the interlayer Coulomb interaction 
is assumed to be weak enough, a = /vp ^ 1, such that 
Pd is determined by the lowest-order perturbation theory 
(Flensberg and Hu, 1994; Flensberg et ai, 1995; Jauho 
and Smith, 1993; Kamenev and Oreg, 1995; Zheng and 
MacDonald, 1993) leading to Eq. (15). 

The drag resistivity (21) depends on the carrier den¬ 
sity, temperature, and interlayer spacing. In semiconduc¬ 
tor devices, only temperature can be varied with relative 
ease, and even then only within the limit /ii( 2 ) S> T [with 
the notable exception of Pillarisetty et al. (2002)]. 

Lifting one or more of the above assumptions leads 
to signihcant changes in the drag effect. Recently drag 
measurements were performed in a system of two parallel 
graphene sheets (Gorbachev et ai, 2012; Kim et ai, 2011; 
Kim and Tutuc, 2012; Titov et ai, 2013a). This system 
offers much greater flexibility compared to prior experi¬ 
ments in semiconductor heterostructures. The graphene- 
based system allows one to scan a wide range of chemi¬ 
cal potentials (by electrostatically controlling carrier den¬ 
sity) from the Fermi-liquid regime to the charge neutral¬ 
ity (or Dirac) point pi — 0. Moreover, using hexagonal 
boron nitride as a substrate (Ponomarenko et al., 2011; 
Titov et al., 2013a), one can decrease disorder strength in 
the system and reach the regime, where transport proper¬ 
ties of the two layers are dominated by electron-electron 
interaction, t ^ Tee- In addition, modern technology al¬ 
lows for a controlled growth of boron nitride yielding de¬ 
vices with a relatively wide range of the interlayer separa¬ 
tions, which can be as low as d = Inm (corresponding to 
only three atomic layers!). While the experiments (Gor¬ 
bachev et ai, 2012; Kim et ai, 2011; Kim and Tutuc, 
2012; Titov et ai, 2013a) were performed at relatively low 
temperatures T < Vg/d (vg is the quasiparticle velocity in 
graphene), the range of temperatures available for these 
measurements (typically, 4 — 240K) is much wider than 
in earlier studies. In a parallel development. Coulomb 
drag measurements in graphene double ribbon structures 
were reported in Chen and Appenzeller (2013). 

In graphene one can reach parameter regimes, which 
were inaccessible in semiconductor samples, see Fig. 18: 
(i) near charge neutrality, the chemical potential may be¬ 
come smaller then temperature, Mi ( 2 ) ^ T] the electronic 
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FIG. 18 (Color online) Summary of the parameter regimes 
and the resulting drag coefficient in graphene for identi¬ 
cal layers, p rain{T/a,v/ d), and uncorrelated disorder. 
Bottom row (below the curve 2, <C /p)’. solutions 
to the quantum kinetic equation, see Sec. IV.B. Curve 
1 (r“^ = o? ^ jT) separates the two regimes in Eq. (104). 
Middle row (cP'T <C t~^ <C T): the region where the QKE 
approach overlaps with the perturbation theory (Narozhny 
et al., 2012). The third-order contribution p® = 0{a^) yield¬ 
ing nonzero drag at p = 0 is shown in red. Upper row 
(r”'^ ^ T): the diffusive regime, where saturates for 
T~^ S> Tjo?). [Reproduced from Schiitt et al. (2013).] 


system becomes nondegenerate; (ii) low-energy excita¬ 
tions in graphene are characterized by the linear Dirac- 
like dispersion; there is no Galilean invariance in the sys¬ 
tem and transport properties are strongly affected by 
electron-electron interaction (Katsnelson, 2012; Schiitt 
et al., 2011). Moreover, electrons interact by means of 
3D nonrelativistic Coulomb interaction, which breaks the 
Lorenz invariance of the Dirac Hamiltonian. 

Nondegenerate systems were considered in thee early 
work on frictional drag (Boiko and Sirenko, 1988; Ja- 
coboni and Price, 1988; Pogrebinskii, 1977; Price, 1983) 
in the context of semiconductors, where elementary ex¬ 
citations are typically modeled by quasiparticles with 
parabolic dispersion. In that case, electron-electron in¬ 
teraction plays a subleading role in single-layer transport 
(due to Galilean invariance). In contrast, in ultra-clean 
graphene near the Dirac point single-layer transport is 
dominated by electron-electron interaction (Narozhny 
et al., 2012; Schiitt et al, 2011). 

The low-temperature degenerate regime T ^ p can be 
achieved in doped graphene, which is expected to ex¬ 
hibit the same qualitative behavior as the semiconductor 
devices. Indeed, in ballistic samples and under the addi¬ 
tional assumption of the small screening length >cd 3> 1, 
one recovers (Narozhny et ai, 2012; Tse et ai, 2007) the 
standard expression for the drag resistivity (21), albeit 
with an extra factor N = 4 reflecting higher degeneracy 
of the single-particle spectrum in graphene. However, 
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this regime might be outside of the experimentally acces¬ 
sible parameter range of drag measurements in graphene- 
based samples (Gorbachev et al., 2012; Kim et al., 2011; 
Kim and Tutuc, 2012; Titov et al, 2013a). 

For weaker doping, the assumption of the small screen¬ 
ing length is invalid and the standard result (21) has 
to be modified (Narozhny et al., 2012). A perturbative 
treatment can still be developed as long as the transport 
properties of both layers are dominated by disorder (i.e., 
Tee t). If electron-electron interaction is weak enough 

min(l, r//Xi) <C 1, (72) 

then the drag conductivity is given by the standard ex¬ 
pression (25). Close to the Dirac point {fii T), this 
yields pu oc At intermediate densities (p ~ T), the 

drag coefficient reaches a maximum and then decays to¬ 
wards the asymptotic limit (21). This decay is character¬ 
ized by a long crossover from the logarithmic behavior at 
Pi > T to the standard result (21) that is only achieved 
for small screening lengths, Kd 1. Asa result, the den¬ 
sity dependence of poid-i ^ I’g/d) cannot be described by 
a power law. Partially due to this fact, several conflicting 
results for po have been reported in literature (Amorim 
and Peres, 2012; Carrega et al., 2012; Hwang et al., 2011; 
Katsnelson, 2011; Lux and Fritz, 2012; Peres et al., 2011; 
Song and Levitov, 2012, 2013; Tse et al., 2007). 


A. Perturbative regime in ballistic samples 

The perturbation theory is valid when transport prop¬ 
erties of the sample are dominated by potential disorder. 


such that T <C Tee, sec Eq. (72). In ballistic samples the 
mean-free path is large compared to the interlayer sep¬ 
aration £ ^ d. For experimentally relevant temperature 
range T < Vg/d, the latter condition is compatible with 
the more standard condition for ballistic transport in dis¬ 
ordered systems Tt ^ 1. The resulting parameter range 
occupies the middle row of the “phase diagram” shown 
in Fig. 18 between the line Tt ~ 1 and curve 2. 

Perturbative calculations in the ballistic regime can be 
performed using either the diagrammatic (see Fig. 2), or 
kinetic-equation approach (see Sec. II.B). In both cases, 
one arrives at the expression similar to (15), where the 
nonlinear susceptibility and screened interlayer interac¬ 
tion (and hence the polarization operator) have to be 
specified for Dirac fermions in graphene. 


1. Nonlinear susceptibility in graphene 


In contrast to the theory reviewed in Sec. II, here we 
are interested in a wide range of chemical potentials in¬ 
cluding the Dirac point /i = 0. The nonlinear susceptibil¬ 
ity and polarization operator in graphene for arbitrary p, 
and T were derived in Narozhny et al. (2012). Assuming 
the long, energy-independent impurity scattering time t 
and neglecting intralayer interaction, the nonlinear sus¬ 
ceptibility has the form 
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where [with W = uj/{2T), Q = Vgq/{2T), and x = p/T] 
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Under the same assumptions, the polarization operator is given by 
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The perturbative calculation amounts to using the polar¬ 
ization operator (76) to determine the effective interlayer 
interaction (7) and then evaluating the drag conductivity 
[cf. Eqs. (15) and (25)] (Tse et al, 2007) 


ocB 

(Jn = 
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16ttT ^ 
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(77) 


using the nonlinear susceptibility (73). The drag resis¬ 
tivity is then given by Eq. (13). For arbitrary /x and T 
this calculation has to be performed numerically (Lux 
and Fritz, 2012; Peres et al., 2011). At the same time, all 
qualitative features of the drag effect can be elucidated 
by using simple limiting values. 

The nonlinear susceptibility (73) decays exponentially 
for q ^ max(/i,T). In the vicinity of the Dirac point, 
T ^ /i, the integral that determines the function 
(/(IT, Q, x) cannot be evaluated in terms of elementary 
functions. It can be shown, however, that in this case the 
nonlinear susceptibility is proportional to /x/T (Narozhny 
et al, 2012) 


g{x < 1) oc /x/T, (78) 


which could be expected since drag is supposed to vanish 
- or, more precisely, to change sign - at the Dirac point. 

In the degenerate limit T ^ /x, the dimensionless func¬ 
tion g{W,Q,x) may be approximated by 


41T 

g{x > 1, jlTj < Q) « 1 - 
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sinhc 


(32 cosh Q + cosh X' 
(79) 

Furthermore, for g » vpq ^ OJ ^ T (or Q ^ IT) the 
nonlinear susceptibility becomes similar to the standard 
(Kamenev and Oreg, 1995) “Fermi-liquid” expression for 
the ballistic regime (17) 


g{x > 1, jlTj < Q) « 4a;/ (vgq), 


(80) 


where the extra factor of 4 corresponds to extra degen¬ 
eracy of Dirac fermions in graphene (Amorim and Peres, 
2012; Tse et al, 2007). 

The relation (18) between F and Imll is not satisfied in 
graphene. This follows from a direct comparison between 
their respective integral representations. In particular, 
the nonlinear susceptibility (73) vanishes at the Dirac 
point due to exact electron-hole symmetry, T{g = 0) = 0 
(Tse et al, 2007), while the polarization operator (76) 
remains finite, Imn(/x = 0) 7 ^ 0 (Schiitt et al, 2011). 

Similarly to the usual Lindhard function (Giuliani and 
Vignale, 2005; Lindhard, 1954), the polarization operator 
in doped graphene has the simple static limit 

n^ito = T = 0) = 2kF/{TrVg). (81) 


2. Lowest-order perturbation theory 


We now use the above approximations to find the limit¬ 
ing expressions for the drag resistivity in the perturbative 
regime (Narozhny et al, 2012). 

In the simplest limit Nag <C T, the perturbative ap¬ 
proach is justified automatically. In this case, the single¬ 
layer conductivity is determined by weak impurity scat¬ 
tering and has the form 


(To = e^Tr ho (g/T), 


(83a) 


where 
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In this limit screening is ineffective and for /x^, T <C Wg/d 
the interlayer spacing drops out of the problem. Then 
we may use the “bare” Coulomb potential (5), while the 
frequency and momentum integration in Eq. (77) are de¬ 
termined by the nonlinear susceptibility (73). 

Close to the double Dirac point, gi <C T, the nonlinear 
susceptibility can be approximated by Eq. (78), while the 
remaining integration is dominated by frequencies and 
momenta of order temperature, uj, Vgq ~ T, yielding a 
dimensionless coefficient. The resulting drag resistivity 
is given by 

PD{g^ < T) « l.Ala^{n/e^)igig 2 /T^). (84a) 


If only one of the layers is tuned close to the Dirac 
point, /xi <C T <C /X 2 , the drag conductivity (77) is in¬ 
dependent of the properties of the second layer, as the 
integration in Eq. (77) is still determined by the region 
LO,Vgq ~ T. The single-layer conductivity in the second 
layer is still determined by g 2 , see Eq. (83). As a result, 

Pnihi < T < /i2) « 5.8a^(fi/e^)(/ii//X2). (84b) 


In the opposite limit g^ T, the nonlinear susceptibil¬ 
ity is given by Eq. (79). Now the momentum integral in 
Eq. (77) is logarithmic and is dominated by large values 
of momentum Q ^ IT. The ratio of the hyperbolic func¬ 
tions in Eq. (79) is similar to the step function: it’s equal 
to unity for Q x and vanishes at larger values of mo¬ 
mentum Q ^ X. Therefore x effectively acts as the upper 
cut-off and the momentum integral can be approximated 
by a logarithm 
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At the Dirac point, the result is somewhat different 
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Consequently the drag coefficient is similar to the stan¬ 
dard results of Sec. II 
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This is to be expected, since at low temperatures T <C /i 
the phase-space argument yielding the dependence is 
justified and the electron-hole asymmetry determines the 
dependence on the chemical potential. The logarithmic 
factor is beyond such qualitative estimates (the result 
(84d) was calculated with logarithmic accuracy). 


a. Static screening for vanishing interaction strength 

For slightly stronger interaction (i.e. smaller dielectric 
permittivity of the insulating substrate) or slightly larger 
interlayer spacing the condition Nagii T, gii Vg/d 
breaks down and one needs to take into account static 
screening. Static screening corresponds to the approxi¬ 
mation (81) to the polarization operator. If the interac¬ 
tion strength is still small a —)■ 0, then the interaction 
can be described by 
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(85) 


where = 4 is due to spin and valley degeneracy. The 
additional constant in the denominator affects the loga¬ 
rithmic integral (84c). As the chemical potential is be¬ 
ing increased away from the Dirac point, the following 
regimes may be gradually achieved [here we discuss these 
regimes for the case of identical layers; generalization to 
the case of two inequivalent layers is straightforward]: 

(i) Na^ T fr Vg/d. This regime is identical to 
the above arguments leading to Eq. (84d). 

(ii) Nafi <^T Vg/d fL. If the chemical potential 
is increased beyond the inverse interlayer spacing, then 
the momentum integration in Eq. (84c) is cut off by Vg/d 
instead of qi. The logarithmic behavior of the drag con¬ 
ductivity will be modified and (Td no longer depends on 
the chemical potential 


(Td ^ ln[z;g/(Td)]. (86a) 


(iia) T Nafi ^ /r Vg/d. In this case one finds in¬ 
stead of Eq. (86a) 

(T£) ~ ln[I/(A^a)]. (86b) 


(hi) T «C Nafj, <C Vg/d fj,. Increasing the chemical 
potential further leads to the regime where the static 
screening can no longer be neglected. Now the lower in¬ 
tegration limit in Eq. (84c) is effectively given by the 
inverse screening length rather than the frequency. The 
upper limit is still determined by the interlayer spacing. 
Therefore the drag conductivity again depends logarith¬ 
mically on the chemical potential (Katsnelson, 2011) 

cr£) ^ ln[ng/(A^a^d)], (86c) 


but now this is a decreasing function, indicating the exis¬ 
tence of the absolute maximum of the drag conductivity 
as a function of the chemical potential. 



FIG. 19 (Color online) Sketch of the drag conductivity (in 
the units of a^e^r^) as a function of the chemical poten¬ 
tial illustrating the results (86). The blue line shows the 
quadratic dependence (84a) in the vicinity of the Dirac point. 
If T ^ NaVg/d, then the region (iia) should be replaced by 
(ii): the logarithmic dependence r^ln(l/Q) should be re¬ 
placed by T'^ ln[vg/{Td)], and the limits Vg/d and T/{Na) 
should be exchanged. [Reproduced from Narozhny et al. 
( 2012 ).] 


(iv) T Vg/d ^ iVof/i ^ ji. Finally, if the chemical 
potential is so large that the screening length becomes 
smaller than the interlayer spacing the momentum inte¬ 
gral in Eq. (84c) is no longer logarithmic. As the inte¬ 
gration is now dominated by momenta large compared to 
T, the nonlinear susceptibility may be approximated by 
Eq. (80), leading to the standard Fermi-liquid result 


(Td 


C(3) 
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XT = 'fakp, 


(86d) 


which differs from that of Kamenev and Oreg (1995) only 
by the factor reflecting valley degeneracy in graphene 
(Amorim and Peres, 2012; Katsnelson, 2011; Tse et al., 
2007). The above results are illustrated in Fig 19. 


b. Static screening for intermediate interaction strength 

The results (84) and (86) rely on the interaction weak¬ 
ness. For stronger interaction, Na >1, (i) the approxi¬ 
mation (85) might be unjustified and the full expression 
(7) for the interaction propagator should be used; (ii) 
the four regimes (86) may not exist, since it might hap¬ 
pen that T/(iVQ:) ^ T <Vgl[Nad) <^Vg/d. In this case, 
perturbative analysis can still be justified in the degener¬ 
ate regime, /j, ^ T, where there are two distinct regimes, 
(a) ^ Vg/d, and (b) fi Vg/d (Narozhny et al, 2012); 

the latter regime is usually identified with the Fermi- 
liquid result (86d). As the single-layer conductivity is 
still large and dominated by disorder, the condition (72) 
can be somewhat relaxed: 

Tee S> T ^ T~^ :s> /fi => a^Tr fi/T. (87) 

Proceeding under the assumptions of static screening and 
the ballistic regime (i.e., the dominant contribution to the 
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effect comes from large momenta Vgq > to), the result of and can be assumed independent of w. The frequency 
momentum integration is determined by the upper limit and momentum integrals factorize and neglecting WjQ 

under the square root in Eq. (79) one finds 


aD=a e T r /o -;a; — 
. 7 w, 
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where 


d{x) = Nax/2. (89) 


The results for weaker interaction, Eqs. (86), can be re¬ 
covered from Eq. (88) by neglecting terms proportional 
to d^ in the denominator [which corresponds to approx¬ 
imating the interlayer interaction (7) by Eq. (85)]. In 
the limit Vg/d, the function /q depends on a single 
parameter 


fo{x\ > 1) « /o(4Ad), 
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The function (90b) describes the crossover between the 
regimes (hi) and (iv) of Eqs. (86) (see Fig. 19). This can 
be seen by evaluating the integral in the two limits (here 
7 o « 0.577216 is the Euler’s constant) 

fo{y<- 1) « -32/3 (Iny-h 7 o-k 11/6), (91a) 


/o( 2/>1)«64C(3)2/-^ (91b) 

Numerically, this crossover spans a large interval of values 
of the chemical potential such that the Fermi-liquid result 
(86d) is practically unattainable in graphene-based drag 
measurements (Tutuc et al., 2009), see Fig. 20. 

In experiment one typically measures carrier density 
rather than the chemical potential (Gorbachev et al., 
2012; Kim et al., 2011; Kim and Tutuc, 2012; Titov et al., 
2013a). In graphene, the electron density is given by 


OO 



tanh — tanh 


2T 


— OO 


Using the asymptotic expressions 

1 \Y, 

n = —?r < 

TTvl \(41n2)AiT, 




(92a) 


(92b) 


one can obtain the qualitative dependence of po on n, 
see Table I. 


_ dQQ^e-^^Q _ sinh^ x 

[(Q -I- d(x))2 — d(x)2e“^^'3]^ (coshQ + cosha;)^ ’ 


parameter region drag coefficient 

/i < T po ~ nT-^ 

T p, Vg/d Pd ~ T'^n~^ \Tx{aN'n}^^d/vg) 

p'>Vg/d pD ^ pY^ ~ T^n~^d~'^ 


TABLE I Asymptotic expressions for the drag coefficient to 
the leading order of perturbation theory assuming “realistic” 
interaction strength aN > 1, identical layers ni = n 2 = n, 
and the experimentally relevant situation T < Vg I d. In the 
opposite regime T ^ Vg/d all results for pn should be divided 
by Td/vg (Lux and Fritz, 2012; Narozhny et al., 2012). 


The strongly doped, Fermi-liquid regime has attracted 
the most attention in literature. Most authors report the 
standard po ^ T^n~^d~'^ behavior (Amorim and Peres, 
2012; Carrega et al., 2012; Hwang et al., 2011; Katsnel- 
son, 2011; Narozhny et al., 2012; Tse et al., 2007) assum¬ 
ing the energy-independent impurity scattering time. 


3. Energy-dependent scattering time 


In graphene, the impurity scattering time strongly de¬ 
pends on the type of disorder and on energy (Katsnel- 
son, 2012). In particular, for Coulomb scatterers (Ando, 
2006; Cheianov and Fal’ko, 2006; Nomura and MacDon¬ 
ald, 2006, 2007) or strong short-range impurities (Ostro¬ 
vsky et al, 2006) 

T(e)=T-o|e|- (93a) 

For weak short-ranged disorder (Shon and Ando, 1998) 

r(e) = 7 /|e|. (93b) 

Moreover, quenched disorder in graphene experiences 
logarithmic renormalization (Aleiner and Efetov, 2006). 

Drag in the presence of Coulomb impurities was first 
considered in Hwang et al. (2011) and Peres et al. (2011). 
Both papers reported a stronger dependence of the drag 
coefficient on the carrier density and interlayer separa¬ 
tion, pd ~ T'^n~‘^d~^. This result was later disputed 
in Amorim et al. (2012); Carrega et al. (2012); and 
Narozhny et al. (2012). These authors showed that the 
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FIG. 20 (Color online) Results of the numerical evaluation 
of the drag coefficient. The squares represent the calcula¬ 
tion of Eq. (25) with the only approximation that the po¬ 
larization operator in the screened interlayer interaction (7) 
was evaluated in the absence of disorder. The red line cor¬ 
responds to the same calculation, with the polarization op¬ 
erator replaced by Eq. (81). The blue line was calculated 
with the approximate expression (90b), valid for fj, ';$> Vg/d. 
Left panel: a = 0.01, Td/vg = 0.1. Right panel: a = 0.3 and 
Td/vg = 0.2; log-log scale. The straight green line represents 
the Fermi-liquid result (86d). [Reproduced from Narozhny 
et al. (2012).] 

energy (or momentum) dependence of t is qualitatively 
irrelevant for the asymptotic behavior oi po- In the de¬ 
generate limit, microscopic calculations lead to the same 
results with T{p.) substituted in place of r. Close to the 
neutrality point, the drag coefficient acquires an addi¬ 
tional logarithmic factor 

PoiPr < T) ~ a^{h/e^)ipip2/T^)lnTT{T). (94) 


4. Plasmon contribution 


The dynamically screened interaction propagator con¬ 
tains plasmon poles, that may (see Section II.B) af¬ 
fect the resulting drag resistivity. Theoretically, plas- 
mons were studied in graphene monolayers (at T = 0) 
in Hwang and Das Sarma (2007); Schiitt et al. (2011); 
and Wunsch et al. (2006) and in double-layer graphene 
systems in Badalyan and Peeters (2012); Profumo et al. 
(2012); and Stauber and Gomez-Santos (2012). Renor¬ 
malization of the plasmon spectrum due to electron- 
electron interaction was considered in Abedinpour et al. 
(2011). Experimentally, plasmons were observed in 
graphene on Si02 substrate (Fei et al., 2011, 2012), 
graphene-insulator stacks (Yan et al., 2012a), and in 
graphene micro-ribbon arrays (Ju et al., 2011). Bound 
states of plasmons with charge carriers, the so-called 
“plasmarons” were observed in Bostwick et al. (2010) and 
Walter et al. (2011). Plasmons subjected to high mag¬ 
netic field were studied in Yan et al. (2012b). For reviews 
of graphene plasmonics see Grigorenko et al. (2012) and 
Luo et al. (2013). More recently, plasmonic excitations 
in Goulomb coupled Af-layer graphene structures were 
studied in Zhu et al. (2013). 


Within the above perturbative approach, i.e., Eq. (72), 
and for low enough temperatures, p,,T<^Vg/d, the 
plasmon contribution to drag is subleading (Narozhny 
et al., 2012). The plasmon pole appears in the region 
w > VgQ. Similarly to the situation in semiconductor de¬ 
vices (Sec. II.B), double-layer graphene systems admit an 
acoustic {ui ^ q) and an optical (w ^ plasmon modes 
(Das Sarma and Madhukar, 1981; Principi et al., 2012). 
In the case where the plasmon decay rate is small (as 
determined be either weak Goulomb interaction or weak 
disorder), one can use the (^-function approximation to 
the interlayer interaction propagator (Flensberg and Hu, 
1995). The corresponding contribution to the drag con¬ 
ductivity contains a small factor g^\D\'^ ~ for small 
momenta Vgq ^ aT (or for Vgq ^ T). If the energy 
dependence of the scattering time is taken into account, 
the small parameter is a^Tr, see Eq. (72). 

The above conclusion is illustrated in Figs. 20 show¬ 
ing a comparison between the full numerical evaluation 
of the perturbative drag coefficient using Eqs. (7), (73), 
(76), and (77) and the same calculation within the ap¬ 
proximation of static screening (85). Numerical modeling 
of experimental samples, see Figs. 21 and 22 below, in¬ 
cludes the contribution was automatically by using the 
dynamically screened interaction propagator (7). 

At the same time, quantitative description of ex¬ 
periments, especially in devices with wider interlayer 
spacing, might be significantly affected by such aspects 
as inhomogeneous dielectric background (Badalyan and 
Peeters, 2012; Garrega et al, 2012) and hybridization be¬ 
tween phonon and plasmon modes (Amorim et al., 2012). 
Plasmon-mediated drag between graphene wave-guides 
was suggested in Shylau and Jauho (2014). 


5. Drag between massless and massive fermions 


Graphene-based double-layer devices can be used to 
observe Goulomb drag between massless and massive par¬ 
ticles by coupling Dirac fermions in monolayer graphene 
to quasiparticles with parabolic dispersion in either bi¬ 
layer graphene (Scharf and Matos-Abiague, 2012) or a 
usual 2DEG (Principi et al., 2012; Scharf and Matos- 
Abiague, 2012). Experimental realizations were reported 
in Fisichella et al. (2014) and Gamucci et al. (2014). 

Theoretical analysis of Principi et al. (2012) and Scharf 
and Matos-Abiague (2012) is based on the standard ex¬ 
pression (77). Both works focus on the low-temperature, 
degenerate regime T p. As expected, in the case of 
strong screening Kd ^ I, both works reproduce the stan¬ 
dard result (21). For >cd <C I, the resulting drag coeffi¬ 
cient is still quadratic in temperature, but contains also 
a logarithmic factor reminiscent of Eqs. (86). Principi 
et al. (2012) report a d-independent drag in the spe¬ 
cial case kg = k 2 D (which implies a density mismatch 
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between the layers due to the difference in the degen¬ 
eracies of single-particle states). In the low-density limit 
n —)• 0 this yields po similarly to the results of 

Carrega et al. (2012), see also Table I. On the other 
hand, Scharf and Matos-Abiague (2012) report po oc 
in the limit d —>■ 0 and for Ug = 2n2D- Such discrepancies 
in the asymptotic behavior of pu may appear due to the 
complicated structure of the nonlinear susceptibility in 
graphene, see Eqs. ( 86 ) and Fig. 19. 

The predicted dependence is observed in experi¬ 
ment (Gamucci et ai, 2014) in the 10K< T < 40K range, 
although with the smaller magnitude. At higher temper¬ 
atures, a violation of Onsager reciprocity was observed. 
This was attributed to the interlayer current. Most inter¬ 
estingly, at lower temperatures T < lOK, the measured 
drag shows a marked upturn that may indicate a phase 
transition at Tc ^ 10 — lOOmK, see Sec. VII.B. 

The system of coupled Dirac and Schrodinger quasipar¬ 
ticles was also considered in Balram et al. (2014), where 
it was found that interspecies interaction plays a signifi¬ 
cant role in determining collective (plasmon) modes. 


6. Numerical evaluation of the drag coefficient 


The above discussion demonstrates that already at the 
perturbative level, the drag conductivity (77) exhibits 
multiple asymptotic dependencies. Consequently, virtu¬ 
ally every paper on the subject presents results of numer¬ 
ical evaluation of Eq. (77). In contrast to the earlier work 
on semiconductor devices (see Sec. II), most authors fo¬ 
cus on the density (or chemical potential) dependence 
rather than on the T-dependence. The overall shape 
of Pain) curves is qualitatively the same in all calcula¬ 
tions. At the Dirac point, drag vanishes, ^^(n = 0) = 0 
(this conclusion does not agree with the experiments of 
Gorbachev et al. (2012) and Titov et al. (2013a), see 
below). Deep in the degenerate (or low-temperature) 
regime, T p,, >cd 1, the drag coefficient reaches the 
standard decaying result (21). Therefore, for interme¬ 
diate densities there has to be a maximum, roughly at 
p T. The corresponding shape is shown in Fig. 20. 

Peres et al. (2011) presented detailed numerical cal¬ 
culations aimed at describing the experimental findings 
of Kim et al. (2011), see Fig. 21. This calculation in¬ 
cluded electrostatic modeling of the device (which in¬ 
cluded two insulators, Si 02 and AI 2 O 3 ), dynamically 
screened (within RPA) electron-electron interaction, and 
the realistic model of Coulomb impurities. For doped 
graphene layers, the results of the calculation show ex¬ 
cellent agreement with the data. 

Theoretical modeling of ultra-clean graphene double¬ 
layers (using boron nitride as a substrate as well as insu¬ 
lating spacer) based the theory of Narozhny et al. (2012) 
was performed by Titov et al. (2013b), see Fig. 22. In 



FIG. 21 (Color online) Results of the numerical evalua¬ 
tion (lines) of the drag coefficient and comparison with the 
data (symbols) of Kim et al. (2011). The interlayer spacing 
{d = 14nm) and dielectric constants of the insulating mate¬ 
rial were chosen to represent the experimental device. Inset: 
the relation of the carrier densities and gate voltage, obtained 
from the electrostatic model of the sample. [Reproduced from 
Peres et al. (2011).] 




FIG. 22 (Color online) Results of the numerical evaluation 
(Titov et al, 2013b) of the drag coefficient (red line) and the 
experimental data (green dots) (Ponomarenko, 2013). The 
values of T = 240K and d are taken from the experiment. 
The only fitting parameter is the energy-independent impu¬ 
rity scattering time (once the value of ag is chosen). The 
polarization operator was calculated at T = 240K and in the 
presence of disorder (in the ballistic regime). 


this calculation, the polarization operator was calculated 
at the experimental temperature in the presence of dis¬ 
order, in contrast to the T = 0, free-electron expression 
(76). The use of full, dynamically screened Coulomb in¬ 
teraction ensured that all plasmon-related features were 
taken into account automatically. Choosing realistic val¬ 
ues (Kozikov et al., 2010; Peres et al., 2011) for the effec¬ 
tive coupling constant, the only fitting parameter in this 
calculation was the impurity scattering time t, which was 
taken to be energy-independent similarly to the above 
discussion. Such calculation was also able to reproduce 
the data (Ponomarenko, 2013) in the doped regime. 
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The results shown in Figs. 21 and 22 confirm the ap¬ 
plicability of the perturbative approach to Coulomb drag 
in doped graphene. In contrast to similar calculations 
aimed at semiconductor devices (see Sec. II), these the¬ 
ories are able to reach quantitative agreement with the 
experimental data with the minimum of htting param¬ 
eters. This implies that frictional drag in graphene is 
dominated by Coulomb interaction, with phonons play¬ 
ing only a subleading role. The latter conclusion can 
be expected, given that electrons in graphene are phys¬ 
ically conhned to move in a two-dimensional plane and 
the rigidity of the crystal lattice (Katsnelson, 2012). 

B. Hydrodynamic regime 

The perturbation theory outlined in Sec. IV.A can be 
justified either in the case of weak interaction a ^ 1 or 
in the degenerate regime T, see Eq. (72). At the 
same time, the applicability condition (72) involves the 
impurity scattering time r: the perturbation theory fails 
if the system is “too clean”, or in other words, if elec¬ 
tronic transport is dominated not by disorder, but rather 
by electron-electron interaction. The latter affects trans¬ 
port properties of graphene due to the absence of Galilean 
invariance: the velocity of Dirac fermions v = v^pje is 
independent of the absolute value of the momentum and 
therefore total momentum conservation does not prevent 
velocity (or current) relaxation. As a result, electron- 
electron scattering is characterized by its own transport 
relaxation time, which may become smaller than the scat¬ 
tering time due to potential disorder. Tee r. 

“Ultra-clean” graphene double-layers were discussed in 
Narozhny et al. (2015); Schiitt et al. (2013); and Titov 
et al. (2013a) within the framework of the quantum ki¬ 
netic equation. In principle, solving the kinetic equation 
in a strongly interacting system is a formidable problem, 
that cannot be solved in general terms using presently 


available analytic methods. However in graphene, one 
can take advantage of the kinematic peculiarity specihc 
to Dirac fermions. Indeed, scattering of particles with 
almost collinear momenta is enhanced since the momen¬ 
tum and energy conservation laws coincide. This re¬ 
stricts kinematics of the Dirac fermions (Fritz et al., 
2008; Kashuba, 2008; Schiitt et al, 2011) and leads to 
the singularity in the collision integral. This singular¬ 
ity leads to the fast thermalization of particles within 
a given direction and allows one to derive macroscopic 
- or hydrodynamic - equations that generalize Eq. (2) 
for interacting Dirac fermions. In monolayer graphene, 
this approach was discussed in Foster and Aleiner (2009); 
Fritz et al. (2008); Kashuba (2008); Muller and Sachdev 
(2008); and Svintsov et al. (2012). An alternative macro¬ 
scopic approach to Coulomb drag in graphene^^ has been 
suggested in Song et al. (2013); Song and Levitov (2012, 
2013). 


1. Collinear scattering singularity 


Singular behavior of the collision integrals in the case of 
collinear scattering of the Dirac fermions (Arnold et al., 
2000; Fritz et al., 2008; Kashuba, 2008; Muller et al., 
2009; Schiitt et al., 2013) is central to the hydrodynamic 
approach to transport in graphene. 

The general form of the kinetic equation in layer i is 
given by Eq. (9) with the addition of the intralayer col¬ 
lision integral. If the system is weakly perturbed from 
equilibrium, then the distribution function can be writ¬ 
ten in the form (10). Weak deviations from equilibrium 
are associated with the smallness of the nonequilibrium 
correction h, allowing one to linearize the collision inte¬ 
grals (Lifshitz and Pitaevskii, 1981). The linearized form 
of the collision integrals is given by 


1 . 1 ', 2 ' 


= E «^12.1'2'/j?/r2’ 1 - /j.l' 1 - /S' [hj,!' + K2' - - KM 


(95a) 


where the function 


w^i.2;i',2' = |(1) 2|C/|l',2')S27r)^5(ei -|- £2 - ei' - € 2 ') S{ki + k2- k[- k^, 


(95b) 


determines the probability of scattering from states 1', 2' imation). Here (1, 2|[/|1', 2') is the interaction matrix 
into states 1,2 (within the Fermi Golden Rule approx- element. The indices i,j = 1,2 denote the two layers^^. 

In graphene, the interaction matrix elements are most 


The theory of Song et al. (2013); Song and Levitov (2012, 2013) 

relies on correlations of the disorder potential in the two layers, perturbative approach of Sec. II.B, the kinetic equation 


see Sec. IV.F. 


(9) contained only the interlayer collision integral. Therefore, 
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conveniently expressed in the basis of the eigenstates of 
the Dirac Hamiltonian |e, e„) labeled by their energy e 
and the unit vector = v/vg pointing in the direction 
of velocity (for a given spin and valley projection): 


, , 1 . 1 , _( 2 ) ( 2 ') 
|(l,2|C/|l',2')r = \U{q)\^ 1 + e e. 1 + e e. ^ 


(96) 

Here q ~ ki — is the transferred momentum and the 
two fractions are the “Dirac factors” (Katsnelson, 2012). 
Now one can separate quantities related to the initial and 
final states in the function Wi 2 ,i' 2 ' by using the identities 


(5(ei + £2 — Cl 


4) 


J dujS{ei 


s'l + Lj)S(e2 — ^2 — 


S(ki+k2 — k[ — k2) = J tfqS{ki — k[+q)S{k2 — k2 — q). 

The (5-functions yield e'l = Vg\ki + q\ and hence allow one 
to sum over the states 1' and 2' in the collision integral 
(95). Each of these sums result in a diverging factor^^. 

^ / (97) 

1 ' _^;2 

One can see that the divergence corresponds to collinear 
scattering by examining the angle ipkiq at the light cone: 

cos (pk^qiuJ = Vgq) = 1 ^ (fikiq ^ 0 {or tt). 

Hence, the argument of one of the above 5-functions van¬ 
ishes: e( = Cl + (jj. Similar conclusion follows for the mo¬ 
mentum fc 2 - Thus, all momenta are collinear. 

Physically, the divergence (97) represents the fact that 
for the linear spectrum the energy and momentum con¬ 
servation laws coincide. Consequently, any relaxation 
rate obtained by integrating the collision integral (95) 
over the state 2 will be logarithmically divergent. In 
order to regularize this divergence, one has to go be¬ 
yond the Golden-Rule approximation and take into ac¬ 
count renormalization of the spectrum (Abrikosov and 
Beneslavskii, 1971; Gonzalez et al, 1999; Son, 2007). 
This leads (Arnold et al, 2000; Fritz et al, 2008) to the 
appearance of a large factor |ln(a)| ^ 1 in generic re¬ 
laxation rates in graphene. In disordered graphene, this 
singularity is also cut off by disorder-induced broadening 
of the momentum-conservation delta-function (Narozhny 
et al, 2012). 


one could associate the states 1 and 2 with the active and passive 
layers and avoid extra layer indices. 

In Sec. IV.A, the nonlinear susceptibility (73) did not exhibit 
this divergence due to an accidental cancellation that is specific 
to the particular case of energy-independent impurity scattering 
time. In a more general situation the cancellation does not occur 
and as a result the rate contains an extra logarithmic factor, 
see Eq. (94). 


2. Macroscopic linear-response theory in graphene 


The collinear scattering singularity (97) allows for an 
approximate, yet nonperturbative solution of the kinetic 
equation in graphene (Fritz et al, 2008; Kashuba, 2008; 
Narozhny et al, 2015; Schiitt et al, 2013). The idea 
is to find zero modes of the collision integral and build 
macroscopic equations for the corresponding currents. 

The standard perturbative description of Coulomb 
drag is based on the energy-independent approximation 
for the nonequilibrium distribution function (10). This 
constant solution of the kinetic equation describes the 
single zero mode of the intralayer collision integral corre¬ 
sponding to conservation of the electric charge (or the 
number of particles). Macroscopic charge flow is de¬ 
scribed by the electric current (12). Integrating the ki¬ 
netic equation, one finds the macroscopic equation for j 
equivalent to the Drude theory, see Eqs. (2). Such solu¬ 
tion is justified by the condition (72), which means that 
the collision integral in the kinetic equation is dominated 
by disorder. 

In contrast, in “ultra-clean” graphene the collision in¬ 
tegral is dominated by Goulomb interaction. Using the 
collinear scattering singularity (i.e., for | ln(Q;)| ^ 1), one 
can neglect all but the zero modes of X [treated as an inte¬ 
gral operator acting on hi{e)]. In practice, this means re¬ 
taining only those terms in the power series of the distri¬ 
bution function hi which correspond to either zero modes 
of the collision integral, or to its eigenmodes with nondi- 
vergent eigenvalues. Fritz et al (2008); Kashuba (2008); 
and Schiitt et al (2013) have developed the following 
two-mode approximation 

hi = (oq ^ -I- aj*^e) V. (98) 

The vectors can be expressed in terms of the two 
macroscopic currents in graphene, the electric current 
(12) and the energy current 

Qi = '^ev6f^. (99) 

The appearance of inequivalent currents is the essential 
feature of graphene physics. In general, the collision inte¬ 
gral has three nondecaying eigenmodes; hence Narozhny 
et al (2015) have used the three-mode approximation: 

hi = (a^*^ -I- aWsign(e) -I- aj*^e) v. (100) 

The sign(e) mode is described by the imbalance current 
(Foster and Aleiner, 2009) 

Pt = ^sign{e)vSf^. ( 101 ) 

Integrating the kinetic equation with the help of either 
of the above approximations for the nonequilibrium dis¬ 
tribution function, one obtains macroscopic equations for 
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the currents j, Q, and P, that generalize the Ohm’s law 
for graphene Solutions of these equations yield linear 
response transport coefficients. Note, that this approach 
does not rely on the Kubo formula. In particular, the 
drag coefficient can be obtained without the use of the 
perturbative expressions (15) or (77). 

The simplest macroscopic equation describes the en¬ 
ergy current. In an infinite sample, where all quantities 
are homogeneous, the equation reads (Schiitt et ai, 2013) 

ev^gUE + (vl/c) [j xB]= Q/t, (102) 

where n is the carrier density in graphene (92). The col¬ 
lision integral does not contribute to Eq. (102) due to 
energy conservation. In the limit ^ 3> T, all currents are 
equivalent, such that j{fi T) « {e/ ^ T), and 
Eq. (102) becomes equivalent to the Ohm’s law (26). In 
this limit, the Galilean invariance is restored, all relax¬ 
ation rates due to electron-electron interaction vanish, 
and all three macroscopic equations become equivalent. 

At the charge neutrality point n = 0, the equation 
(102) yields (Vg/c) [j x B] = Q/t. This simple-looking 
relation illustrates all the essential qualitative features of 
linear response transport in graphene. Firstly, in the ab¬ 
sence of disorder, t —>■ oo, the equation becomes sense¬ 
less, at least when the system is subjected to external 
magnetic field. Physically, this means that in the ab¬ 
sence of disorder the assumption of the steady state in an 
infinite system becomes invalid: under external bias, the 
energy current increases indefinitely. Secondly, if the sys¬ 
tem is stabilized by disorder, but B — 0, then one finds 
<5 = 0. Finally, if the system is subjected to external 
magnetic field, the electric and energy currents are or¬ 
thogonal, j E Q. This leads to appearance of classical, 
positive magnetoresistance (Muller and Sachdev, 2008; 
Narozhny et ai, 2015) 

SR{B-, fj. = 0)(x ivgT/c^){B'^/T^), (103) 

as well as magnetodrag in graphene, see Sec. IV.D below. 
These results are in sharp contrast with the standard 
Drude theory, see Eqs. (3). 


3. Coulomb drag in weakly disordered graphene. 


Close to charge neutrality and in the presence of weak, 
uncorrelated disorder a^Tr ^ 1 (i.e. the 

drag resistivity in the absence of magnetic field was found 
in Schiitt et ai (2013) and has the form 


PoitJ-i < T) 


2.87 


a^piP2 


/rf -I- /i2 + 0.49T/ [o^t) 


(104) 


The full three-mode equations are too cumbersome to reproduce 
here, the interested reader is referred to Narozhny et al. (2015). 



FIG. 23 (Color online) Leading-order drag coefficient in the 
ballistic regime as a function of carrier densities (in units of 
10^^cm“^) for d — 9nm. Left: po sX T = 250K; the upper 
panel refers to ultra-clean graphene = 0.5K; the lower left 
panel shows the evolution of po with increasing disorder from 
T~^ = 0 to = 50K. Right: po for — 50K; the lower 
panel shows pu for T = 150, 200, 250, and 300K. [Reproduced 
from Schiitt et al. (2013).] 


As long as any (even infinitesimal) disorder is present, po 
vanishes at the double Dirac point poiPi = P 2 = 0) = 0, 
and grows sharply in its immediate vicinity, see Fig. 23. If 
only one of the layers is tuned to the Dirac point (median 
lines in Fig. 23), the drag resistivity always vanishes 

Pd{pi = 0 ,/ i 2 7 ^ 0 ) = pd{pi 7 ^ 0 ,/ r 2 = 0 ) = 0 . 

If one varies the carrier density in one of the layers 
through the Dirac point, then the drag resistivity changes 
sign. In the color maps in Fig. 23 this is represented by 
the color change between neighboring quadrants. The 
same sign pattern of the drag resistivity (in zero magnetic 
field) were observed in experiments of Gorbachev et al. 
(2012); Kim et al. (2011); and Kim and Tutuc (2012). 

At the double Dirac point in the absence of disorder, 
one finds (for B = 0) pd{pi = P 2 = 0) ~ a^gh/e^. This 
peculiar feature is shown in Fig. 23 by the black curve 
in the lower left panel. It is however unlikely that this 
result is relevant to the nonzero drag resistivity at the 
Dirac point observed in Gorbachev et al. (2012). A pos¬ 
sible explanation for this observation is provided by the 
higher-order effects (Schiitt et al., 2013). 

For intermediate disorder strength, a^T <^T, 

the applicability region of the hydrodynamic approach 
overlaps with that of the conventional perturbation the¬ 
ory reviewed in Sec. IV.A and one recovers perturbative 
results, see Fig. 18. 

Finally, let us stress the novel qualitative feature of 
the hydrodynamic approach: the electron-hole asymme¬ 
try does not play a definitive role in the drag effect. In¬ 
deed, the “drag rate” dominates the observable ef¬ 
fect only under the standard assumptions of the Fermi- 
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liquid behavior in the two layers. On the contrary, in the 
vicinity of the Dirac point in graphene, another scatter¬ 
ing process, the interplay of fast interlayer energy and 
current relaxation which is insensitive to the electron- 
hole asymmetry, becomes important. Further examples 
to such novel behavior are presented in Section IV.D. 

C. Diffusive regime 

In strongly disordered graphene samples or, equiva¬ 
lently, at the lowest temperatures, Tt I, the electron 
motion becomes diffusive. In this regime, the standard 
perturbative approach based on Eq. (77) is applicable. 
In particular, the polarization operator has the standard 
form (30). The nonlinear susceptibility can be found us¬ 
ing the argument leading to Eq. (28). In graphene close 
to the Dirac point, /i <C T <C the derivative of the 
longitudinal conductivity with respect to the carrier den¬ 
sity is independent of the precise nature of disorder and 
is given by (Schiitt et al., 2013) 

da jdn ~ 

In contrast to the theory reviewed in Sec. II, in graphene 
the Thomas-Fermi screening length is much longer than 
the interlayer spacing xrd ^ 1; hence one finds the fol¬ 
lowing expression for the drag resistivity 

Pd (/.tj < T < T~^) ^ {h/e^)a^p,ip2TT^, (105) 

vanishing at /i^ = 0 due to the electron-hole symmetry. 

In the degenerate regime, p,^ T, one recovers the 
usual quadratic temperature dependence of the drag re¬ 
sistivity. The behavior of pu in the diffusive regime is 
summarized in Fig. 18 (the upper row). The “Fermi- 
liquid” result (34) is only recovered in the academic limit 
of strong screening xd ^ 1. This regime is not shown in 
Fig. 18 since in graphene it can be reached only at the 
extreme values of the chemical potential, see Fig. 20. 

Calculations of the lowest-order drag resistivity in the 
diffusive regime are essentially the same in any system, 
see Sec. II.C. As shown in Fig. 18, the behavior of po 
at the lowest temperatures may be dominated by higher- 
order drag effects. 

D. Giant magneto-drag in graphene 

Although the effect of classical magnetoresistance in 
multi-band systems is well known in semiconductor 
physics (Seeger, 2002), the equivalent effect in Coulomb 
drag was only recently observed in graphene-based de¬ 
vices (Gorbachev et al., 2012; Titov et ai, 2013a). One 
of the reasons is that the majority of earlier drag mea¬ 
surements were performed in double-well semiconductor 
heterostructures. Then each of the layers is represented 




FIG. 24 (Color online) Left: Off-diagonal resistivity in 
magnetic field, measured in a graphene-based double-layer de¬ 
vice. The two graphene sheets are kept at “opposite” carrier 
densities ni = —n 2 = n and T = 150K. [Reprinted by permis¬ 
sion from Macmillan Publishers Ltd: Nature Physics, Gor¬ 
bachev et al. (2012)] Right: Mechanism of magnetodrag at 
charge neutrality. Upper panel: in an infinite system quasi¬ 
particle currents in the two layers (denoted by Pi) flow in 
the same direction, leading to positive p),^. Lower panel: in a 
thermally isolated system no net quasiparticle flow is possible; 
the quasiparticle currents in the two layers have opposite di¬ 
rections yielding negative p),^. [Reproduced from Titov et al. 
(2013a).] 

by a two-dimensional electron gas that is formed by elec¬ 
trons occupying the lowest level in the quantum well at 
the interface between two semiconductors in the device. 
In contrast in graphene, the conductance and valence 
bands touch at the Dirac point and as a result, both 
electrons and holes participate in transport phenomena 
at low doping. 

The experimental data on magnetodrag in graphene 
(Gorbachev et al., 2012) is shown in Fig. 24 There 
are two outstanding features in Fig. 24. At high carrier 
densities (or in the degenerate regime), the effect of the 
magnetic field is relatively weak. This observation is con¬ 
sistent with the expectation, that transport properties of 
doped graphene are dominated by one of the two bands 
(the contribution of the other being exponentially sup¬ 
pressed) . In the vicinity of the Dirac point both types of 
carriers contribute to transport. Moreover, the leading 
contribution to drag at zero field vanishes right at the 
neutrality point due to exact electron-hole symmetry of 
the Dirac spectrum. Once the magnetic field is applied, 
the system develops a drag signal which is no longer de¬ 
termined by electron-hole asymmetry. As a result, the 
drag resistivity near the Dirac point in the presence of 
weak magnetic field is much higher than the maximum 


Notice, that Gorbachev et al. (2012) adopted an alternative def¬ 
inition of the drag resistivity pj?„, = E 2 xljix. Therefore in this 
Section we will discuss the off-diagonal resistivity = E 2 x/jix 
rather that po that is defined in the rest of the paper with the 
opposite sign. 
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FIG. 25 (Color online) Off-diagonal resistivity in mag¬ 
netic field, measured in a graphene-based double-layer device. 
Both graphene sheets are kept at the same carrier density 
ni = 712 = n and at T = 240K. Solid symbols represent the 
experimental data. [Reproduced from Titov et al. (2013a).] 


value in zero field, see e.g. Fig. 22. 

The classical, two-band mechanism of magnetodrag in 
graphene at charge neutrality can be readily illustrated in 
the case, where the system size is much larger than any 
characteristic length scale, such that the two graphene 
sheets may be considered effectively infinite. In this case 
(see Fig. 24), the driving current in the active layer, 
corresponds to the counter-propagating flow of electrons 
and holes with zero total momentum (due to the exact 
electron-hole symmetry). Once the weak magnetic field is 
applied, electrons and holes are deflected by the Lorentz 
force and drift in the same direction. The resulting quasi¬ 
particle flow. Pi, carries a nonzero net momentum in the 
direction perpendicular to Ji. This momentum can be 
transferred to the passive layer by the interlayer Coulomb 
interaction inducing the quasi-particle current, P 2 , in the 
same direction as Pi. The Lorentz forces acting on both 
types of carriers in the passive layer drive the charge flow 
in the direction opposite to Ji. If the passive circuit is 
open, this current is compensated by a finite drag voltage, 
yielding a positive drag resistivity (Titov et al., 2013a). 

This mechanism of magnetodrag at charge neutrality is 
closely related to the anomalous Nernst effect in single¬ 
layer graphene (Muller and Sachdev, 2008; Wei et al., 
2009; Zuev et al., 2009). Indeed, the quasi-particle cur¬ 
rent is proportional to the heat current at the Dirac point. 
The fact that the Lorentz force in the electron and hole 
bands has the opposite sign is also the reason for the 
vanishing Hall effect at charge neutrality. 

Despite being qualitatively clear, the above description 
of magnetodrag yields the induced drag voltage which has 
the sign opposite to that observed in experiment (Gor¬ 
bachev et al, 2012; Titov et al, 2013a), see Fig. 24. In 
fact, the negative drag in Fig. 24 can only appear if the 
quasiparticles currents in the two layers Pi and P 2 have 
opposite directions. According to Narozhny et al. (2015) 
and Titov et al. (2013a), this is what happens in small, 


mesoscopic samples used in experiment. 

Consider the continuity equation for the quasiparticle 
current Pi, including relaxation by electron-hole recom¬ 
bination (Foster and Aleiner, 2009; Titov et al., 2013a) 

VPi = -{pi - po)/Tph - {pi - P2)/(2tq), (106) 

where pi are the quasiparticle densities in the two layers, 
po = ttT^/{ 3Vg) is the equilibrium quasiparticle density 
at the Dirac point, Tph describes the energy loss from the 
system dominated by phonon scattering, and tq charac¬ 
terizes quasiparticle imbalance relaxation due to inter¬ 
layer Coulomb interaction. The equation for the pas¬ 
sive layer can be obtained by interchanging layer indices. 
In the absence of quasiparticle recombination, hard-wall 
boundary conditions at the sample boundaries allow only 
for the trivial solution. In contrast, taking into account 
inelastic processes, one finds the nontrivial solution illus¬ 
trated in Fig. 24: Pi = —P 2 . 

Combining the continuity equation (106) with the hy¬ 
drodynamic description of linear response transport in 
graphene (with the additional gradient terms that ac¬ 
count for inhomogeneity of physical quantities in finite- 
size systems), one can describe the negative drag ob¬ 
served in experiment (Narozhny et al., 2015; Titov et al., 
2013a), see Fig. 25. The exponential collapse of theo¬ 
retical curves at high carrier density is an artifact of the 
two-mode approximation adopted in Titov et al. (2013a). 
The more accurate three-mode approximation (Narozhny 
et al., 2015) includes thermoelectric effects formulated 
in terms of energy currents; the corresponding hydro- 
dynamic description yields only the power-law decay of 
the magnetodrag at /r, ^ T, in contrast to the exponen¬ 
tial collapse shown in Fig. 25. As compared to lower- 
temperature data (see Fig. 24), the results shown in 
Fig. 25 exhibit qualitatively new features which can be 
attributed to higher efficiency of relaxation processes at 
higher temperature. 

E. Hall drag in graphene 

Hall drag measurements in graphene were reported 
in Titov et al. (2013a). These experiments were per¬ 
formed at relatively high temperatures T = 160 — 240K, 
where macroscopic coherence is not expected to exist. 
While disorder effects in graphene are often attributed to 
Coulomb scatterers characterized by mean free time that 
is linear in energy, the measured Hall drag resistivity is 
not small as would follow from a mechanism similar to 
that suggested in Hu (1997) and von Oppen et al. (2001). 

Instead, double-layer graphene samples demonstrate a 
much simpler, yet still strong effect based on the coex¬ 
istence of electron and hole liquids in each layer (Foster 
and Aleiner, 2009). Consequently, the observed Hall drag 
resistance, Fig. 26, is large when one of the layers is close 
to the neutrality point and vanishes if two layers have 
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FIG. 26 (Color online) Left panel: Hall drag resistivity in 
graphene as a function of gate voltages controlling carrier den¬ 
sities in the two layers. White diagonal area corresponds to 
vanishing Hall drag for ni = —n 2 . Lines track positions of 
maxima in single-layer resistivity in top (open symbols) and 
bottom (solid symbols) layers. Right panel: Hall drag resis¬ 
tivity as a function of carrier density for ni = n 2 = n. Blue 
squares represent the experimental data. The red curve rep¬ 
resents the theoretical prediction. [Reproduced from Titov 
et al. (2013a).] 
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FIG. 27 (Color online) Schematic view of the drag resistivity 
at low temperatures. The dashed line illustrates the third- 
order drag effect. Left panel: The black solid line represents 
the lowest order contribution to drag. The arrows indicate the 
tendency of the two terms with the decrease of temperature 
T —>■ 0. Right panel: The purple solid line represents the 
contribution of correlated disorder. [Reproduced from Schiitt 
et al. (2013).] 


possible values of a^j. 


the same charge densities with opposite signs (a white 
line running from the top left to bottom right corner in 
the left panel Fig. 26). 

Hall drag effect in graphene can be understood with the 
help of the hydrodynamic theory (Narozhny et al., 2015; 
Titov et al, 2013a). Indeed, given the presence of two 
noncollinear currents in the model, it is not surprising 
to see the nonzero Hall drag away from the Dirac point, 
where both the conventional Hall effect and Hall drag 
change sign together with the carrier density and thus 
have to vanish Hall drag also has to vanish in the 
degenerate regime where only one band contributes to 
transport and the standard single-band theory (2) holds. 
However, this regime lies outside of the parameter range 
of the experiment (Titov et al, 2013a). Thus, some Hall 
drag signal is observed at all densities, but p^y decays 
to rather small values as the density increases beyond 
n ~ 1 X 10^^cm“^. Interestingly enough, the data show 
a sign change of p^y at n « ±2 x lO^^cm”^. This rather 
weak effect requires a more accurate consideration. 

The right panel of Fig. 26 shows the results of the hy¬ 
drodynamic theory alongside experimental data. This 
calculation was performed without any fitting (Titov 
et al., 2013a). The value of impurity scattering time r 
was determined from the measured single-layer resistiv¬ 
ity. The effective interaction parameter was estimated by 
the most plausible value for graphene on hBN, Ug « 0.2 
[see e.g. Kozikov et al. (2010) and Reed et al. (2010) for 
general considerations and the experimental evidence for 


A similar two-fluid model was used in Song et al. (2013) to ex¬ 
plain Hall drag in terms of the “energy-driven drag mechanism”. 
Indeed, if one omits the interlayer frictional force, one would still 
And nonzero Hall drag due to the interlayer energy relaxation. 


F. Higher-order effects in graphene 

All theories of Coulomb drag in graphene discussed so 
far were concerned with the leading-order contribution of 
the interlayer interaction. Indeed, even the nonperturba- 
tive results of the hydrodynamic approach were obtained 
by solving the kinetic equation with the collision inte¬ 
gral (95), where the transition probability was estimated 
using the Fermi Golden Rule. All such theories predict 
vanishing drag at the point of exact electron-hole sym¬ 
metry [with the exception of the academic case of pure 
graphene, see Eq. (104) and Fig. 23]. 

However, measurements (Gorbachev et al, 2012; Titov 
et al, 2013a) reveal nonzero drag at the double Dirac 
point, see Fig. 24. At the time of writing, there is no 
consensus in the community regarding the origin of this 
effect. At the same time, higher-order processes (see 
Sec. H.D) are known to be insensitive to the electron-hole 
symmetry and thus may provide a plausible explanation 
(Titov et al, 2013a). 


1. Third-order drag in graphene 


The third-order drag effect in graphene was considered 
in Schiitt et al (2013). The principle results are shown 
in Fig. 18 in red. A schematic illustration of the relative 
strength of the second- and third-order contributions is 
given in the left panel in Fig. 27. 

The third-order drag resistivity in the diffusive regime 
can be found similarly to the conventional case, see 
Sec. H.D. All microscopic details are masked by the 
diffusive nature of electronic motion. However, due to 
the relatively weak screening and the possibility to tune 
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the carrier density to the Dirac point, one finds a richer 
physical picture with multiple parameter regimes. 

The standard “Fermi-liquid” regime (Levchenko and 
Kamenev, 2008b) corresponds to the condition 

Nk ^ max , 

where = 4 describes spin and valley degeneracy of 
quasiparticle states in graphene. Here the temperature- 
independent result (40) is reproduced, although with the 
extra factors of N 

pf ~ {h/e^)N-^g-^Kd)-^. (107) 

At higher temperatures, one can achieve a different, 
high-temperature regime with 


30mK (Tan et al., 2007) [possibly due to the specific char¬ 
acter of impurities in graphene (Ostrovsky et al, 2007)]. 

For weak disorder or higher temperature the diffusive 
approximation fails. Drag in vicinity of the Dirac point 
can then be described by the quantum kinetic equa¬ 
tion approach. The previously reviewed results, e.g., 
Eq. (104) were obtained by approximating the collision 
integrals with the help of the Fermi Golden Rule, see 
Eq. (95b). However, taking into account next-order ma¬ 
trix elements yields a nonzero contribution, similar to the 

(3) 

above third-order result p)j . 

Taking into account the second-order matrix element, 
one can generalize the Golden Rule expression (96) by 
using the combination 

[C/i^ + C/W|2 |[/(i)|2 + 2Re{u[l\u[fr}. (112) 


rf-i <€.Nk<€. V¥/D. 


In this case, the resulting drag resistivity decays rapidly 


d) h I 1 

^ e2 g3 {Nxdy \ T ) 


(108) 


Since oc a and (x all relaxation rates will 
now get an additional contribution of the order of a^. In 
particular, the “drag rate” gets a contribution that 
is independent of the carrier density 

~ a^N{pLlTf + a^NT, (113) 


The experiment of Gorbachev et al. (2012) was per¬ 
formed on samples with the small interlayer spacing. In 
the limit xd 1 , one finds three different temperature 
regimes. 

Close to the Dirac point and at lowest temperatures, 
the drag resistivity is temperature-independent: 

p%\pL«.T-,TT<.a^)^hle^. (109) 

Note, that this result is also independent of the strength 
of the Coulomb interaction al 

At somewhat higher temperatures (or, equivalently, for 
slightly weaker disorder strength), the third-order contri¬ 
bution decays as function of temperature 

pgV < T < T-i < a-^T) - (Ve^)(a^Tr)-3/2. 

( 110 ) 

These results are illustrated in the right panel in Fig. 27. 

Away from the Dirac point, the third-order contribu¬ 
tion decays as a function of the chemical potential (or 
equivalently, carrier density) and quickly becomes sub¬ 
leading, see the left panel in Fig. 27: 


^/iT » max l,a ^(Tt)^/^ ^ 


h 1 

e2 (pt)3 ■ 


( 111 ) 


As a result, p\) may only be detectable at low tempera¬ 
tures and in vicinity of the Dirac point. 

While estimating = 0), the single-layer conduc¬ 

tivity was assumed to be of the order of the quantum 
conductance CT ~ e^//i, i.e. discarding localization effects. 
Indeed, single-layer measurements on high-quality sam¬ 
ples show temperature-independent conductivity down to 


which dominates near the Dirac point. In this case, one 
may neglect the conventional, second-order drag contri¬ 
bution; the result is (Schiitt et al., 2013) 


h a^T + 

^ T + a^p^rN ’ 


p < Tt > 1. 


Exactly at the Dirac point this yields 


PD^(h/e^)a^. (114) 


This result is illustrated in the right panel of Fig 27 by 
the horizontal asymptote at Tt ^ 1. 


2. Interlayer disorder correlations 


Within the conventional theory, charge carriers in each 
layer scatter off an independent disorder potential. This 
picture is clearly applicable to the cases where impu¬ 
rities are mostly concentrated in the substrate insula¬ 
tor sandwiching the double-layer structure. In the case 
of the standard double-well heterostructures (Eisenstein, 
1992; Gramila et al., 1991; Hill et al, 1997; Lilly et al., 
1998; Solomon et al, 1989), the random potential origi¬ 
nates in the delta-doped layers providing charge carriers. 
These layers are typically located on the outer sides of the 
double-well structure. In graphene, disorder potential is 
often attributed to the insulating substrate, in particular 
to Si 02 . Indeed, in graphene-based samples of Kim et al. 
(2011) and Kim and Tutuc (2012), graphene monolayers 
are exfoliated onto a thick Si 02 dielectric, while the in¬ 
terlayer spacer consists of 14mm-thick AI 2 O 3 . In such a 
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structure, the impurity potential created by the silicon 
oxide is likely to affect only the nearest monolayer. 

In contrast, the samples of Gorbachev et al. (2012) 
consist of graphene-hexagonal-boron-nitride heterostruc¬ 
tures, where the interlayer spacer contains only few 
atomic layers of the same insulator (boron nitride) that 
is used as a substrate. In this case, impurity potential 
originating from the interlayer spacer would be equally 
felt by carriers in both graphene layers. Another sce¬ 
nario for disorder correlation (Gorbachev et al, 2012; 
Song and Levitov, 2012) involves interactions between 
charge-density inhomogeneities forming due to impurity 
potential in the two layers. 

The effect of the correlated disorder in the drag mea¬ 
surements is insensitive to the electron-hole symmetry 
(Gornyi et al, 1999; Hu, 2000a), and thus may also pro¬ 
vide an explanation (Schiitt et al, 2013; Song and Levi¬ 
tov, 2012) for the observed nonzero drag in graphene at 
the Dirac point (Gorbachev et al, 2012). 

At high temperatures, Tt ^ 1, the effect of the corre¬ 
lated disorder can be described by the skeleton diagram 
similar to the third-order drag contribution, see the right 
panel of Fig. 7. Interlayer disorder correlations can be 
incorporated into the scattering amplitude, but now in¬ 
stead of the second-order matrix element in Eq. (112), 
one has to introduce an interlayer disorder scattering rate 
l/(rri 2 ). The resulting “drag rate” is given by 

l/rl,”'- ~ a^T/{TTi2) = aVn2, 
corresponding to the drag resistivity 
~ aV(Tri2), 

which overcomes the third-order drag contribution pjj' ~ 
at l/ri 2 > aT. This happens in the perturbative 
regime (1/r > a^T for moderately correlated disorder, 
ti 2 ~ r), where the correlated-disorder contribution can 
be calculated diagrammatically. 

Macroscopic inhomogeneities can be described in terms 
of macroscopic spatial fluctuations Spi in chemical poten¬ 
tials of the two layers (Song and Levitov, 2012), charac¬ 
terized by the correlation function 

F‘>f{r - r') = {dfi,{r)6pj{r')) ^ 0. (115) 

Assuming the spatial scale of the fluctuations to be much 
larger than all characteristic scales related to the parti¬ 
cle scattering, one can solve the hydrodynamic equations 
locally, yielding the local drag rate 

l/TD{r) ^ a'^NfXi{r)p2{r)/T. (116) 

Averaging over the small fluctuations of the correlated 
chemical potentials, one arrives (Schiitt et al, 2013) at 
the correction to the universal third-order result (114), 

ATO(f. = 0 ) ~ 


Finally, in the ultraclean limit 

l/r < a^NFl^^/T, (118) 


one can approximate the local drag resistivity by an ana¬ 
log of Eq. (104): 


ApD{r;p. = 0 ) 


h 2 5ai6u2 

— (X -. 

SpiSpi + 6p,2Sp2 


(119) 


In particular, for perfectly correlated chemical poten¬ 
tials, 6p,i{r) = 5/^2(r), the fluctuations drops out from 
Eq. (119) and the local resistivity turns out to be inde¬ 
pendent of r. In a more general case, the averaging over 
fluctuations becomes nontrivial, but this can only affect 
the numerical prefactor in the final result. Thus, the cor¬ 
related large-scale fluctuations of the chemical potentials 
in the layers in effect shift the curve 1 in Eig. 18 upwards, 
extending the validity of the fully equilibrated result. 


PD ~ {h/e^)a^, (120) 


to the case of finite disorder, Eq. (118), at the Dirac 
point. This implies that in the case of correlated inhomo¬ 
geneities the disorder-induced dip in the lower left panel 
of Fig. 23 develops only for sufficiently strong disorder. 


V. COULOMB DRAG AT THE NANOSCALE 

The effects of Coulomb interaction are especially pro¬ 
nounced at the nanoscale. In quantum dot devices one 
can utilize the Coulomb-modified Fano resonance to de¬ 
tect the electric charge (Johnson et al, 2004). Two- 
level pulse technique was used to detect individual elec¬ 
tron spin (Elzerman et al, 2004). Quantum dots were 
also used as high-frequency noise detectors (Onac et al, 
2006). Transport measurements on adjacent but elec¬ 
trically isolated quantum point contacts (QPCs) exhibit 
a counterflow of electrons [i.e. detector current flowing 
in the direction opposite to the driving current (Khra- 
pai et al, 2007)]. In nanosize CdSe-CdS semiconductor 
tetapods (Mauser et al, 2010), Coulomb drag-like effects 
lead to photoluminescent emission. 

Theoretically, Coulomb drag in a system of two elec¬ 
trically isolated QPCs was considered in (Levchenko and 
Kamenev, 2008a). Within the linear response the drag 
mechanism was found to be similar to that in the bulk 
2D electron systems. Remarkably, already for seemingly 
modest drive voltages (much smaller than temperature) 
the system crosses over to the nonlinear regime, where 
the effect is dominated by the excess shot noise of the 
drive circuit. Nonlinear transport was also found to be 
crucial for drag effects in a system of parallel quantum 
dots (Moldoveanu and Tanatar, 2009). An exciting new 
development is the proposal to use the drag effects to 
study transport properties of polaritons in optical cavi¬ 
ties and, in particular, their superfluidity (Berman et al, 
2010 a,b). 
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A. Quantum dots and quantum point contacts 

Interactively coupled mesoscopic and nanoscale cir¬ 
cuits, such as quantum wires (Debray et ai, 2000, 2001; 
Laroche et al., 2011; Morimoto et al., 2003; Yamamoto 
et al., 2006), quantum dots (Aguado and Kouwenhoven, 
2000; Onac et al, 2006) or point contacts (Khrapai 
et al., 2006, 2007), provided new fruitful ways of study¬ 
ing Coulomb drag phenomena and revealed a plethora of 
interesting physics. These devices typically have dimen¬ 
sions smaller than the temperature length Lt = vf/T 
and voltage-related length scale Ly = vp/ieV), and dif¬ 
fer substantially form their two-dimensional quantum- 
well counterparts in several important ways, (i) The 
strength of Coulomb interaction is naturally enhanced by 
reducing system size that should lead to more profound 
dragging effect, (ii) Transmission across the device in 
the drag (drive) circuit or both can be efhciently con¬ 
trolled by the gate voltages that allows to open quantum 
conduction channels one by one. (iii) The electron-hole 
symmetry in such devices is broken much stronger than in 
bulk systems. In mesoscopic devices this is due to a ran¬ 
dom configurations of impurities, while in the quantum 
nanocircuits the effect is due to the energy dependence of 
transmission coefficients in Landauer picture of transport 
(Biittiker et al., 1985; Landauer, 1957, 1970). (iv) Be¬ 
cause of the above reasons, the quantum circuits may be 
easily driven out of the linear response domain and cor¬ 
responding voltage scale is parametrically smaller than 
the temperature, (v) Ultimately, the mere mechanism of 
drag in the nonlinear regime is different and governed by 
the quantum noise fluctuations. 

The most peculiar feature of the observed Coulomb 
drag in such systems was that the drag current exhibited 
maxima for specific values of the gate voltage, where the 
drive circuit was tuned to an opening of another con¬ 
ductance channel, see Fig. 28 for the illustration. This 
hinted the importance of the electron shot noise in the 
drive circuit, which was known to exhibit a qualitatively 
similar behavior (Lesovik, 1989; Reznikov et al., 1995). 
Indeed, electron current shot noise power is proportional 
to the product of the transmission and reflection coeffi¬ 
cients that is peaked between the conductance plateaus. 
It was argued early on that drag may be interpreted as a 
rectihcation of nearly equilibrium classical thermal fluc¬ 
tuations in the drive circuit (Kamenev and Oreg, 1995). 
The extension of this idea to rectification of the quantum 
shot noise was plausible and happened to be correct in 
a certain regime. The subtlety of this picture was that 
such a rectification is only possible due to electron-hole 
asymmetry in both circuits, otherwise drag currents of 
electrons and holes cancel each other. The mismatch 
between transmission probabilities of electron and hole 
excitations is maximal at the verge of an opening of 
the new conduction channel, which implies that spikes 
of drag conductance may originate form the asymmetry 




Plunger gate voltage (mV) 


FIG. 28 Left panel inset-(a) represents scanning electron mi¬ 
crograph of the gate structure defined on top of the semi¬ 
conductor heterostructure. The gates highlighted by dashed 
lines are used to define a quantum dot (QD) on the left and a 
quantum point contact (QPC) on the right. Inset-(b) shows 
current Iqd versus plunger gate voltage whereas inset-(c) dis¬ 
plays QPC conductance Gqpc as a function of the gate volt¬ 
age. In such device the QPC is used as a noise generator and 
the QD as a detector. Right panel shows current through the 
QD, as a function of the plunger gate voltage, under the influ¬ 
ence of shot noise generated by the QPC with characteristic 
peaks. [Reproduced from Onac et al. (2006).] 


alone rather than shot noise. 

In order to get an insight into these delicate details con¬ 
sider the linear response regime when drag conductance 
go can be expressed as follows (Levchenko and Kamenev, 
2008a) 


go 


dw |Zi2(a;)P ri(a;)r2(u;) 
SttT sinh^(w/2T)' 


( 121 ) 


Here ^ 12 ( 0 ;) is the interactively-induced trans-impedance 
relating local fluctuating currents and voltages between 
the circuits (Geigenmiiller and Nazarov, 1991). The cor¬ 
responding rectihcation coefficients are given explicitly 

by 

Q n 

( 122 ) 

where Rq = 27r/e^ is the quantum of resistance, e± = 
e±a;/2, /(e) is the Fermi distribution function, and |ti„p 
is the energy dependent transmission coefficient in the 
transversal channel n of the circuit i = 1,2. This expres¬ 
sion admits a transparent interpretation: potential huc- 
tuations with frequency uj, say on the left of the quantum 
point contact, create electron-hole pairs with energies e± 
on the branch of right moving particles. Consequently 
the electrons can pass through the quantum point con¬ 
tact with the probability |ti„(e+)p , while the holes with 
the probability |ti„(e_)p. The difference between the 
two gives the net current flowing across the contact while 
the Fermi functions in Eq. (122) take care of the statisti¬ 
cal occupation of participating scattering states. Notice 
that unlike in the Landauer formula for conductance of a 
single quantum point contact where transmissions can be 
treated as being energy independent, the energy depen- 
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dence of these probabilities in the drag formula is crucial 
in order to have the asymmetry between electrons and 
holes, and thus nonzero rectification ri(a;). A particu¬ 
lar functional dependence of T on frequency depends on 
a model and details of device circuitry. It is instructive 
to focus on a limit of a single partially open channel in 
a smooth adiabatic quantum point contact. One may 
think then of the potential scattering barrier across it as 
being practically parabolic. In such a case its transmis¬ 
sion probability is given by 

|tm(e)P = [exp[(eVgi - e)/Ai] -k 1]“^ (123) 

where is an energy scale associated with the curva¬ 
ture of the parabolic barrier in the point contact i, while 
gate voltages Vgi move the top of the barrier relative 
to the Fermi energy within each of the point contact. 
This form of transmission was used to explain quantum 
point contact conductance quantization (Glazman et al., 
1988) and it turns out to be useful in application to the 
Coulomb drag problem. For the low temperature limit 
T ^ Ai using Eq. (123) in Eq. (122) and carrying out 
energy integration yields 

cos\v{eVgi//S.i) + cosh(a;/Ai) 
cosh(eVgi/Ai) -I- 1 

(124) 

In the opposite limit when T ^ Ai one should replace 
Ai —)■ T. One should notice that for small frequency 
Fi oc whereas trans-impedance ^ 12 ( 0 ;) is practically 
independent of frequency in this limit since its charac¬ 
teristic scale is typically set by the inverse i?C'-time of 
circuits. Assuming that T <C max{Ai,r^p} one arrives 
at 

go 1 

5 Q 6 A1A2 cosh^(eVgi/Ai) cosh^(el42/A2) ’ 

(125) 

where u = Zi 2 {Q)/Rq. The resulting expression for the 
drag conductance exhibits peak as a function of gate 
voltage in drag or drive quantum point contact. Yet at 
this level it has nothing to do with the shot noise peaks, 
but rather reflects rectification of near-equilibrium ther¬ 
mal fluctuations (hence proportionality to T^) along with 
the electron-hole asymmetry (hence a nonmonotonic de¬ 
pendence on Vgi). However, one should realize that the 
crossover to the nonlinear regime of transport in such 
devises can occur at rather low voltages eV* ~ T^/A^ 
such that Eq. (121) becomes inapplicable already at 
V > V*. More general considerations by (Chudnovskiy, 
2009; Levchenko and Kamenev, 2008a; Sanchez et ai, 
2010 ) revealed that for the out of equilibrium nonlinear 
regime the drag current is due to the rectification of the 
quantum shot noise and hence proportional to the Fano 
factor ~ |tniP]- K again exhibits a generic 

nonmonotonic behavior of drag with multiple peaks but 
for the entirely different reason independent of asymme¬ 


try factor. Nonlinear transport was also found to be cru¬ 
cial for drag effects in a system of parallel quantum dots 
(Moldoveanu and Tanatar, 2009). 

Drag phenomena in quantum circuits can be naturally 
connected to our earlier discussion of drag in mesoscopic 
systems in Sec. III. Indeed, one or both circuits may 
be represented by a multichannel quasi-one-dimensional 
(or two-dimensional) mesoscopic sample. In this case 
J2n |tn(e)P = 5 (e) is a dimensionless (in units Rq^) con¬ 
ductance of the sample as a function of its Fermi energy. 
As discussed above, such conductance exhibits universal 
fluctuations, that is g{e) = g + Sg{e), where ^ 1 is 
an average conductance and Sg{e) ~ 1 is a sample and 
energy-dependent fluctuating part. Since the character¬ 
istic scale of the energy dependence of the fluctuating 
part is the Thouless energy Et one naturally finds from 
Eq. (122) that corresponding mesoscopic fluctuations of 
the rectification coefficient are of the order 

riuj)^±-^^ (126) 

This result ultimately leads to the an estimate of the 
variance of drag in the form of Eq. (58) 

Quantum Coulomb drag circuits provide rich platform 
to explore nanoscale transport far beyond ideas of using 
them for high-frequency noise sensing. In particular, a 
different drag effect may also be observed in the absence 
of any drive current if one nanocircuit is made hotter 
than the other - the cold circuit is expected to rectify 
the thermal charge fluctuations of the hot circuit (Soth- 
mann et al., 2012). Furthermore, interactively coupled 
devices provide unique tools to test nonlinear fluctuation- 
dissipation relations and its closely related Onsager sym¬ 
metry relations in the far from equilibrium conditions 
when detailed balance is explicitly broken (Bulnes Cue- 
tara et ai, 2013; Sanchez et al., 2010). 

Other intriguing examples include nanosize CdSe-CdS 
semiconductor tetapods (Mauser et al., 2010) where 
Coulomb drag-like effects lead to photoluminescent emis¬ 
sion. As an alternative to optical probes, electrical 
read-out of a single electron spin becomes possible in a 
Coulomb drag-like devices of interactively coupled QPC 
and QD (Elzerman et al., 2004). 

B. Optical cavities 

Coulomb interaction is not exclusive to electrons and 
can couple any charges. Moreover, even neutral, com¬ 
posite objects may interact with charges by means of 
an effective “polarization” or “charge-dipole” interac¬ 
tion (Margenau and Kestner, 1969), which ultimately 
stems from the Coulomb interaction between an external 
charge and individual charged constituents of the com¬ 
posite object. In particular, long-ranged interactions be¬ 
tween spatially separated electrons and polaritons may 


i' ^ 2eAi 

Taw) = In 
no 
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lead to interesting drag effects (Berman et ai, 2010a,b; 
Kulakovskii and Lozovik, 2004) that can be used, e.g., for 
designing electrically controlled optical switches (Berman 
et ai, 2014). 

Two-dimensional excitonic polaritons have been a sub¬ 
ject of intensive research (Amo et al, 2009; Balili et ai, 
2007; Kasprzak et al, 2006; Snoke, 2008). These ex¬ 
citations appear as a result of resonant exciton-photon 
interaction in a system consisting of an optical micro¬ 
cavity and a quantum well embedded within the cavity. 
The lower polariton branch is characterized by extremely 
small effective mass raising the possibility of achieving 
the Bose-Einstein condensation and superfluidity at rel¬ 
atively high temperatures (Balili et ai, 2007; Littlewood, 
2007). 

The optically excited excitons in microcavities should 
not be confused with the spontaneously formed excitons 
in double quantum wells discussed in Sec. VILA. In par¬ 
ticular, these excitons can be excited by laser pumping 
in the single quantum well embedded within the cavity. 
Now, if a second quantum well is added to the device 
(Berman et al, 2010a), then Coulomb interaction binding 
electrons and holes into excitons can be screened (Finkel- 
stein et al, 1995; Gubarev et al, 2000) by a 2DEG popu¬ 
lating the second well. As a result, the excitonic binding 
energy is reduced and as the density of the 2DEG ap¬ 
proaches a certain critical value, the excitons may disap¬ 
pear altogether. The excitonic collapse manifests itself 
through disappearance of the corresponding line in the 
photoluminescence spectrum. 

Keeping the electron density below the above critical 
value, one obtains a system containing coexisting, spa¬ 
tially separated excitons and electrons. The effective in¬ 
teraction between electrons and excitons was considered 
in Lozovik and Nikitkov (1999). This interaction leads 
to mutual friction between the two systems that can be 
observed by selectively exciting one of them by external 
probes. 

By focusing laser pumping on a particular region 
within the cavity, one can generate a gradient of exciton 
and polariton densities. These gradients induce a flow 
of both polaritons and excitons. The long-range interac¬ 
tion between the excitons (or the exciton component of 
the polaritons) may transfer energy and momentum to 
the electronic system in the second quantum well, gener¬ 
ating an electric current or inducing voltage, similarly to 
the standard drag effect discussed in Sec. 11. 

Alternatively, one can drive a current through the 
2DEG. In this case, the mutual friction will lead to the 
appearance of the exciton flow. These excitons are entan¬ 
gled with cavity photons and their flow will create a flow 
of polaritons. In other words, the long-ranged electron- 
exciton interaction allows one to effectively “move” the 
cavity photons by applying electric current to the 2DEG 
(Berman et al, 2010a,b). The drag effects in microcavi¬ 
ties are schematically illustrated in Fig. 29. 


Bragg mirror Bragg mirror 
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FIG. 29 Left: quasiparticle flow in the cavity polariton sub¬ 
system induced by the electric current in the 2DEG at low 
temperatures. Right: electric current in the 2DEG induced 
by the optically excited flow in the polariton subsystem. [Re¬ 
produced from Berman et al (2010b).] 
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FIG. 30 (Color online) Schematic of the wedge-shaped micro¬ 
cavity formed by two distributed Bragg reflectors (DBR) that 
encompasses the embedded quantum wells The excitons are 
located in the quantum well (gray) between the reflectors. A 
metal layer deposited on the upper DBR creates a Y-shaped 
potential energy landscape for the polaritons. The driving 
current runs perpendicularly to the stem of the channel in 
the quantum well (green). Reprinted with permission from 
Berman et al. (2014). Copyright (2014) American Chemical 
Society.] 


Recently, Berman et al (2014) have proposed to use 
the drag effect in optical cavities for building an elec¬ 
trically controlled optical switch, see Fig. 30. The po¬ 
laritons are assumed to be created at a constant rate by 
external laser pumping. The wedge-like shape of the mi¬ 
crocavity is chosen in order the create a force driving 
the polaritons along the cavity towards the Y-junction. 
Without the drag effect, the polariton flux is distributed 
equally between the two branches of the junction. Driv¬ 
ing an electric current through a second quantum well 
(shown in green in Fig. 30) results in a drag force in the 
junction region that effectively redistributes the polari¬ 
tons flux between the branches. Berman et al (2014) 
find that for realistic parameters of the device one can 
achieve 90% accuracy of the switching of the polariton 
flow. 
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VI. COULOMB DRAG BETWEEN PARALLEL 
NANOWIRES 

It is well-known that physics of electrons confined to 
one spatial dimension (ID) is dominated by interactions. 
Coulomb drag between two closely spaced but electri¬ 
cally isolated quantum wires was used to observe Wigner 
crystallization (Yamamoto et ai, 2002, 2006, 2012) and 
Luttinger-liquid effects (Debray et ai, 2001; Laroche 
et ai, 2008, 2014). The effect was also used to study ID 
sub-bands in quasi-lD wires (Debray et al, 2000; Laroche 
et al, 2011). 

Early theoretical work on drag between ID systems 
(Gurevich and Muradov, 2000, 2005; Gurevich et ai, 
1998; Hu and Flensberg, 1996; Raichev and Vasilopou- 
los, 2000, 1999) was based on the Fermi-liquid approach 
and targeted multiple-channel wires at high enough tem¬ 
peratures, where electron correlation effects (other than 
screening) are not important. Tanatar (1998) consid¬ 
ered the role of disorder. It is however well-known, that 
the Fermi-liquid theory fails for purely ID systems, i.e. 
single-channel wires (Giamarchi, 2004), quasi-lD wires 
with single ID subband occupancy (Laroche et al., 2014), 
and systems comprising a small number of coupled ID 
channels. Coulomb drag between two Luttinger liquids 
with point-like interaction region was discussed in Flens¬ 
berg (1998) and Komnik and Egger (2001). Nazarov and 
Averin (1998) considered two independent Luttinger liq¬ 
uids coupled by interwire backscattering. Schlottmann 
(2004a,b) used Bethe-Ansatz methods to solve the prob¬ 
lem of two wires coupled by a particular d-function po¬ 
tential. Especially interesting is the prediction of the 
Mott-insulator-type state corresponding to formation of 
two interlocked charge density waves (CDW) in quantum 
wires (Fuchs et al, 2005; Klesse and Stern, 2000). 

A theory of Coulomb drag based on the Tomonaga- 
Luttinger liquid (TLL) theory (Haldane, 1981a,b; Lut¬ 
tinger, 1963; Tomonaga, 1950) predicts a behavior that 
qualitatively deviates from that in higher dimensions. 
Below a certain crossover scale T*, the drag resistivity 
between infinitely long quantum wires of equal electron 
density is predicted to increase exponentially with de¬ 
creasing temperature (Klesse and Stern, 2000) 

Pd ~ pTexp(A/r). (127) 

The energy gap A and crossover temperature T* are 
complicated functions of the interwire distance d, width 
of wires w, effective Bohr radius of the host mate¬ 
rial, and electron density n. For widely separated wires 
{kpd 1) they are exponentially suppressed 

A^T* ^ Epexpl-kFd/il-K)], (128) 

and the drag resistivity exhibits the high-temperature 


power-law behavior^® 

PD - {h/e^)kF\\T/EFf^-^, (129) 

for all practically relevant scales. Here K is the TLL 
interaction parameter in the relative charge sector deter¬ 
mined by the difference of the small-momentum intra- 
and interwire couplings and A is dimensionless interwire 
backscattering potential strength. 

The physical picture behind Eq. (127) is that at low 
temperatures x <T* the electrons in both wires form 
zigzag-ordered interlocked charge density wave (GDW). 
Then a relative charge displacement can be created 
only by overcoming a potential barrier, which ultimately 
translates into transport via activation, and consequently 
into Arrhenius-like behavior of drag. 

For short wires, Klesse and Stern (2000) report a qual¬ 
itatively different behavior. Here the CDW in one wire 
may slip as a whole relative to the CDW in the other wire. 
These instantaneous slips may stem from either thermal 
fluctuations or tunneling events. The latter leads to the 
drag resistance that tends to a finite, but exponentially 
large (in the wire length L) value as T —>■ 0. In contrast, 
Ponomarenko and Averin (2000) find a vanishing drag 
resistance, pu ^ T'^, regardless of whether the CDW is 
formed or not. 

For wires with different electron densities, Fuchs et al. 
(2005) find that the drag resistance (127) is suppressed 
by an additional exponential factor exp{—\Sp,\/T), where 
6p, = pi — p 2 is the difference between the chemical po¬ 
tentials Pi in the two wires. The high-temperature result 
(129) also becomes exponentially suppressed as soon as 
\dp\ exceeds the temperature. 

Allowing for a spin degree of freedom adds extra com¬ 
plexity to the problem, since the system might be un¬ 
stable towards gap opening in the spin sectors^^. If this 
does not happen (or at temperatures exceeding the spin 
gaps), the system shows the same qualitative behavior as 
above, but the exponent in Eq. (129) changes to 2K — 1. 
However, if the single wires develop spin gaps, the drag 
resistivity vanishes at T = 0 (Klesse and Stern, 2000). 

At temperatures above T* , the charge sector is gapless 
and the system can be described as two coupled wires 
in the TLL phase. For quasiparticles with linear dis¬ 
persion the only process contributing to drag is the in¬ 
terwire backscattering characterized by large momentum 
transfers q ~ 2kF. This process can be described by the 
usual drag formula (15), where one typically assumes the 
nonlinear susceptibility to be proportional to the imagi¬ 
nary part of the polarization operator (Fiete et al., 2006; 


In ID, p£) = — L)(dV 2 /dIi), is the drag resistivity 

per unit length (Klesse and Stern, 2000; Pustilnik et al., 2003). 
For a comprehensive discussion of ground state properties of ca- 
pacitively coupled ID systems see Carr et al. (2013) and Gia¬ 
marchi (2004). 



45 


Pustilnik et ai, 2003): 


h f dqduj (fVq [Imn(q,w)]^ 

J 47 r3 n?T sinh^(w/2r)’ 


(130) 


where Vq describes the interwire interaction. In the 
limit qd^ 1, the asymptotic form of Vq is given by 
Vq = le)yj2TTl(qd) exp{—qd). The polarization oper¬ 
ator for the TLL model is known (Giamarchi, 2004). For 
spinless fermions, the spectral weight of 2kF density fluc¬ 
tuations is given by 


Imn(g±,a;) 


sin ttK / 2'KaT 
Att'^u \ u 


(131) 
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where a ^ is the short-distance cut-off of the TLL 
theory, q± = q ±2kF, u is the renormalized Fermi ve¬ 
locity, and B{x,y) is the Euler bet a-function. Using 
Eq. (131) in the expression (130) one recovers Eq. (129), 
which was obtained by Klesse and Stern (2000) by means 
of a renormalization group analysis. In the perturba¬ 
tive approach, the interaction parameter A in Eq. (129) 
is given by A = V^kF^F- This leads to the exponential 
dependence oi pf) on distance separating the wires [since 
V2kF ^ exp{—2kFd)]. The regime of spin-incoherent Lut- 
tinger liquid and effect of disorder modify temperature 
dependence of Eq. (129) (Eiete et al, 2006). In the 
weakly interacting limit {K ~ 1) the drag resistivity is 
expected to grow linearly with temperature (Gurevich 
et al, 1998; Hu and Elensberg, 1996). 

In recent years, a lot of the attention was devoted to ID 
liquids with nonlinear dispersion [for reviews on this topic 
see Deshpande et ai (2010); Imambekov et ai (2012); 
and Matveev (2013) ]. In the TLL theory, the curvature 
of the quasiparticle spectrum is described by an irrelevant 
operator (in the renormalization group sense). However, 
at high enough temperatures it might lead to important 
effects and even mask the pure Luttinger behavior. In 
the context of Goulomb drag (Aristov, 2007; Dmitriev 
et al, 2012; Pereira and Sela, 2010; Pustilnik et al, 2003; 
Rozhkov, 2008, 2009), this is particularly important since 
nonlinearity of the band kinematically allows drag with 
small momentum transfer, q ^ T/vf ^ kF- 

Analytic calculation of the dynamical structure factor 
Imn((7, w) for arbitrary interactions and nonlinear dis¬ 
persion is a major challenge. However, such calculation 
is readily available in the case of weakly interacting elec¬ 
trons. At finite temperatures, but with the accuracy of 
the order T ^ mvp, the one-loop diagram yields 


Imn((7, uj) 


m _ sinh(a;/2r) _ 

4k cosh(qvF^+/2T) cosh(qvF^-/2T) ’ 


where = 2m{uj — VFq)/q^ ± 1- It is now tempting to 
follow the conventional path and use this result for ImH 
in the expression for the drag (130) to obtain 

PD {hkF/e‘^){Vo/vFf{T/EFf, (133) 


with the conclusion that curvature effects restore the 
Eermi-liquid behavior of drag in ID wires; further¬ 
more, the contribution (133) would dominate over the 
backscattering component (129) already at temperatures 
above T > [in Eq. (133) Vq should be un¬ 

derstood as Vq^T/vp]- even higher temperatures, 
{v/d) <T < Ef, the same approach yields a saturating 
drag resistivity, po ^ {hkF/e^){VQ/vFY{v/dEF^, fol¬ 
lowed by a falloff pu oc at T > Ef- For non¬ 

identical wires there appears an additional energy scale, 
Ts = kF\6v\ describing splitting between symmetric and 
antisymmetric plasmons modes in the double-wire sys¬ 
tem, which is determined by the difference between Fermi 
velocities in the wires Sv = vfi — vf 2 - Assuming that 
Ts <^Td, Eq. (133) holds only for Ts < T < T^, whereas 
below Ts drag resistivity due to forward scattering de¬ 
creases as pd oc T^ with lowering temperature. 

However, as shown by Dmitriev et ai (2012) the above 
conclusions about the forward scattering contribution to 
drag may be premature. The reason is subtle: the ex¬ 
pression (130) was derived under the tacit assumption 
that the intralayer relaxation processes due to electron- 
electron interaction are faster than the interwire momen¬ 
tum transfer. Now, in purely ID systems relaxation is 
determined by three-body collisions (Levchenko et al, 
2011; Lunde et al, 2007; Micklitz et al, 2010; Rieder 
et al, 2014) as inelastic two-body interaction is forbid¬ 
den by energy and momentum conservation. Same kine¬ 
matic restrictions require that intrawire backscattering 
responsible for equilibration involves states deep at the 
bottom of the band. Because of the Pauli statistics, the 
probability to find such a state unoccupied is exponen¬ 
tially small. Gonsequently, the equilibration rate in 
ID is exponentially suppressed, oc and the 

expressions (130) and (133) are difficult to justify. 

At the same time, interwire backscattering with small 
momentum transfer q ~ TIVf ^ If is also allowed in 
ID systems with nonlinear spectrum. This process in¬ 
volves a pair of scattering states: one near the Eermi level 
and another at the bottom of the band. Dmitriev et al 
(2012) found a solution of two coupled kinetic equations 
[cf. Eqs. (9)] yielding the drag resistivity in the form 


hkF f ^o\ 1 Td I Ef -2EfIt 


(134) 


By comparing the exponential factors in Eqs. (134) and 
(129), one can see that backscattering-induced drag fric¬ 
tion due to soft collision (namely collisions with small 
momentum transfer) dominate over direct backscattering 
with 2kF momentum transfer at temperatures T > Td- 
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FIG. 31 (Color online) Left: measurement schematic (top) 
and single-wire conductance quantization (bottom). Right: 
temperature dependence of the drag signal for three different 
samples (the magnitude of po in sample 2-L is divided by 200 
for visibility). For samples 2-L and 2-C, the temperature de¬ 
pendence was taken with no more than one ID subband occu¬ 
pancy in each wire, whereas the number of ID subband occu¬ 
pied in sample 3-R is known to be bounded by 0 < Ndrive < 2 
and 0 < Ndrag < 3. [From Laroche et al. (2014). Reprinted 
with permission from AAAS.] 


This is despite the fact that the contribution of soft col¬ 
lisions being strongly suppressed compared to Eq. (133). 
In addition, depending on the delicate interplay between 
the rates of two-particle interwire backscattering with 
small momentum transfer and triple-body intrawire chi¬ 
rality changing soft collisions, one may identify a param¬ 
eter regime where pD oc is suppressed with only 

one exponential factor, namely (Dmitriev et al., 2012). 

Coulomb drag between true ID systems was recently 
observed by Laroche et al. (2014) in a system of verti¬ 
cally integrated quantum wires where each wire has less 
than one ID subband occupied. The most striking the¬ 
oretical prediction, i.e., the upturn in the temperature 
dependence was revealed below the crossover tempera¬ 
ture T* ~ 1.6K, see Fig. 31. However, a quantitative 
comparison between the data and above theoretical re¬ 
sults proved to be difficult. Using the experimental esti¬ 
mates for the carrier density riio = y/n^D and interwire 
distance d ~ 40nm, one arrives at kpd ~ 2. Then, from 
Eq. (128) one finds the values for the Luttinger parameter 
K ~ 0.1 — 0.2 (for samples 3-R and 2-C) corresponding 
to very strong interaction that is beyond the applicabil¬ 
ity of the bosonization theory of Klesse and Stern (2000). 
On the other hand, fitting the high-temperature data to 
the power-law behavior (129) yields K ~ 1.5. This esti¬ 
mate, however, should be approached with caution, since 
Eq. (129) was derived for identical wires which is not the 
case in experiment, where electronic densities in the par¬ 
ent 2D layers differ by about 20%. In that case one ex¬ 
pects (Fuchs et al, 2005) an exponential suppression of 
Pd- All these issues remain to be clarified both theoreti¬ 
cally and experimentally. 


VII. NOVEL MANY-BODY STATES IN DOUBLE-LAYER 
SYSTEMS 

When a double-layer system is subjected to a strong 
magnetic field, the standard theoretical description of 
Coulomb drag (Jauho and Smith, 1993; Kamenev and 
Oreg, 1995; Zheng and MacDonald, 1993) fails: in con¬ 
trast to naive expectations, numerous experiments (Feng 
et al, 1998; Hill et al., 1996, 1998; Jorger et al, 2000c; 
Lok et al, 2001a,b; Patel et al, 1997; Pillarisetty et al, 
2003; Rubel et al, 1997a, 1998, 1997b) show significant 
dependence of the measured drag resistivity pD on the 
applied held, especially in the extreme quantum regime 
(Lilly et al, 1998; Murphy et al, 1994; Nandi et al, 
2012 ). 

Further experiments revealed the existence of novel 
quantum Hall states that are specihc to bilayer systems 
and have no analog in single-layer samples. Early work 
in this direction was reviewed in Eisenstein (1992, 1997). 
Remarkably, the bilayer many-body states exhibiting the 
quantum Hall effect (Murphy et al, 1994) may at the 
same time support a condensate of indirect (or interlayer) 
excitons (Finck et al, 2010; Nandi et al, 2012; Wiersma 
et al, 2007). An interlayer exciton is a bound pair of an 
electron from one layer and a hole from another layer of 
the device. The exciton carries no electric charge. Nev¬ 
ertheless, exciton transport (especially in the superhuid 
state) leads to interesting electrical effects. The exper¬ 
imental situation in the held is reviewed in Eisenstein 
(2014) and Eisenstein and MacDonald (2004). Here we 
focus on the manifestations of this exciting new physics 
in the drag measurements. 


A. Quantum Hall Effect in double-layer systems 

In a seminal paper, Halperin (1983) has suggested a 
generalization of the Laughlin wave function for the anal¬ 
ysis of multi-component systems. The simplest example 
of an extra degree of freedom that can be accounted for 
using this approach is the electron spin. A double-layer 
system provides another example, which is similar to the 
spin-1/2 in some respects and is signihcantly different in 
other. The two possible values of the layer index can 
be represented by the two orientations of a pseudo-spin 
(Stern et al, 2000; Yang et al, 1994). However, unlike 
the real spin, the double-layer system does not possess 
the SU{2) symmetry due to the difference between the 
intra- and interlayer matrix elements of the Coulomb in¬ 
teraction. Consequently, in the double-layer system the 
energy eigenstates do not have to be eigenstates of the 
total spin operator Sp (Girvin and MacDonald, 1997). 
As a result, states described by Halperin’s wave func¬ 
tions that are not eigenstates of St may be realized in 
double-layers (Eisenstein et al, 1992; Suen et al, 1992). 

In this review we are mostly interested in double-layer 
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systems where tunneling between the two layers is neg¬ 
ligible. Such systems support novel many-body quan¬ 
tum Hall states, that are specific to bilayers and arise 
due to interlayer Coulomb interaction (Chakraborty and 
Pietilainen, 1987; Haldane and Rezayi, 1987). Yoshioka 
et al. (1989) investigated a wide class of such states us¬ 
ing Halperin’s two-component wave functions (Halperin, 
1983). The ground state of the system crucially depends 
on the ratio of the interlayer separation and magnetic 
length d/^o- For a given filling factor, the magnetic 
length £o is proportional to the average separation be¬ 
tween electrons in one layer. Therefore, the ratio d/£o 
parametrizes the relative strength of intra- and interlayer 
Coulomb interaction. Assuming truly two-dimensional 
layers (i.e., setting aside complications that arise due to 
the finite width of the quantum wells in GaAs samples), 
one hnds that the interlayer many-body states are stable 
for d/1'0 ~ 1. At large d, the interlayer Coulomb interac¬ 
tion is inefficient and then the system behaves as if one 
connects two quantum Hall samples in parallel (Eisen- 
stein, 1997). This observation can be illustrated with 
the help of the typical phase diagram shown in Fig. 32 
for the case of the total hlling factor vt = £ (Murphy 
et al, 1994). In the opposite limit, d/lo —> 0, the sys¬ 
tem approaches the S'17(2)-symmetric point, and thus 
the Halperin states that are not eigenstates of St are 
expected to collapse (Eisenstein, 1997). 

Double-layers at the total hlling factor vt = I and with 
large interlayer separation (experimentally, d/^o ~ 2 — 4) 
behave as two weakly coupled systems of composite 
fermions (i.e. each layer is at = 1/2) while exhibit¬ 
ing strongly enhanced drag as compared to the zero-held 
case, see Sec. H.G. As the ratio d/fo is decreased, ex¬ 
periments (Kellogg et al, 2003, 2002b) show a gradual 
development of the Hall drag signal and a non-monotonic 
behavior of the longitudinal drag resistivity po- As d/£o 
approaches the transition into the strongly correlated, 
many-body state po shows strong enhancement, fol¬ 
lowed by a decrease. In the strongly coupled interlayer 
i/T = I state pd practically vanishes. At the same time, 
the Hall drag resistance develops a quantized plateau, see 
Fig. 33. Similar behavior was observed in Tutuc et al 
(2009, 2004); and Wiersma et al (2004) 

Early theoretical work on drag in quantum Hall states 
was focused on the non-dissipative drag (Duan, 1995; 
Renn, 1992; Yang, 1998; Yang and MacDonald, 2001). In 
contrast to the case of weak magnetic held (see Sec. H.G), 
a strong, quantized Hall drag has been identihed as a sig¬ 
nature of the interlayer correlated states. The 2x2 Hall 
resistivity matrix (for the two layers) has been shown 


Since the transition between the weakly and strongly coupled 
quantum Hall states is still poorly understood, one should be 
speaking in terms of the transition region instead of the precise 
critical value of d/£o. 
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FIG. 32 Phase diagram of the quantum Hall effect at = 1 
in double-layer systems. Asas is the tunnel splitting and 
e^/ (e£) is the Coulomb energy. Each symbol corresponds to a 
particular double-layer sample. Only the samples represented 
by solid symbols exhibit a quantized Hall plateau at ut ~ I- 
The interlayer quantum Hall state exists also in the absence 
of tunneling. [Reproduced from Murphy et al (1994).] 

(Renn, 1992; Yang, 1998) to be proportional to the Gram 
matrix (Gonway and Sloane, 1988; Read, 1990) describ¬ 
ing topological order in the quantum Hall state (Wen, 
1995): 

p// = {h/e'^)Kij ^ Pil = nhje?, n > 0. (135) 

Similar conclusion was reached in Yang and MacDonald 
(2001) on general topological grounds. 

Kim et al (2001) have suggested to use the drag resis¬ 
tivity to distinguish between various quantum Hall states 
in double-layer systems at p't = 1- For the compressible 
(weak-coupling) state at large interlayer separation, the 
Hall drag resistivity vanishes, while the longitudinal drag 
is determined by gauge-field fluctuations and is given by 
Eq. (53). The compressible state exhibits a strong pair¬ 
ing instability (Bonesteel, 1993; Greiter et al, 1991). If 
Landau-level mixing is substantial (as it often is in exper¬ 
imental samples), the paired state may be described by 
the (3,3,—1) Halperin wave function. This state resem¬ 
bles apx + ipy superconductor of composite fermions. As 
a result, it is expected to exhibit the quantized Hall drag 
resistivity (135) with n = —1. 

For smaller interlayer separation [d £q) the system 
undergoes a transition into an incompressible, corre¬ 
lated “quantum Hall ferromagnet” state described by the 
(1,1,1) Halperin wave function. This state possesses a 
gapless neutral mode and is characterized by the Hall 
resistivity (135) with n = 1. 

The nature of the transition between the compress¬ 
ible, weak coupling state at large interlayer separation 
and the incompressible, strong coupling state at d ~ £0 
is not completely understood (Eisenstein, 2014; Finck 
et al, 2010). Numerical evidence (Burkov et al, 2002; 
Schliemann et al, 2001) suggests a first order transition 
at r = 0, which contradicts the experimental observa¬ 
tion of gradual development of the quantized Hall drag 
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(Kellogg et ai, 2003, 2002b), see Fig. 33. Stern and 
Halperin (2002) suggested a phenomenological descrip¬ 
tion of the drag resistivity in the transition region. Pos¬ 
tulating that in the transition region the system is split 
into regions of the strong-coupling (1,1,1) phase and re¬ 
gions of the weak-coupling compressible phase, they de¬ 
scribe the transition as the point where the fraction / 
of the sample occupied by the (1,1,1) phase reaches the 
percolation threshold fc = 1/2. 

In a system of identical layers, the linear response the¬ 
ory can be formulated in terms of symmetric and anti¬ 
symmetric states. Denoting the 2x2 resistivity matri¬ 
ces corresponding to symmetric and antisymmetric cur¬ 
rents by p® and p“, one finds the drag resistivity as 

pD = (pa _ pS)/2. 

In the weak-coupling phase at £o, the drag resis¬ 
tivity is very small, <C Neglecting p^, one may 

approximate the resistivities as (in units of h/e^) 


p“(d»4) = p"(d»4) 



(136) 


where e = l/{kp£tr) ^ 1 (within the composite fermion 
model), kp = Ann is the Fermi wave vector, n is the elec¬ 
tronic density, and £tr is the transport mean free path. 
For T < IK, the experimentally measured values of p^ 
are almost two orders of magnitude less than e. 

The strong-coupling (1,1,1) phase exhibits features of 
the quantum Hall state for the symmetric currents 

Po(rf^4)= o) ’ 


while for the antisymmetric currents it is a superfluid 
(Stern and Halperin, 2002) pg(c? < £o) = 0. 

Analyzing the system close to the transition as a com¬ 
posite system comprising regions of both phases. Stern 
and Halperin (2002) have found a phenomenological ex¬ 
pression for the drag resistivity 


P 8e/(l-/)(!-2/) 

e2 -h 4(1 - 2/)2 


(138) 


As / increases from zero, this drag resistivity grows from 
zero [or rather, the very small value in the compressible 
state that is neglected in Eq. (138)] reaching a maximum 
at /* « 1/2 — e/4 (for small e) and again vanishing at 
the percolation threshold, in qualitative agreement with 
the non-monotonic drag observed in Kellogg et al. (2003, 
2002 b). 

Furthermore, using the semicircle law (Dykhne and 
Ruzin, 1994), it can be shown that to the lowest order 
in e the drag resistivities satisfy the relation (Stern and 
Halperin, 2002) 

(pf, + 1/2)2+ (pfj2^ 1/4^ (139) 




FIG. 33 (Color online) Left: Hall (open dots) and longitudi¬ 
nal (closed dots) drag resistance a,t vt ~ 1 and T = 50mK as 
a function of the ratio d/£o. The two lower panels show the 
temperature dependence of the location and the half-width 
of the peak in The lines are guides for the eye. [Re¬ 

produced from Kellogg et al. (2003).] Right: Drag resistivity 
and Hall drag resistivity in units of hje^ for different temper¬ 
atures. The end points represent the strong- (p]/^ = — 1) and 
weak-coupling {p^y = 0) regimes. The dashed line represents 
Eq. (139). [Reproduced from Tutuc et al. (2009).] 


yielding vanishing Hall drag for the compressible state 
(/ = 0) and the quantized value (135) with n = — 1 for 
the (1,1,1) state at / ^ 1/2. Between the two extremes 
the negative p^y varies monotonously. The apparent 
discrepancy in the sign of p^y obtained by Stern and 
Halperin (2002) and Kim et al. (2001) seem to stem 
from the alternative definition of drag resistivities. Sim¬ 
ilar predictions for transport coefficients, in particular 
Eq. (139), but without the explicit phase separation were 
obtained by Simon et al. (2003). An alternative model 
invoking the coexistence of the two phases was suggested 
by Spivak and Kivelson (2005). 

The semicircle relation (139) was experimentally tested 
in Tutuc et al. (2009), see Fig. 33. Instead of comparing 
a number of double-well devices with different interlayer 
separations (Kellogg et ai, 2003), Tutuc et al. (2009) var¬ 
ied the electron density and observed the transition be¬ 
tween the strong-coupling state ai vp = ^ and the weakly 
coupled state at I'p = 2. The data at intermediate tem¬ 
peratures T « 300K are in a good quantitative agreement 
with the theory. At the same time, Eq. (139) is only ap¬ 
proximate and is expected to hold if the drag resistivity 
is much larger than the symmetric bilayer resistivity at 
all fillings. Drag resistivity in the weak-coupling state is 
also neglected. Given these approximations, the agree¬ 
ment between the data and the phenomenological theory 
of Stern and Halperin (2002) is satisfactory. 


B. Interlayer exciton formation 

Further experiments revealed the most intriguing fea¬ 
ture of the strong-coupling quantum Hall state at vp = 1-. 
the presence of the exciton condensate capable of neutral 
superfluid transport (Eisenstein, 2014; Eisenstein and 
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MacDonald, 2004). Originally envisioned for optically 
generated excitons in bulk semiconductors (Blatt et ai, 
1962; Keldysh and Kopaev, 1964; Keldysh and Kozlov, 
1968; Moskalenko, 1962), the phenomenon has been pre¬ 
dicted also for indirect excitons in double-layer systems 
(Lozovik and Yudson, 1976; Shevchenko, 1976). 

The quantized Hall effect along with the vanishing lon¬ 
gitudinal resistivity at vt = 1 indicate a gapped spec¬ 
trum of charged excitations. In these measurements 
(Eisenstein, 2014; Eisenstein et aL, 1992; Suen et ai, 
1992), electrical currents in the two layers flow in the 
same direction. In contrast, the condensate couples to 
antiparallel or counterflowing currents (Kellogg et al., 
2004; Tutuc et al., 2004; Wiersma et ai, 2004) and man¬ 
ifests itself through vanishing Hall voltage. The simplest 
explanation for this observation is based on charge neu¬ 
trality of excitons: as neutral objects, excitons do not 
experience the Lorentz force and hence no Hall voltage 
develops when equal, counter-propagating currents are 
flowing through the two layers. 

Another spectacular manifestation of the exciton con¬ 
densate is the Josephson-like tunneling anomaly (Finck 
et al, 2008; Spielman et al.^ 2000; Tiemann et al, 2008a; 
Wiersma et al, 2006, 2007; Yoon et al, 2010) that the¬ 
oretically was predicted in Park and Das Sarma (2006) 
and Wen and Zee (1992) and later discussed in Dolcini 
et al (2010). 

Finally, the latest experiments revealing the existence 
of the exciton condensate utilized the multiple connected 
Corbino geometry (Finck et al, 2011; Nandi et al, 2012; 
Tiemann et al, 2008a,b). For a theoretical discussion 
of the superfluid flow in the Corbino geometry see Su 
and MacDonald (2008). The advantage of the Corbino 
samples is that they support the exciton flow through the 
bulk (in contrast to the Hall bar samples where transport 
is dominated by the edges). 

Coulomb drag has played an important role in discov¬ 
ering the interlayer correlated state (Eisenstein, 2014). 
Quantized Hall drag measured in the simply connected 
square geometry (Kellogg et al, 2002b) was one of 
the first indications of anomalous in-plane transport in 
double-layer systems at i/t = 1. Remarkably, the quan¬ 
tized Hall voltage has been found to be the same in both 
layers. At first glance, this contradicts the boundary con¬ 
ditions of the drag measurement: drag experiments in¬ 
volve passing current through one of the layers and mea¬ 
suring the induced voltage in the other, where no current 
is allowed to flow. The absence of the current seems to 
yield the absence of the Lorentz force and hence lead to 
the standard conclusion that no Hall voltage should be 
induced in the passive layer, see Sec. H.G. However, this 
argument does not take into account collective effects. In 
the presence of the condensate, the driving current can 
be decomposed into the symmetric and antisymmetric 
parts (Stern and Halperin, 2002). While the symmetric 
current carries the electric charge, the antisymmetric - 



FIG. 34 (Color online) Corbino Coulomb drag. Solid lines 
show the drive and drag currents. The measurement was per¬ 
formed at vt = 1, T — 17mK, and d/£o = 1-5. Dashed lines 
represent the results of simulations incorporating estimated 
series resistances and measured Corbino conductivity. The 
inset shows the measurement schematic. The resistances Ri 
comprise both external circuit resistors and the resistances in¬ 
trinsic to the device. [Reprinted by permission from Macmil¬ 
lan Publishers Ltd: Nature, Nandi et al (2012).] 

or counterpropagating - current is equivalent to the con¬ 
densate flow. In the passive layer the two currents cancel 
each other thus satisfying the boundary condition. At 
the same time, it is the symmetric, charge-carrying cur¬ 
rent that can couple to the magnetic field. This current 
is shared between the layers, yielding the identical quan¬ 
tized Hall voltage across both layers. 

Similar arguments lead to the expectation of “perfect” 
longitudinal drag (Su and MacDonald, 2008): the sym¬ 
metric current shared between the layers should be re¬ 
sponsible not only for the Hall, but also for the lon¬ 
gitudinal voltage in the passive layer. This prediction 
was tested in a dedicated experiment by Nandi et al 
(2012) using Corbino samples. Deviating from the stan¬ 
dard setup, Nandi et al (2012) have closed the electric 
circuit in the passive layer and measured the induced cur¬ 
rent, rather than the voltage. In this case, perfect drag 
means that the induced current should be same in mag¬ 
nitude as the driving current passed through the active 
layer while flowing in the opposite direction. This is ex¬ 
actly what has been observed by Nandi et al (2012), at 
least for small driving currents, see Fig. 34. 

The above arguments neglect the impact of disorder 
that might affect the presumed dissipationless excitonic 
transport (Su and MacDonald, 2008) across the bulk of 
the device (Fertig and Murthy, 2005; Fil and Shevchenko, 
2007; Huse, 2005; Lee et al, 2011). Assuming a phe¬ 
nomenological resistance Rg of the excitonic system, one 
still finds [neglecting the Corbino conductance (Nandi 
et al, 2012)] perfect drag h = h = V/iRi + R 2 + Rs), 
where Ri represent the net resistances in series with the 
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FIG. 35 (Color online) Coulomb drag in a grapliene-2DEG 
vertical heterostructure. Left: measured drag resistivity. The 
dashed line represents the best fit for the standard temper¬ 
ature dependence Rd = aT^, a = (5.8 ± 0.3) x 
Middle: a fit of the low-T upturn based on Eq. (140). The 
critical temperature found from the fit is Tc ~ 10 — lOOmK. 
Right: the low-T upturn in a bilayer graphene-2DEG het¬ 
erostructure. The fit based on Eq. (140) yields ~ lOOmK. 
[Reprinted by permission from Macmillan Publishers Ltd: 
Nature Communications, Gamucci et al. (2014).] 

Corbino sample, see the inset in Fig. 34. As the mag¬ 
nitude of Ri -I- i?2 is expected to always exceed 2h/e^ 
(Pesin and MacDonald, 2011; Su and MacDonald, 2008), 
the ability of the experiment to detect small values of 
Rs is limited. The issue of dissipation in the excitonic 
system might be clarified by future multi-terminal mea¬ 
surements. 

So far we have discussed experiments on the exciton 
physics in double-layer systems comprising similar elec¬ 
tronic layers (Eisenstein and MacDonald, 2004). It is 
also possible to create devices with oppositely doped lay¬ 
ers, the so-called electron-hole bilayers (Das Gupta et ai, 
2011; Keogh et ai, 2005). Coulomb drag measurements 
in these systems (Croxall et al., 2008; Seamons et ai, 
2009) do not provide a direct evidence of interlayer co¬ 
herence, but nevertheless demonstrate an upturn in pu as 
the temperature is lowered below IK. The upturn is seen 
only in devices with smaller (20nm) interlayer separation 
suggesting exciton formation. 

A microscopic theory of Coulomb drag in proximity 
to a phase transition was suggested by Hu (2000b) and 
Mink et al. (2012, 2013). As the system approaches the 
transition temperature Tc from above, the drag resistivity 
was found to exhibit a logarithmic divergence 

PD=Po + A\n[Tc/iT-Tc)], (140) 

where po and A are two fitting parameters (Gamucci 
et al., 2014). While qualitatively resembling the upturn 
observed in electron-hole bilayers (Croxall et al., 2008; 
Seamons et al., 2009), the theory accounts neither for a 
subsequent downturn at the lowest temperatures, nor the 
apparent violation of Onsager reciprocity (Croxall et al., 

2008) [although the latter might be related to heating 
effects (Seamons et ai, 2009)]. The theory also does not 
make falsifiable predictions regarding the dependence of 
Pd on carrier densities in the two layers (Morath et al., 

2009) [at higher temperatures, where the data show the 


standard dependence, the density dependence oi po 
is stronger than expected on the basis of the Fermi-liquid 
many-body calculations (Hwang and Das Sarma, 2008b)]. 

The logarithmic temperature dependence (140) fits 
well with the upturn in the drag resistivity observed in 
Gamucci et al. (2014) in hybrid devices comprising ei¬ 
ther a monolayer or bilayer graphene sheet and a GaAs 
quantum well, see Fig. 35. In fact, the search for exciton 
physics was one of the main motivations for experimental 
studies of Coulomb drag in double-layer graphene-based 
structures (Gorbachev et al., 2012; Kim et ai, 2011). 

Exciton condensation in graphene has attracted con¬ 
siderable theoretical attention (Abergel et al., 2013; 
Aleiner et al., 2007; Efimkin and Lozovik, 2011; Fil and 
Kravchenko, 2009; Kharitonov and Efetov, 2008, 2010; 
Lozovik et ai, 2012; Lozovik and Sokolik, 2008; Min 
et al., 2008; Pikalov and Fil, 2012; Sodemann et al., 2012; 
Suprunenko et al., 2012; Zhang and Joglekar, 2008). Sev¬ 
eral contradicting values of the transition temperature in 
double-layer graphene systems have been reported. The 
initial estimate (Mink et al., 2012; Zhang and Joglekar, 
2008) of Tc close to room temperature appeared to be 
too optimistic. Screening effects (Kharitonov and Efe¬ 
tov, 2008, 2010) were shown to lead to extremely low 
values under Imk [Tc ~ Ep). More recent investiga¬ 
tions involving detail analysis of screened Coulomb inter¬ 
action (Abergel et ai, 2013; Lozovik et al., 2012; Sode¬ 
mann et al., 2012), multiband pairing (Lozovik et al., 
2012; Mink et al., 2012), and pairing with nonzero mo¬ 
mentum (Efimkin and Lozovik, 2011) suggest somewhat 
higher values of Tc, making the transition experimentally 
accessible. 

High-temperature coherence and superfluidity has also 
been suggested in thin films of topological insulators 
(Efimkin et al., 2012; Mink et al., 2012, 2013; Seradjeh 
et al., 2009). 

The effect of exciton condensation on Coulomb drag 
has been investigated in graphene numerically by Zhang 
and Jin (2013) and in topological insulator films analyt¬ 
ically in Efimkin and Lozovik (2013). The latter work 
focused on the drag effect at temperatures exceeding Tc, 
where the pairing fluctuations are expected to play an 
important role. In addition to the Maki-Thompson-type 
contribution (Hu, 2000b; Mink et al, 2012, 2013) to the 
drag resistivity, Efimkin and Lozovik (2013) have ana¬ 
lyzed the Aslamazov-Larkin-type contribution and found 

Sp^^ (X [\n{T/Tc)]-\ (141) 

Far away from the transition, the result (141) decays log¬ 
arithmically, similarly to Eq. (140), but close to the tran¬ 
sition exhibits a stronger divergence 6pp^ oc (T — Tc)~^. 






51 


VIII. OPEN QUESTIONS AND PERSPECTIVES 

The physics of the Coulomb drag in double-layer sys¬ 
tems is well understood if both layers are in the Fermi- 
liquid state (Flensberg et al, 1995; Kamenev and Oreg, 
1995). The current in the passive layer is created by ex¬ 
citing electron-hole pairs (each pair consisting of an oc¬ 
cupied state above the Fermi surface and an empty state 
below) in a state characterized by finite momentum. The 
momentum comes from the electron-hole excitations in 
the active layer created by the driving current. The mo¬ 
mentum transfer is due to the interlayer Coulomb inter¬ 
action. Therefore it follows from the usual phase-space 
considerations that the drag coefficient is proportional 
to the square of the temperature po oc T^. Remark¬ 
ably, this simple argument is sufficient to describe the 
observed low temperature dependence of pn- Deviations 
from the quadratic dependence at higher temperatures 
are primarily due to the effect of phonons and plasmons 
(Rojo, 1999). 

The universality of the Landau Fermi-liquid theory 
(Altshuler and Aronov, 1985; Lifshitz and Pitaevskii, 
1981) can be traced to the linearization of the quasi¬ 
particle spectrum. Within this approximation all details 
of the microscopic structure of the system are contained 
in a limited number of parameters, such as the Fermi 
velocity and density of states (DoS) at the Fermi level. 
Many observable quantities (e.g. the electronic contribu¬ 
tion to the specific heat, spin susceptibility, period of the 
De Haas-van Alphen oscillations, etc.) can be expressed 
in terms of these parameters and thus exhibit the “uni¬ 
versal” behavior (as a function of temperature or external 
fields). Same arguments can be applied to elementary ex¬ 
citations in strongly doped graphene (/i ^ T), where the 
Fermi-liquid theory is expected to be applicable. 

Coulomb drag belongs to a different class of observ¬ 
ables. In conventional semiconductor devices, it reflects 
the degree of electron-hole asymmetry in the system 
vanishing in the approximation of linearized spectrum 
(Kamenev and Oreg, 1995). The drag coefficient is deter¬ 
mined by the subleading contribution taking into account 
the curvature of the quasi-particle spectrum. Indeed, in 
the passive layer the momentum is transferred equally to 
electrons and holes so that the resulting state can carry 
current only in the case of electron-hole asymmetry. Like¬ 
wise, this asymmetry is necessary for the current-carrying 
state in the active layer to be characterized by nonzero 
total momentum. The electron-hole asymmetry mani¬ 
fests itself (Narozhny et al, 2001; von Oppen et ai, 2001) 
in the energy (or chemical potential) dependence of such 
quantities as the density of states, single-layer conductiv¬ 
ity, and diffusion coefficient. Within the Fermi-liquid the¬ 
ory (Kamenev and Oreg, 1995), the asymmetry is weak, 
dui/dpi « (Jijpii leading to the drag effect, that is much 
weaker than the single-layer conductivity. 

Coulomb drag in non-Fermi-liquid systems is much 


more interesting. In particular, it has been used to study 
novel strongly-correlated, many-body states in double 
quantum wells (Eisenstein, 2014), graphene (Gamucci 
et ai, 2014), quantum wires (Laroche et ai, 2014), and 
optical cavities (Berman et al., 2014), where practical 
applications in optical switches have been suggested. In 
these systems, drag measurements have proved to be an 
invaluable tool to study the microscopic structure of com¬ 
plex, interacting many-body systems. 

At the same time, our understanding of many of these 
systems is incomplete. In contrast to the Fermi-liquid 
theory, many aspects of the strongly-correlated many- 
body states lack a detailed theoretical description. Con¬ 
sequently, their transport properties, including Coulomb 
drag, can be evaluated only with the help of heuristic 
or phenomenological models. One can only hope that a 
proper microscopic theory of these effects will eventually 
be developed. 

This brings us to the list of unresolved questions re¬ 
lated to the theory reviewed in this paper and possible 
direction of the held in the near future. 

(i) At low enough temperatures and especially in 
strong magnetic helds, double-layer systems may host 
excitonic condensates (Eisenstein, 2014). In monolayer 
graphene, such condensation is also possible, but for rea¬ 
sonably weak interactions the condensation temperature 
appears to be rather low (Aleiner et ai, 2007; Kharitonov 
and Efetov, 2008; Mink et al, 2012). Nevertheless, a pos¬ 
sibility of interlayer correlated states in graphene-based 
systems [and possibly in hybrid devices involving other 
materials, (Geim and Grigorieva, 2013)] is rather exciting 
and certainly requires theoretical attention. 

(ii) The hydrodynamic approach of Sec. IV.B should 
be extended to include thermoelectric effects in graphene- 
based double-structures as well as in monolayer and bi¬ 
layer graphene. As pointed out in Foster and Aleiner 
(2009) and Narozhny et al. (2015), the quasiparticle im¬ 
balance in graphene may play a decisive role in ther¬ 
mal transport. Another promising direction may be 
opened by generalization of the macroscopic linear- 
response equations to a true, nonlinear hydrodynamics. 
The relation between the quantum kinetic equation of 
Zala et al. (2001) and the hydrodynamic approach [both 
in graphene (Narozhny et ai, 2015) and in 2DEG (An¬ 
dreev et al., 2011; Apostolov et al., 2014)] is also of cer¬ 
tain theoretical interest. 

(iii) Dirac fermions can be found as low-energy exci¬ 
tations not only in graphene, but also in topological in¬ 
sulators (Bernevig and Hughes, 2013; Shen, 2013). An 
extension of the present theory of Coulomb drag to var¬ 
ious possible system configurations involving topological 
insulators and/or hybrid devices involving topological in¬ 
sulators, graphene, etc. appears to be very promising 
(Mink et al., 2012). 

(iv) Novel aspects of Luttinger liquid physics and role 
of equilibration processes on drag can be further explored 
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with the edge states of quantum Hall systems or topo¬ 
logical edge liquids of quantum spin Hall effect. Some 
theoretical predictions have already being made (Zyuzin 
and Fiete, 2010) and recent experimental advances (Al- 
timiras et al, 2010; Du et al, 2015; Konig et ai, 2013; 
Roth et ai, 2009) bring these exciting perspectives within 
reach. 

(v) Mesoscopic fluctuations of Coulomb drag in ballis¬ 
tic samples should be further analyzed on the basis of 
the microscopic theory. The theory should be further ex¬ 
tended to the cases of Dirac fermions in graphene and 
composite fermions at the half-filled Landau level. Ex¬ 
perimental work in this direction has been already initi¬ 
ated in Kim et al. (2011) and Price et al. (2010). 

(vi) The third-order drag effect (see Sec. IV.F) bears a 
certain resemblance to the well-known Altshuler-Aronov 
corrections to single-layer conductivity (Altshuler and 
Aronov, 1985; Zala et al, 2001). In Zala et al. (2001) 
it was shown that the dominant contribution to conduc¬ 
tivity at low (diffusive regime) and high (ballistic regime) 
temperatures technically comes from different diagrams 
describing conceptually similar, but at the same time dis¬ 
tinct interference processes. Similarly, we expect that the 
third-order drag contribution in ballistic regime might be 
governed by scattering processes which are distinct from 
those considered in Levchenko and Kamenev (2008a). 

We would like to close this review by pointing out the 
surprising richness of the Coulomb drag problem. The 
original suggestion of a way to observe interwell inter¬ 
actions in semiconductor heterostructures has developed 
into a vibrant field of research where technological ad¬ 
vances go hand in hand with theoretical developments. 
New experiments with novel materials keep being devised 
and stimulate new avenues for theoretical thinking. We 
should be expecting to see further intriguing discoveries 
being made related to frictional drag in the foreseeable 
future. 
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